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Response to a step force at FT. CMC 
uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.9 500 

IV-21a Depth vs. Time. Response to a step 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.8 1)0'J 
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IV-21b Pitch vs. Time. Response to a step 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.8 50'S 

IV-21c Stem Plane Angle vs. Time. Response 

to a step force at FT. CMC uses optimal 

design in combined controllers with 

B=800, C=10, E=l. Parameter X=0.8 50“^ 



IV-21d Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. 

Parameter X=0.8 5) O 

IV-22a Depth vs. Time, Response to a step 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.7 5 A j 

IV-22b Pitch vs. Time. Response to a step 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.7 ^^2. 

IV-22c Stem Plane Angle vs. Time. Response 

to a step force at FT. CMC uses optimal 

design in combined controllers with 

B=800, C=10, E=l. Parameter X=0. 7 o 5 

IV-22d Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. 



Parameter X=0.7 3 \ ^ 

IV-23a Depth vs. Time. Response to a step 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.05 

IV-23b Pitch vs. Time. Response to a step 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.05 3 I G 
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IV-23C 



IV-23d 



IV-24a 



IV -24b 



IV -24c 



IV-24d 



IV-25a 



IV- 25b 



IV-25C 



Stem Plane Angle vs. Time. Response 
to a step force at FT. CMC uses optimal 
design in combined controllers with 
B=800, C=10, E=l. Parameter X=0.05 

Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. 

Parameter X=0.05 

Depth vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 
C=0.001, E=0.001. Parameter X=0.9 

Pitch vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 
C»0.001, E=0.001. Parameter X=0.9 

Stem Plane Angle vs. Time. Response 
to a step force at FT. CMC uses optimal 
design in combined controllers with 
B=164, C=0.001, E=0.001. Parameter 
X=0.9 



Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, E=0.001. 
Parameter X=0.9 

Depth vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 
C=0.001, E=0.001. Parameter X=0.5 

Pitch vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 
C=0.001, E=0.001. Parameter X=0.5 

Stem Plane Angle vs. Time. Response 
to a step force at FT. CMC uses optimal 
design in combined controllers with 
B=164, C=0.001, E=0.001. Parameter 
X=0.5 - 



5^7 

8 

3 1 8 
8^0 

82 ! 

822 

883 

884 

828 
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IV-25d 



IV-26a 



IV-26b 



IV-26C 



IV-26d 



IV-27a 



IV-27b 



IV-27C 



Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, 

E=0.001. Parameter X=0.5 

Depth vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 

C=0.001, E=0.001. Parameter X=0.4 

Pitch vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 

C=0.001, E=0.001. Parameter X=0.4 2)28 

Stem Plane Angle vs. Time. Response 

to a step force at FT. CMC uses optimal 

design in combined controllers, with 

B=164, C=0.001, E=0.001. Parameter 

X=0.4 228 

Fairwater- Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=164, C=0.001, 



E=0.001. Parameter X=0. 4 2)50 

Depth vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 

C=0.001, E=0.001. Parameter X=0.3 

Pitch vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.'001, E=0.001. Parameter X=0.3 5o2 



Stem Plane Angle vs. Time. Response 

to a step force at FT. CMC uses optimal 

design in combined controllers, with 

B=164, C=0.001, E=0.001. Parameter 

X=0.3 - - 3ob 
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IV-27d 



IV-28a 



IV-28b 



IV-28C 



IV-28d 



IV-29a 



IV-29b 



IV-29c 



IV-29d 



Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=164, C=0.001, 

E=0.001. Parameter X=0.3 

Depth vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 

C=0.001, E=0.001. Parameter X=0.05 

Pitch vs. Time. Response to a step 
force at FT. CMC uses optimal design 
in combined controllers, witi B=164, 

C=0.001, E=0.001. Parameter X=0.05 

Stem Plane Angle vs. Time. Response 
to a step force at FT. CMC uses optimal 
design in combined controllers, with 
B=164, C=0.001, E=0.001. Parameter 
X=0.05 5^1 

Fairwater Plane Angle vs. Time. 

Response to a step force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=164, C=0.001, 



E=0.001. Parameter X=0.05 

Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.7 

Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.7 3^0 

Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses optimal 
design in combined controllers, with 
B=800, C=10, E=l. Parameter X=0.7 



Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. 

Parameter X=0.7 5^2. 



2k 



IV-30a Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.6 

IV-30b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.6 

IV-30c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses optimal 
desing in combined controllers with 
B=800, C=10, E=l. Parameter X=0.6 

IV-30d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. 

Parameter X=0.6 

IV-31a Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.55 5^7 

IV-3lb Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0.55 5 A B 

IV-31c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses optimal 
design in combined controllers with 
B=800, C=10, E=l. Parameter X=0.55 

IV-31d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. _ 

Parameter X=0.55 

IV-32a Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.5 
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IV-32b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.5 

IV-32c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses optimal 

design in combined controllers, with 

B=800, C=10, E=l. Parameter X=0.5 5 “DO 

IV-32d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=800, C=10, E=l. 



Parameter X=0.5 

IV-33a Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.45 DtiD 

IV-33b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.45 2)'D'S> 

IV-33c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses optimal 
design in combined controllers, with _ 

B=800, C=10, E=l. Parameter X=0.45 2)D7 



IV-33d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=800, C=10, E=l. 



Parameter X=0.45 2)D8 

IV-34a Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, _ 

C=10, E=l. Parameter X=0.3 

IV-34b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.3 3GO 
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IV-34c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses 
optimal design in combined controllers, 
with B=800, C=10, E=l. Parameter 
X=0.3 - 

IV-34d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=800, C=10, E=l. 



Parameter X=0.3 ^ G2. 

IV-35a Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.05 'hG'b 

IV-35b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers with B=800, 

C=10, E=l. Parameter X=0. 05 



IV-35c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses 
optimal design in combined controllers, 
with B=800, C=10, E=l. Parameter 
X=0.05 - 5G5 

IV-35d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=800, C=10, E=l. 



Parameter X=0.05 

IV-36a -Depth vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.005 567 

IV-36b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses optimal design 
in combined controllers, with B=800, 

C=10, E=l. Parameter X=0.005 56 S 
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IV-36c 



IV-36d 



IV-37a 



IV-37b 



IV-37C 



IV-37d 



IV-38a 



IV-38b 



Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers, 
with B=800, C=10, E=l. Parameter 

x=o.oo5 

Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=800, C=10, E=l. 

Parameter X=0.005 570 

Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 

C=0.001, E=0.001. Parameter X=0.5 57 ] 

•Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.001, E=0.001. Parameter X=0.5 57 2. 

Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers, 
with B=164, C=0.001, E=0.001. 

Parameter X=0.5 

Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, 

E=0.001. Parameter X=0.5 



Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.001, E=0.001. Parameter X=0.3 575 

Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.001, E=0.001. Parameter X=0.3 



577 



575 
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IV-38c 


Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers 
with B=164, C=0.001, E=0.001. 

Parameter X=0.3 ?)H 


IV-38d 


Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, 

E=0.001. Parameter X=0.3 1)1^ 


IV-39a 


Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 
C=0.001, E=0.001. Parameter X=0.2 


IV- 39b 


Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C-0.001, E=0.001. Parameter X=0.2 380 


IV-39C 


Stem plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers, 
with B=164, C=0.001, E=0.001. 

Parameter X=0.2 38^ 


IV-39d 


Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, 

E=0.001. Parameter X=0. 2 


IV -40a 


Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 

C=0.001, E=0.001. Parameter X=0.15 383 


IV -40b 


Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.001, E=0.001. Parameter X=0.15 3 88 
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IV -40c 


Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers, 
with B=164, C=0.001, E=0.001. 

Parameter X=0.15 ^85 


IV-40d 


Fairwater Plane Angle vs. Time 
Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers, with B=164, C=0.001, 

E=0.001. Parameter X=0.15 


IV-41a 


Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 
C=0.001, E=0.001. Parameter X=0.1 


IV -41b 


Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers, with B=164, 
C=0.001, E=0.001. Parameter X=0.1 


IV-41C 


Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers 
with B=164, C=0.001, E=0.001. 

Parameter X=0.1 ^88 


IV-41d 


Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, ^ ^ 

E=0.001. Parameter X=0.1 


IV-42a 


Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.001, E=0.001. Parameter X=0. 005 88] 


IV -42b 


Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with B=164, 

C=0.001, E=0.001. Parameter X=0.005 88 2. 
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IV-42c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses 
optimal design in combined controllers 
with B=164, C=0.001, E=0.001. 

Parameter X=0.005 

IV-42d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with B=164, C=0.001, 



E=0.001. Parameter X=0.005 5^^ 

IV-43a Depth vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 

loop CE. Parameter X=0.7 595 

IV -43b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 

loop CE. Parameter X=0.7 59)G 



IV-43c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses BPOC 
with B=800, C=10, E=1 and SOPC with the 
same closed loop CE. Parameter X=0.7 59 T 

IV-43d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses BPOC with B=800, C=10, E=1 and 
SOPC with the same closed loop CE. 



Parameter X=0.7 598 

IV-44a Depth vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, C=1 and SOPC with the same closed 
loop CE. Parameter X=0.5 399 

IV-44b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 

loop CE. Parameter X=0.5 i^OO 



IV-44c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses BPOC 
with B=800, C=10, E=1 and SOPC with the 
same closed loop CE. Parameter X=0.5 \ 
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IV-44d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses BPOC with B=800, C=10, E=1 and 
SOPC with the same closed loop CE. 

Parameter X=0.5 A02 

IV-45a Depth vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 
loop CE. Parameter X=0.4 40 5 

IV-45b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 
loop CE. Parameter X=0.4 

IV-45c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses BPOC 
with B=800, C=10, E=1 and SOPC with the 
same closed loop CE. Parameter X=0.4 

IV-45d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses BPOC with B=800, C=10, E=1 and 
SOPC with the same closed loop CE. 



Parameter X=0.4 ■4o6 

IV-46a Depth vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed ^ 

loop CE. Parameter X=0.3 ^01 

IV-46b Pitch vs. Time. Response to a pulse 

force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 

loop CE. Parameter X=0.3 40 8 



IV-46c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. CMC uses BPOC 
with B=800, C=10, E=1 and SOPC with the 
same closed loop CE. Parameter X=0.3 

IV-46d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses BPOC with B=800, C=10, E=1 and 
SOPC with the same closed loop CE. 

Parameter X=0.3 410 
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IV-47a 



IV-47b 



IV-47C 



IV-47d 



IV-48a 



IV-48b 



IV-48C 



IV-48d 



IV-49a 



Depth vs. Time. Response to a pulse 
force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 
loop CE. Parameter X=0.1 AU 

Pitch vs. Time. Response to a pulse 
force at FT. CMC uses BPOC with B=800, 

C=10, E=1 and SOPC with the same closed 
loop CE. Parameter X=0.1 A ^ 2. 



Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses BPOC 
with B=800, C=10, E=1 and SOPC with the 
same closed loop CE. Parameter X=0.1 Al^ 



Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses BPOC with B=800, C=10, E=1 and 
SOPC with the same closed loop CE. 

Parameter X=0.1 AJA 

Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with C=10, 

D=3000, E=l. Parameter X=0.4 A ^ !) 

Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with C=10, 

D=3000, E=l. Parameter X=0. 4 A^G 

Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers 
with C=10, D=3000, E=l. Parameter 

x=o .4 An 



Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with C=10, D=3000, E=l. 

Parameter X=0.4 •A B 

Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with C=10, 

D=3000, E=l. Parameter X=0.5 



An 



IV -49b 


Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controll'ers with C-10, 

D=3000, E=l. Parameter X=0. 5 4 5.0 


IV-49C 


Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers 
with C=10, D=3000, E=l. Parameter 

x=0. 5 


IV-49d 


Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. CMC 
uses optimal design in combined con- 
trollers with C=10, D=3000, E=l. 

Parameter X=0. 5 4 2.5, 


IV-50a 


Depth vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with C=10, 

D=3000, E=1.0 . Parameter X=0,6 445 


IV-50b 


Pitch vs. Time. Response to a pulse 
force at FT. CMC uses optimal design 
in combined controllers with C=10, 

D=3000, E=1.0 .Parameter X=o.6 


IV-50C 


Stem Plane Angle vs. Time. Response 
to a pulse force at FT. CMC uses 
optimal design in combined controllers 
with C=10, d=3000, E=1.0 .Parameter X=0.6" 


IV-50d 


Fairwater Plane Angle vs. Time. 

Response to a pulse force of FT. CMC 
uses optimal design in combined con- 
trollers with C=10, D«3000, E=1.0 .X=0.6- mSG 


V-la 


Root locus for SOPOC with B=800, 

C=10, E=1.M > 0 4 41 


V-lb 


Root locus for SOPOC with B=800, 

C=10, E=1.M<0 4 5? 


V-2a 


Depth vs. Time. Response to a pulse 
force at FT. SOPC equal to 0.3* SOPOC 
(B=800, C=10, E=l) 4 
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V-2b Pitch vs. Time. Response to a pulse 

force at FT. SOPC equal to 0.3-SOPOC 
(B=800, C=10, E=l) 

V-2c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. SOPC equal to 
0.3-SOPOC (B=800, C=10, E=l) 

V-3a Depth vs. Time. Response to a pulse 

force at FT. SOPC equal to 0.6’SOPOC 
(B=800, C=10, E=l) 

V-3b Pitch vs. Time. Response to a pulse 

force at FT. SOPC equal to 0.6*SOPOC 
(B=800, C=10, E=l) 

V-3c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. SOPC equal to 
0.6-SOPOC (B=800, C=10, E=l) 

V-4a Depth vs. Time. Response to a pulse 

force at FT. SOPC equal to 0.8-SOPOC 
(B=800, C=10, E=l) 

V-4b Pitch vs. Time. Response to a pulse 

force at FT. SOPC equal to 0.8*SOPOC 
(B=800, C=10, E=l) 

V-4c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. SOPC equal to 
0.8-SOPOC (B=800, C=10, E=l) 4^7 



V-5a Root locus for SOPC with CE same as 

for BPOC (B=800, C=10, E=1).M>0 -- 

V-5b Root locus for SOPC with CE same as 

■ for BPOC (B=800, C=10, E=1).M<0 -- 

V-6a Root locus for SOPOC with B=164, 

C=0.001, E=0.001.M>0 4 40 

V-6b Root locus for SOPOC with B=164, 

C=0.001, E=0.001.M<,0 

V-7a Root locus for SOPC with CE same as 

for BPOC (B=164, C=0.001, E=0.001). , ^ 

M>0 




V-7b Root locus for SOPC with CE same as 

for BPOC (B=164, C=0.001, E=0.001). 

M<.0 AA?) 

V-8a Root locus for SOPOC with D=3000, 

C=10, E=1 M> 0 MMH 

V-8b Root locus for SOPOC with D=3000, 

C=10, E=1 M40 

V-9a Root locus for SOPC with CE same as 

for BPOC (D=3000, C=10, E=l) M>0 

V“9b Root locus for SOPC with CE same as 

for BPOC (D=3000, C=10, E=l) M<0 

V-lOa Root locus for SOPC using 5si of 

BPOC (D=3000, C=10, E=1).Mt>0 

V-lOb Root locus for SOPC using Ssi of 

BPOC (D=3000, C=10, E=1).MicO 

V-lOc Root locus for SOPC using Ssi of 

BPOC (B=800, C=10, E=l).Mi >0 

V-lOd Root locus for SOPC using §51 of 

BPOC (B=800, C=10, E=l). Mi 40 

V-lOe Root locus for SOPC using SsL of 

BPOC (B=164, C=0.001, E=0.001).Mi >0 A 52 

V=10f Root locus for SOPC using 5si of 

BPOC (B=164, C=0.001, E=0.001).Mi <0 455 

V-lla Depth vs. Time. Response to a pulse 

force at FT. Both planes submarine 
with 6sco-rq = • f(k) . X=0.5 ^*0 A 

V-llb Pitch vs. Time. Response to a pulse 

force at FT. Both planes submarine 
with 6scoTr) =6si • f (k) . X=0.5 A55 

V-llc Stem Plane Angle vs. Time. Response 

to a pulse force at FT. Both planes 
submarine with Si,corq =osi* f(k) . 

X=0.5 A 5G 
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V-lld Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. 

Both planes submarine with 
6 scovn =<Ssi-f(k). X=0.5 

V-12a Depth vs. Time. Response to a pulse 

force of FT. Both planes submarine 
with 6 scott^ = 6 si* f(k) . X=-0.5 

V-12b Pitch vs. Time. Response to a pulse 

force of FT. Both planes submarine 
with 6500 x ^1 =osi-f(k). X=-0.5 A 50 

V-12c Stem Plane Angle vs. Time. Response 

to a pulse force of FT. Both planes 
submarine with^scoi^ = 6 si • f(k) . 

X=-0.5 A60 

V-12d Fairwater Plane Angle vs. Time. 

Response to a pulse force of FT. 

Both planes submarine with 

6 scOTt-] = 6 si.f(k). X=-0.5 

V-13a Depth vs. Time. Response to a pulse 

force at FT. Both planes submarine 
with 6 scoxx-) = 65 ) • f (k) . X=-o-75 0 € ^ 

V-13b Pitch vs. Time. Response to a pulse 

force at FT. Both planes submarine 
with 6 scorn =6sl* f(k) . X=-0.75 

V-13c Stem Plane Angle vs. Time. Response 

to a pulse force at FT. Both planes 
submarine witbSscoxn = 6 si- f(k) . , ^ , 

X=-0.75 

V-13d Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. 

Both planes submarine with 

&corn X=-0.75 HS'b 

V-14a Depth vs. Time. Response to a pulse 

force at FT. Both planes submarine 
withSscoxT] = 6 si-f(k). X=- 1.0 A*oG 
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V-I4b 

V-I4c 

V-I4d 

V-15a 

V-15b 

V-15c 

V-15d 

B-1 

C-la 

C-lb 

C-lc 

C-2a 



Pitch vs. Time. Response to a pulse 
force at FT. Both planes submarine 
with6sco'riq =6$)L • f (k) . X=-l.O 

Stem Plane Angle vs. Time. Response 
to a pulse force at FT. Both planes 
submarine with Sscorn =6&i • f(k) . 

X=-1.0 - 

Fairwater Plane Angle vs. Time. 

Response to a pulse force at FT. 

Both planes submarine with 

Sscorri = 6s\ . f(k) . x=-l.0 

Depth vs. Time. CMC 

Pitch vs. Time. CMC 

Stem Plane Angle vs. Time. CMC 

Fairwater Plane Angle vs. Time. CMC 

Relation Bemeen Body Axes and Fixed 
Axes 

Depth vs. Time. Stem planes only 
submarine "not in trim" with feedback 
controller (SOA=-0.5, SOB=-7.03, SOC= 
32.5, SOD=360) and zero pitch and depth 
ordered 

Pitch vs. Time. Stem planes only 
submarine "not in trim" with feedback 
controller (SOA=-0.5, SOB=-7.03, SOC= 
32.5, SOD=360) and zero pitch and depth 
ordered 

Stem Plane Angle vs. Time. Stem planes 
only submarine "not in trim" with feed- 
back controller (SOA=-0.5, SOB=-7.03, 
SOC=32.5, SOD=360) and zero pitch and 
depth ordered 

Depth vs. Time. Stem planes only 
submarine "not in trim" with feedback 
controller (SOA=-0 . 0436 , SOB=-3.49, 
SOC=0.0523, SOD=360) and zero pitch 
and depth ordered 

3^ 



AST 
A €3 

ASA 

Mio 

m\ 

uis 



C-2b 



C-2c 



C-3a 



C-3b 



C-3c 



C-4a 



C-4b 



C-4c 



Pitch vs. Time. Stem planes only- 
submarine "not in trim" with feedback 
controller (SOA=-0 . 0436 , SOB=-3.49, 

SOC=0.0523, SOD=360) and zero pitch 

and depth ordered A'B' 

Stem Plane Angle vs. Time. Stem 
planes only submarine "not in trim" 
with feedback controller (SOA=-0 . 0436 , 

SOB=-3.49, SOC=0.0523, SOD=360) and 

zero pitch and depth ordered A '9 '^ 

Depth vs. Time. Stem planes only 
submarine "not in trim" with feedback 
controller (SOA=-0. 0656, S0B=-1.47, 

SOC=26.2, S0D=184.1) and zero pitch 

and depth ordered A83 

Pitch vs. Time. Stem planes only 
submarine "not in trim" with feedback 
controller (SOA=-0 . 0656 , SOB=-1.47, 

SOC=26.2, S0D=184.1) and zero pitch 
and depth ordered 

Stem Plane Angle vs. Time. Stem 
planes only submarine "not in trim" 
with feedback controller (SOA=-0 . 0656 , 

S0B=-1.47, SOC=26.2, S0D=184.1) and 
zero pitch and depth ordered 

Depth vs. Time. Stem planes only 
submarine "in trim" with feedback 
controller (SOA=-0.5, SOB=-7.03, S0C= 

32.5, SOD=360) and zero pitch and 
depth ordered 

Pitch vs. Time. Stem planes only 
submarine "in trim" with feedback 
controller (SOA=-0.5, SOB=-7.03, 

SOC=32.5, SOD=360) and zero pitch 

and depth ordered 

Stem Plane Angle vs. Time. Stem 
planes only submarine "in trim" with 
feedback controller (SOA=-0.5, SOB= 

-7.03, SOC=32.5, SOD=360) and zero 

pitch and depth ordered ® 
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C-5a 



C-5b 



C-5c 



C-6a 



C-6b 



C-6c 



C-7a 



C-7b 



Depth vs. Time. Stem planes only 
submarine "in trim" with feedback 
controller (SOA=-0 . 0436 , SOB=-3.49, 
SOC=0.0523, SOD=360) and zero pitch 
and depth ordered ^ 

Pitch vs. Time. Stem planes only 
submarine "in trim" with feedback 
controller (SOA=-0.0436, SOB=-3.49, 
SOC=0.0523, SOD=360) and zero pitch 
and depth ordered 

Stem Plane Angle vs. Time. Stem 
planes only submarine "in trim" with 
feedback controller (SOA=-0 . 0436 , 
SOB=-3.49, SOC=0.0523, SOD=360) and 
zero pitch and depth ordered 

Depth vs. Time. Stem planes only 
submarine "in trim" with feedback 
controller (SOA=-0.0656, SOB=-1.47, 
SOC=26.2, SOD=184.1) and zero pitch 
and depth ordered 

Pitch vs. Time. Stem planes only 
submarine "in trim" with feedback 
controller (SOA=-0 . 0656 , S0B=-1.47, 
SOC=26.2, S0D=184.1) and zero pitch 
and depth ordered 

Stem Plane Angle vs. Time. Stem 
planes only submarine "in trim" with 
feedback controller (SOA=-0.0656, 
S0B=-1.47, SOC=26.2, S0D=184.1) and 
zero pitch and depth ordered 

Depth vs. Time. Stem planes only 
submarine "not in trim" with SOPOC 
(B=164, C=0.001, E=0.001) and zero 
pitch and depth ordered 

Pitch vs. Time. Stem planes only 
submarine "not in trim" with SOPOC 
(B=164, C=0.001, E=0.001) and zero 
pitch and depth ordered 






Ago 



A 



A97 







A9^ 



A) G 
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C-7c 

C-8a 

C-8b 

C-8c 

C-9a 

C-9b 

C-9c 

C-9d 

C-lOa 



Stem Plane Angle vs. Time. Stem 
planes only submarine "not in trim" 
with SOPOC (B=164, C=0.001, E=0.001) 
and zero pitch and depth ordered 

Depth vs. Time. Stern planes only 
submarine "in trim" with SOPOC 
(B=164, C=0.001. E=0.001) and zero 
pitch and depth ordered 

Pitch vs. Time. Stem planes only 
submarine "in trim" with SOPOC 
(B=164, C=0.001, E=0.001) and zero 
pitch and depth ordered 

Stem Plane Angle vs. Time. Stem 
planes only submarine "in trim" 
with SOPOC (B=164, C=0.001, E=0.001) 
and zero pitch and depth ordered 

Sinusoidal force at AU 

Depth vs. Time. Response to sinusoidal 
force at AU. Stern planes only sub- 
marine "in trim" with SOPOC (B=164, 
C=0.001, E=0.001) and zero pitch and 
depth ordered 

Pitch vs. Time. Response to sinusoidal 
force at AU. Stem planes only sub- 
marine "in trim" v>7ith SOPOC (B=164, 
C=0.001, E=0.001) and zero pitch and 
depth ordered 

Stem Plane Angle vs. Time. Response 
to sinusoidal force at AU. Stem 
planes only submarine "in trim" with 
SOPOC (B=164, C=0.001, E=0.001) and 
zero pitch and depth ordered 

Depth vs. Time. Response to sinusoidal 
force at AU. Stem planes only sub- 
marine "in trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD= 

360) and zero pitch and depth ordered 



“500 

5oi 

go 
go 
go 5 
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C-lOb 



C-lOc 



C-lla 



C-llb 



C-llc 



C-12a 



C-12b 



Pitch vs. Time. Response to sinusoidal 
force at AU. Stem planes only sub- 
marine "in trim" with feedback con- 
troller (SOA=-0.5, SOB=-7.03, SOC=32.5, 

SOD=360) and zero pitch and depth 

ordered G 

Stem Plane Angle vs. Time. Response 
to sinusoidal force at AU. Stem 
planes only submarine "in trim" with 
feedback controller (SOA=-0.5, SOB= 

-7.03, SOC=32.5, SOD=360) and zero 

pitch and depth ordered ‘7)07 

Depth vs. Time. Response to sinusoidal 
force at AU. Stern planes only sub- 
marine "in trim" with feedback con- 
troller (SOA=-0.0436, SOB=-3.49, SOC= 

0.0523, SOD=360) and zero pitch and 
depth ordered 

Pitch vs. Time. Response to sinusoidal 
force at AU. Stem planes only sub- 
marine "in trim" with feedback con- 
troller (SOA=-0.0436, S0B=-3.49, S0C= 

0.0523, SOD=360) and zero pitch and 
depth ordered 

Stern Plane Angle vs. Time. Response 
to sinusoidal force at AU. Stem 
planes only submarine "in trim" with 
feedback controller (SOA=-0 .0436, 

S0B=-3.49, SOC=0.0523, SOD=360) and 

zero pitch and depth ordered 7) lO 

Depth vs. Time. Response to sinusoidal 
force at AU. Stern planes only sub- 
marine "in trim" with feedback con- 
troller (SOA=-0.0656, SOB=-1.47, S0C= 

26.2, S0D=184.1) and zero pitch and __ 

depth ordered 7> ^ ^ 

Pitch vs. Time. Response to sinusoidal 
force at AU. Stem planes only sub- 
marine "in trim" with feedback con- 
troller (SOA=-0.0656, S0B=-1.47, S0C= 

26.2, S0C=184.1) and zero pitch and 

depth ordered S\2. 



42 



C-12c Stem Plane Angle vs. Time. Response 

to sinusoidal force at AU. Stem 
planes only submarine "in trim" with 
feedback controller (SOA=-0 . 0656 , 

SOB=-1.47, SOC=26.2, S0C=184.1) and 
zero pitch and depth ordered 

C-13a Depth vs. Time. Response to sinusoidal 

input at AU. Stem planes only sub- 
marine "in trim" with SOPOC (B=800, 

C=10, E=l) and zero pitch and depth 

ordered *5 ^ ^ 

C-13b Pitch vs. Time, Response to sinusoidal 

input at AU. Stem planes only sub- 
marine "in trim" with SOPOC (B=800, 

C=10, E=l) and zero pitch and depth _ _ 

ordered t) ' ‘O 

C-13c Stem Plane Angle vs. Time. Response 

to sinusoidal input at AU. Stem 
planes only submarine "in trim" with 
SOPOC (B=800, C=10, E=l) and zero 



pitch and depth ordered “5 ' ^ 

D-1 Ballast vs. Speed 

D-2 Xq vs. Speed T) ^ S 

E-1 Kvs.S^co^ and X=3.2 - t) '' 9 

E-2 K vs.§i^Coni and X=1.2 

E-3 K vs.^lcc-u and X=0 . 7 

E-4 K vs.S^can andX=0.5 

E-5 K vs.6Ko"U and X=0.3 ‘O 2^0 

E-6 K vs.6^corr) and X=0.05 9-^ 
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INTRODUCTION 



Operational experience has shown that there is a definite 
need for a new, or at least improved, method for control of 
a near surface submarine. 

"Near surface" generally refers to a submarine which is 
close enough to the surface to be affected by surface waves. 

The purpose of this paper was to devise a depth and pitch 
controller that could improve the sea state capability and 
attain more accurate depth control. In order to decide on 
the proper control strategy, gaining some knowledge of the 
nature of the near surface problem was required. 

In Chapter I, some of the reasons for a submarine to be 
near STirface are listed. This is followed by a section pre- 
senting elements of wave theory. The next section then con- 
siders the effects of the surface excitation, on the submarine. 
Due to these excitations the submarine operatina in the near 
surface environment may be deflected far enough from its 
ordered trajectory to cause position saturation of the for- 
ward planes. The resulting situation is essentially a stern 
planes only control system. This configuration may be unstable. 

To prevent serious problems in this area, the controller 
must be able to perform the following functions : 

1. Identification of plant dynamics. 

2. Evaluation of present position in relation to an opti- 
mum defined and action to drive the system toward the optimum 
by modification of controller parameters. 
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Closed loop optimization is used for plants which exhibit 
substantial changes in their dynamic characteristics and where 
random disturbances may ha.ve considerable influence on the 
plant state. 

Aiming in that direction, two control systems were 
designed. The one is based on stern planes only control and 
is presented in Chapter II. The other is based on having both 
planes available and is presented in Chapter III, in three 
slightly different forms. 

The two separate controllers were designed using optimal 
control theory and their combination led to a configuration 
capable of improving system performance. The submarine wa.s 
subjected to step and pulse forces and the controller tested 
by simulation on the computer. In Chapter IV, the new scheme 
of controller is introduced and the results of the test are 
presented and analyzed. In Chapter V, further investigation 
of the combined scheme was attempted in an effort to better 
understand its function and test the generality of its 
parameters-. 

Conclusions by the author conlude the main body of the 
paper. However , seven Appendices are included. Appendix A 
describes the linearization process of equations of motion. 

In Appendix 3 a derivation of equations of motion from first 
principles is repeated from Ref. 9- A method for relating 
classical feedback to optimal theory is presented in Appendix 
C and a simplified trim analysis in Appendix D. In Appendix 
E, computer programs are described , setting limits on values 
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of parajneters in the combined controller of Chapter IV. La.stly, 
in Appendix F the equations of motion for a submarine in six 
degrees of freedom are repeated from Ref.l and in Appendix G 
all drawings are kept. 
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I. 



THE NEAR SURFACE PROBLEM 



A. MISSIONS OF NEAR SURFACE SUBMARINES 

Knowledge of the reason for a submarine to be near surface 
helps us to gain some insight into control requirements. Some 
noteworthy operations that a near surface submarine in straight 
motion is expected to be capable of are; 

1. Periscope depth for navigational aiding observations. 

2. Pretorpedo run optical sightings. 

3. Transmission. 

4. Missile launching. 

The demand for each one in a combat environment varies and 
satisfactory operation is directly connected with relative 
convenience in choosing speed, particular course with respect 
to the direction of prevailing seas and accepting limits on 
pitch and heave amplitudes. 

The near surface submarine cannot neglect the effects of 
surface waves. Furthermore, the problem is compounded in 
comparison with that encountered by the surface ship, since 

the operator cannot visually detect the approaching 

« 

disturbance . 

The outstanding characteristic of the open ocean is its 
irregularity. The irregular sea can be described by statis- 
tical mathematics on the basis of the assumption that a large 
number of regular waves having different lengths, directions 
and amplitudes are superimposed. Also the irregular motions 
of a ship in a seaway can be described as the superposition 
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of the responses of the ship to all the wave components of 
the seaway. 

B . WAVE THEORY 

In the theory of surface waves, simple two-dimensional 
waves are considered, which can be described completely by 
the motions, pressures, and so on, in any plane perpendicular 
to the crest lines. It is assxamed that (a) the crests are 
straight, infinitely long, parallel and equally spaced and 
that wave heights are constant, (b) water has zero viscosity 
and is incompressible. 

To facilitate the theoretical treatment of ship motions 
in regular waves and the use of the statistical approach to 
the study of behavior in actual sea states, the simple har- 
monic (sinusoidal) wave was introduced. 

The surface wave is the visible manifestation of pressure 
changes and water particle motions, affecting the entire 
body of fluid, theoretically to infinite depth. 

In Fig. I-l, the parameters and properties applicable for 
a sinusoidal wave on the free surface are displayed. 

9, the velocity potential is defined as a function whose nega- 
tive derivatives yield the velocity components of the fluid, 
and from this, all of the desired wave characteristics can be 
derived. 

For a two-dimensional wave in any depth of water 

9 = '^(x - x) 
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where 



X = horizontal coordinate, positive in the direction of wave 
propagation. 

z = vertical coordinate, positive downward. 

"^a = surface wave amplitude 

Lw = wave length 
h = depth of water 
k = wave number 2TT/Lw 
Vw = wave velocity 
t = time . 



For the case of deep water, the ratio 

cosK _ e + e 

V.H 






approaches 



- KX 



and then sm k(>^ 

Then = \A = horizontal component of water velocity 

=■ -xTaVcoe cosK(^X-^«"'^') 



“Ox 



= w = vertical component of water velocity 



= KlaVgj V. (x-V-ji) 



In deep water all particles describe circular paths having 

radii at the surface and decreasing with depth in propor- 
_ vz. 

tion to S. 

2- Q 

It can be shown that Vuj - which defines the 

the velocity of wave in any depth. In deep water 

” K 2, n 

For many problems, the most important aspect of waves is the 



distribution of pressure below the surface. 
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The elevation at lines of equal pressure in a wave 

relative to the still water pressure lines is shown to be 
CO;,K cos and for 

cosh v:h 

deep water COsY 

.xr 

Since e decreases as z increases, the contours of equal 

pressure are attenuated with depth. 

To obtain the surface wave profile, z is taken equal to 
zero. Then 

- Jq C05 v:(x-\/wi') 

-^a COs(kX-wi') where oJ ^ ^ 

Most of the characteristics of the simple harmonic wave-- 
except the precise profile — are the same as those of the more 
exact wave formulation. 

The contours at constant pressure that have been presented 
also indicate the increase or decrease in pressure relative 
to still water at any point in terms of depth or head. That 
is , p - p-c^ (^2.-'^') snd in deep water 'ip= cOS('<X'<^’t) • 

z should be measured to the center of the circular path de- 
scribed by the particle at the point of question. 

Application of the formula shows that under the crest the 
pressures are decreased and under the hollow are increased. 

The energy in a train of regular waves consists of kin- 
etic energy associated with the orbital motion of water 
particles and potential energy resulting from the change of 
water level in wave hollows and crests. 
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It can be shown that for a simple cosine wave the wave 
energy is half potential and half kinetic when averaged over 



a wavelength. The total energy is: 



E = Vsi.Pg^a'-^ or if the average energy per unit 

area of surface is considered. 




Stokes and others formulated the hydrodynamic theory of 
waves of finite amplitude. It corresponds with the observed 
fact that actual waves have sharper crests and flatter 
hollows from the simple cosine wave. The approximate result 
is that the surface profile simple cosine curve is modified 
by a harmonic which is half the length of the fundamental. 
This formulation leads to a limiting wave height from crest 
to trough of 0.14Lw or about l/7Lw. Real waves will break 
as this height is approached. 

C. FACTORS AFFECTING THE CONTROL OF THE NEAR SURFACE 
SUBMARINE 

The forces on a submerged body may be expressed as the 
integral of the pressure taken over the surface of the body. 
Under the assumptions already made (ideal fluid, etc.). 



The last term above is connected with the contours of constant 



it is seen for example, that for a submarine mean depth of 
50 ft. there can still be a significant wave amplitude, 
especially in moderate or high sea states. The problem 




pressure. From 
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examined as one of hydrostatics gives ^ - p c^-V -\^ 

where V is the volume of the body and K is the unit vector 
in the Z direction. Ideally, this passage of a pure 
pressure wave without the accompanying particles velocity 
would create essentially no force, since the integral of the 
pressure would not change but remain constant. 

Next there is a component of the pressure equation which 
results in the submarine being sucked towards the surface. 

Part of this force is usually referred to as a venturi effect, 
although the last is applicable for rigid boundaries and 
steady flow. 

Assuming infinite depth beneath the submarine, there will 
be no constriction of flow around the lower half of the boat. 
Assuming also that the submarine is in a uniform flow field 
with no disturbance, there will be a slight reduction in the 
cross sectional area of flow around the upper half, resulting 
in a higher velocity and thus a lower pressure. 

The reduction in cross sectional area is lessened due 
to the nonrigidity of the boundaries and the wave that the 
si±>marine will create on the surface. In reality, the water 
particles, due to existing surface irregularities, will have 
orbital velocities, another discrepancy of the venturi effect 
steady flow assumption. 

Consideration of a regular sinusoidal component of the 

irregular seaway shows that the velocities of the orbital 

2 

particles increase the magnitude of the V term in the pres- 
sure equation. So: 
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= VsVi'P + '-J-o COs(^'<.X-cO"^"^ 

= V5V\p t U.O QO-^OJ^+O ooKC^Q. d Vb cxx*bi.-\xa.vv^ 

O.’O^ ~ S — \J. 5 U'Z-{ • * S 

Squaring and averaging over time 

X N / ^ / 2. -2.KZ-0 

Vv-cd = Vshlp + V^, "Usiac^ mox .£ 

This effect coupled with the previously mentioned venturi 
effect can create as much as ten to twenty tons of lift on a 
submarine . 

Provided that the sea state remains unchanged, this is a 
time averaged force which remains relatively constant in time. 
Even if the essential irregularity of the sea has been recog- 
nized, it is a fact that over a wide area and often for a 
period of several hours, the sea may maintain a characteristic 
appearance which defies precise description in terms of aver- 
age or typical wavelength and height but which nevertheless 
is constant or steady. At other times or places, the sea will 
be quite different and yet there will again be a characteristic 
appearance. These observations suggested the possibility of 
statistical description of the sea and this was found to be 
quite feasible. The stochastic case can be modeled using 
energy density analysis. The spectrum can provide a mean 
value for the frequency and the wave height. The expected 
height and frequency can be combined to form an average or 
expected wave. This wave can be used to predict the average 
suck due to particle velocities. 

Upon encountering this force, the normally stable sub- 
marine becomes unstable and seeks a new equilibrium position 
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by ascending to the surface unless proper precautions are 
taken to prevent "broaching." 

There are basically two ways to offset suck. Adding 
ballast presents the difficulty in determining the proper 
amount for various sea- states without making the submarine 
dangerously heavy and sluggish. Using an angle on the sail 
planes to create a downward force presents the problem of 
reduction of control authority, saturation of planes and even 
losing the downward force if the orbital particle velocity 
is such that the angle of attack on the planes is reduced to 
zero. In addition, the particles also induce time dependent 
forces and moments, pitch and heave in particular. 

In Ref. 2, a qualitative example of the happenings as a 
submarine passes under a regular sinusoidal wave is given. 

The additional term in the pressure equation is 

usually neglected. 

An overall survey of the excitations experienced by the 
near surface submarine reveals that they may be broken into 
two basic categories; The first is a d.c. component in the 
form of a suck and the second is on a.c. component which is a 
combination of the various oscillatory forces. When the body 
is subjected to a stochastic seaway, accurate simulation and 
prediction becomes quite involved if not impossible. For a 
submarine under a seaway, the following must be noticed also: 

1. While at great depths, some of the partial derivatives 
are intuitively zero, for instance. Near the surface, the 
Z force varies with depth even in calm seas. The mere 



presence of Z^Z alone is sufficient to make the system unstable 
in the vertical plane. 

2. Under a seaway, the frequencies of heave and pitch 
will be identical to those at the encountered excitation. Then 

lexci-l - Tvnax ^\n(cO€>^ + \ ^ 

; ^ (oJe = oi encout^\et,j 

V\e%cd - n’fqQxCOsQuJet -r t J 

for a regular crested wave. When a stochastic excitation is 

introduced and assuming linear superposition, 

”Z SXCll - 2 I S’lO t + sO 

= 1 ^ 
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rilexcit = X, CosfcJe.tt t-C) 

Usb revjeoe consfavi-i 

ond ^ 

c — 

7eccU c 7. SinCooeii; t Z'O 
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These infinite sums are impossible to work, so they have 
to be replaced, if simplification is desired, by some known 
frequency and wave height statistical averages for the exam- 



ined sea state. 



II. STERN PLANES ONLY DEPTH AND PITCH CONTROL SYSTEM 



In order to keep the simulated stern planes only submarine 
on ordered depth and pitch in the presence of disturbances and 
to effect depth changes, a depth and pitch controller is 
designed- 

Figure II-l illustrates the assumed structure for the con- 
troller optimal design. The application of the optimal control 
is very well served by the linearization of the equations of 
motion. For motions in the vertical plane, roll, pitch, heave 
and surge are applicable. The method reduces the problem to 
four equations. 

For a completely submerged submarine, the effects of surge 
are small enough to be neglected and roll is not considered. 

For a submarine near surface, the problem of roll is worthy of 
an entire paper by itself. To simplify matters, it was neg- 
lected. For a submarine in head or following seas this is a 
very valid step. The state vector for the controller design 
was therefore ” i 
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0 

and the model for determining controller gains 
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Details on linearization of the equations of motion are 



given in Appendix A. 

XordL 

Using r = O for the command inputs, then for 

- ©oto. 

I J 

small perturbations x can be replaced by E, the state error 
vector and the problem solved as a linear regulator one, with 
states and control unbounded. 

The cost function to be minimized will be 

-C 5 9 E ^ u fi] dt 
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Forming the Hamiltonian 
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Necessary conditions for optimality are: 
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will finally give 
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It can be shown that the matrix K satisfies the matrix differ- 
ential equation 

k CO -- - K 5 - Q ft 5 " K 

with the boundary conditions: 

^C'tO^ Q Note that K is symmetric. 



For a plant completely controllable, with A, B, R and Q 
constant matrices, K > K (a constant matrix) as if -^00 . 

That is for an infinite duration process, the optimal control 
is stationary. From a practical point of view it is feasible 
in some cases, to use the fixed control law even for processes 
of finite duration. Then the K matrix can be determined also 
by solving the nonlinear algebraic equations (k = o') : 

O ^ - K ft - fT k -Q + \< ft ft" \< 

Using the method followed in Ref. 1 . 
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Finally 
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Integration backwards in time is executed since \^(ii)-^r 
until a steady state solution is obtained. 

Setting T! - t ^ d^I = - c±il . oj ^ '^ - O 
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dlt 



60 




I 



The signs of all equations are reversed and the initial 



conditions are those at ^ . Then 

tv= "Vr lb 1<5. = - Vc I f u fa 1 I. 



or 



D5RD^ 



Vu 



V2\/u 



V5\ 




E 



For the above, small changes are required in the original 
program (program 2 of Ref. 1) . This is also due to the fact 
that the approximate proportionalities observed in Ref. 1 

between the K’s and powers of U (speed) were again found valid. 

2 

^12' ^14' ^23 ^34 inversely proportional to U . 

3 

^22' ^24' ^44 inversely proportional to U . 

The controller gains depend on the selection of weighting 
factors. It became apparent that a method was required that 
would indicate at least the region from which the factors 
could be selected. Such a method and the reasons that led to 
it are described in Appendix C. Using this method, two sets 
of weighting factors were found, the one with C=10, B=800, 

E=l, and the other with C=0.001, B=164, E=0.001, which gave 
to the resulting controllers the desired characteristics. 

The weighting factors C=10, E=l, D=3000 (from Ref. 1) were 
also tested with the submarine at a speed of 15 knots and with 
no ballast added or removed. The controller failed to keep 
the submarine at the ordered depth. It was found that if the 
submarine were brought "in trim," then the depth keeping 
ability of the controller was considerably improved. 

In Appendix D, a simplified trim analysis is presented to 
help in the calculations for the required ballast for the "in 
trim" condition, at different speeds. 



For the stern planes only submarine, small changes were 
required in the rest of the simulation program. Db was set 
equal to zero and the bow plane angle generator was removed. 

A short description, connecting the methods of solution of 
the nonlinear equations of motion followed in Ref. 1 and 5, 
is given below. 

First step is the matrix formulation of the system equa- 
tions. In Ref. 5 
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above 


matrix. 


being a 


matrix 


of constants, is inversed 



and multiplied b'^ the vector in the right hand side of the 
matrix equations, and the resulting six differential, 
linear equations are solved simultaneously. 



non- 



III. BOTH PLANES DEPTH AND PITCH CONTROL SYSTEM 



A. OPTIMAL SOLUTION, UNBOUNDED CONTROLS 

The approach selected for the design of the controller, to 
maintain automatically depth and pitch, was to let the planes 
individually achieve whatever angle was required to meet an 
optimal control based on the minimization of the cost function 



where E = state error vector 
a = control vector 
Q,R = weighting matrices 

That is, it was desired to maintain the state vector close to 
the ordered trajectory, without excessive expenditure of con- 
trol effort. It was assumed that admissible controls and 
states are not constrained by any boundaries. 

This approach was used also in Ref. 1 and for this reason 
only the required changes after the selected linear model for 
stern planes only submarine will be presented. 

The changes originated in the inclusion of term -UL’Q" 
in the state equations and resulted in the following altered 
form of the gain derivative equations: 
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^l'^2 included for the case that different weighting on 

the planes effort would be desired. As used in the simulation 
runs of this paper, C was set equal to C^tC^- 

B. CONSTRAINTS ON ADMISSIBLE CONTROLS 

The same controller in an alternate form can result if 
bounds are defined on the controls and Pontryagin's minimum 
principle applied. 



Under the assumption of independent controls, the quanti- 
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The two controllers are equivalent. The one will be 
referred to as unbounded and the other as bounded. 

C. LINEAR TRACKING CONTROLLER, EXACT DESIGN 

Until now the procedure of Ref. 1 was followed and the 
initial tracking problem was converted by use of simplifying 
approximations to a linear regulator one, where it was desired 
to maintain the state error vector close to the origin. 

From an attempted more exact examination of the linear 
tracking problem, a modified controller resulted whose capa- 
bilities are considered greater. This is because approxima- 
tions such as use of instead of the correct 

are avoided. 

Forming the Hamiltonian [Ref. 5] 

K ^ 11 + '/2 11 U Hn + ft X ^ ft u 

u = - li" £ 

— ^ X “ 

and for p = ^ ^ 

0^--^ 5.— Qo K^te 

i = ^ ^ '5"] 1 ftQ r 

with = H ^ sOftl - - VI 

Integration from t^ to t^ gives K and ft With H = 0 K results 
as before. 

Using the symbols and notations of Ref. 1 then, 

^ =(pT- 1 - Q :C ' again 
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In the expression for U. the part -R 5. ^ X is the same 

as in the regulator case with X the state vector instead of 
E the error vector. So, finally 
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This is 


easily shown from 
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The controller now becomes : 

-pS. c \4.? . OEPTH + • 70D01 + Vc25 • PlTCl-l + -PlOOT +S2 

- VHM 'DEPTW 4-k.2^ 'ZODOT -PlTCU +- •P^OO'T 4- SA 

• E> M V p t, . fea 0 /c J 

aeoB = 4 PA -&95VC2 

etc. 

D. TESTING THE CONTROLLERS 

All tests v;ere run at the forward speed of 15 knots and 
with the values of the weighting factors C=10, D=3000, E=l, 
to provide results for easy comparison with that of Ref. 1. 

To check the depth keeping capabilities of the bounded 
form of controller under the presence of external disturbances, 
a ramp input AU = -2.E-05 -TINE , was applied at the auxiliary 
trim tank (Fig.iV-A ). The results are shown in Fig. III-l 
for the otherwise unballasted submarine. 

The controller kept the submarine at near zero depth until 
the fairwater planes became saturated. Then the "light" 
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submarine started losing depth quickly. The following must 
also be pointed out at this point: In the unbounded controller 

of Ref. 1 and in order to limit the effect of initial depth 
error (corresponding to a large step order) resulting in large 
plane angles, unachievable on a real ship, the depth and pitch 
errors were limited as shown in Fig. III-2. The limiting 
values tested were i 2 feet for depth and i 10° for pitch 
errors. Also to reduce the noise level in the plane positions, 
the plane angle ordered was filtered with a simple low pass 
filter. 

In the test of Fig. III-l, the limiters were preserved in 
the bounded controller but the filters were removed. Thus as 
seen in Fig. III-lc., although the submarine reaches a depth 
of -80 feet, the reaction of the unsaturated stern planes is 
very slow. 

From Fig. III-3 and III-4, the effect of the low pass 
filters can be seen. Fig. Ill- 3 shows the response of the 

submarine with bounded controller and without filters to a 

- 4 

sinusoidal input, AU=2.0xl0 x sin 0.8*t. 

Figure III-4 presents the response to the same input when 
the unbounded controller of Ref. 1 as modified in this section 
is used. Except for noise components, the results are other- 
wise similar. 

To compare the effectiveness of the controllers to depth 
changes, an ordered depth of 10 feet was 'used. Figure III- 5 
shows the response of the submarine with the unbounded con- 
trollers, with limiters and filters. 
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Figure III-6 presents the response of the submarine with 
the bounded controller, with limiters but no filters. Again, 
no appreciable difference is detected, as was expected since 
the controllers are in fact the same in that region of opera- 
tion. Both responses were acceptable. 

In Fig. III-7, the submarine was ballasted with AU= 

8.8X10 FT=6.4X10 AT=-6.4X10 The bounded controller 

was tested for its ability to keep depth. To check this, the 
ordered depth was put to zero. The controller maintained a 
steady state depth error of 0.05 feet and a steady state pitch 
error of 0.0002vod.A steady state fairwater plane angle of 
about -0.7° was required to maintain depth. This small value 
implies that the simulated submarine was almost "in trim." 

The next test was again with the submarine unballasted. 

A depth of 10 feet and a pitch of -0.150 rad (-8.6°) were 
ordered with the tracking controller in use. The error lim- 
iters and noise filters were removed. Responses are shown in 
Fig. III-8. Although the resultant steady state values of 
10.4 feet and -9.8° could be considered acceptable, the sub- 
marine overshoots to a depth of 24 feet before achieving its 
final depth. The negative pitch was ordered in order to 
speed up the depth change but on the other hand this created 
difficulties for the planes in their action to stop the 
pitching submarine when it reached the ordered depth for the 
first time. 

In Fig. III-9, the same maneuver was repeated while a force 

AU=2X10”'^ sin 0 . 8t was also applied on the boat. The results 
were similar. The depth overshoot was not acceptable. 
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since it was desired to be able to use suitable pitch of 
the submarine when ordering depth changes, the method of pre- 
senting the orders to the tracking controller was modified 
and the following statements were added to the program: 

PI I = ■ -0.055 
PI5 = 0.052 

hOCG- P\5S CoCHivo') 

= 5.0. 4- Pi OCG/R P\T £ 

VI ^ a prop ( ao-") 

- p>prn P C ov") 

200P = p • OCH N&/p( 005 

POpD = t^<O^T Cpii ,Pl2.voi 
The inputs were now the desired depth change and the rate of 
depth change (DCHNG and RATE) . Suitable pitch between PI^ and 
is automatically ordered during the depth transition. The 
values of RATE, Pl^, PI^ and the form of XD can easily be modi- 
fied to achieve desired depth change characteristics. 

In Fig. III-9, results are presented for the ballasted 

-4 

submarine, the sinusoidal input AU=2X10 sin 0.8t super- 
imposed and an ordered depth change of 50 ft. The final depth, 
achieved without significant overshooting, was 51.7 feet and 
the pitch 0.019 rad ( ~ 1°) . The values of the parameters 
RATE, Pi^, PI^ , as appear above, were selected to give a response 
similar to that of the submarine in Ref. 1 with the unbounded 
controller, the error limiters and noise filters. This was 
verified by the results. 
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Although setting bounds on the planes may lead the boat 
to unadmissible limit cycles, the nature of the near surface 
problem prohibits the use of the unbounded controller with 
error limiters when improved performance is desired. 

The tracking controller, due to its more exact design, 
is expected to have better performance. With the introduced 
method of ordering depth changes, the error limiters are no 
longer necessary. Further investigation is required on the 
operational characteristics of the tracking controller. 

In the rest of this work, the bounded controller is mainly 
used due to the fact that it was developed earlier. Lack of 
time prevent duplication of the tests with the tracking con- 
troller in use. 

In Fig. III-IO, the block diagrams of plant and controller 
in the linear regulator problem and in the linear tracking 
problem are shown together to depict differences. 

In Appendix C, the method developed to connect linear feed- 
back with optimal control theory is examined for the possibil- 
ity of its application in the case of the both planes submarine 
It may be repeated here that the results are rather prohibitive 
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IV. CONTROLLER FOR THE NEAR SURFACE SUBMARINE 



A. FORMULATION OF THE CONTROLLER. COMBINED MODE 

CONTROLLER (CMC) 

The need for an improved depth keeping performance for the 
near surface submarine under a seaway led to an investigation 
for an "optimal" near surface controller. 

The inadequacy of the already designed optimal unbounded 
and bounded ones is easily and simply demonstrated. A 70-ton 
pulse input is applied at the auxiliary trim tank and the 
unacceptable results of submarine response shown in Fig. IV-5 
and 6 are discussed in the simulation section of this chapter. 

The phenomena governing submarines response to submerged 
turbulence are not completely understood and much of the 
existing experimental data have been classified. Thus many 
simplifications were accepted in the formulation of the prob- 
lem and a great effort was initially devoted in the direction 
of finding a mathematical proof that would suggest an optimal 
solution to the simplified problem. Even this task was proven 
very difficult to the author and a mathematical solution was 
unattainable in the time allocated for this work. 

The next step was toward the examination of a proposed . 
controller scheme, for which existing test results were known 
to be encouraging. In that scheme, two separate control sys- 
tems are designed, one based on having both planes available 
and the other based on stern planes only control. The two 
stern planes commands are then computed and combined as a func- 
tion of the forward plane command. 



Defining = fl (Depth Error, Depth Rate, Pitch Error, 

Pitch Rate) 

= stern plane command in both planes mode 

= f2 (Depth Error, Depth Rate, Pitch Error, 
Pitch Rate) 

= stern plane command in stern only mode 

1^-5 CO'T^ = = f3 (Depth Error, Depth Rate, Pitch Error, 

Pitch Rate) 

= bow plane command in both planes mode. 

ihen the actual stern plane command formed o-^ a linear combi- 
nation of Ss\ is 

(bsco'f^ = VC -6 SI ( 1 - (bs'l coKpt^ 

i \f I <C6-^-nnax 

(i.e. bow planes not saturated) 

V,- CO 'f’l, 6^1^00 Otherwise 

The function K = was used with ^ ^vriCix = 

\ dico^h/ 

and X a variable parameter. 

The problem now was to try to find a "best" value for X 
for most acceptable controller performance, test the general- 
ity of this value in various kinds of controller designs and , 
if possible, also find a mathematical proof of its uniqueness 
or at least a proof of the existence of a restricted region 
of permissible values. For convenience, the following abbre- 
viations will be used in the next. 

S 0 P 0 C. = stern only planes optimal controller 

S 0 P C = stern only planes controller (not optimal) 

B P 0 C = both planes optimal controller 
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CMC 



= combined mode controller (meaning the 
proposed scheme) 

It is again noted here, that the S 0 P-0 C. and B P 0 C 
designed in Chapters II and III are in fact optimal only under 
the assumptions of unbounded states and controls. 

The synthesis of the CMC, was realized at first by the 
use of the already designed S 0 P 0 C. and B P 0 C. Then the 
S 0 P 0 C was replaced in the CMC byaSOPC. corres- 
ponding to a characteristic equation, same with that result- 
ing from the B P 0 C. The formulation of the CMC. from the 
independently designed SOPOC orSOPC and B P 0 C ones 
is indicated in Fig. IV- 1 to 3. 

In the next section of this chapter simulation results, 
mainly of the use of C M C on a submarine excited by external 
forces, are presented. Parallel analysis of the data on indi- 
vidual and comparative bases is attempted. 

Some conclusions appear proper at the end of the chapter. 
The part concerning the question of the generality of a 
unique "best" value or region of values for the CMC. param- 
eter is left for the next chapter. 

B. SIMULATION PESULTS 

The computer used for the simulation was the I.B.M, 360/67 
located at the W. R. Church Computer Facility of the Naval 
Postgraduate School. 

Program 2, developed in Ref. 1 as a modification in D.S.L. 
of the original N.S.R.D.C. digital program of Ref. 3, is the 
main program used. The hydrodynamic coefficient values were 
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kept the same and all runs were made in forward motion, at a 
constant speed of 15 knots (25.33 ft/s) and with the submarine 
approximately "in trim." Appropriate modifications in the 
equations that lead to the control orders to the planes were 
made occasionally, as for example in the case of the S 0 P 0 C. 

In all cases the linear model, derived in Appendix A, was 
used. The trim-controller part of the program was not 
operated. 

To simulate the effect of an external forcing term, water 
was removed or added to the trim tanks whose position as de- 
fined in Ref. 1 is shown in Fig. IV-4. The depth and pitch 
error limiters used in the controller of Ref. 1 had to be 
removed when operation of CMC. was considered. They were 
used initially in the unbounded controller in order to avoid 
excessive values of the planes in response to a large step 
order. Continuation of their use will undesirably alter the 



affect the operation of the CMC. 

At first the narrow pulse force shown in Fig. IV-5a was 
applied on a submarine (forward trim tank) using the B P 0 C. 
with bounded controls and with the limiters removed. The 
controller had weighting factors B=800, C=10, E=l. 

From Fig. IV-5b to IV-5e, it is seen that the controller, 
after the saturation of the bow planes, acted slowly to sta- 
bilize the submarine with a steady state depth error of about 
-65 ft and a steady state pitch error of -0 . 0 86 Y’ad,.The bow 
planes never recovered in the time interval of 200 sec. the 




and consequently will 
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simulation was held. In the case of the bounded controller 
with B=164, C=0.001, E=0.001, application of the same input 
was destructive as shown in Fig. IV-5f to IV-5i. 

The controller with B=800, C-10, E=1 was also tested with 
the limiters replaced. The controller succeeded now in re- 
turning the boat to the ordered zero depth and pitch but with 
an overshoot of about -30 ft in depth (Fig. IV-6a to 6d) . 

When the weighting factors were changed to C=10, E=l, D=3000, 
the controller failed completely again (Fig. IV-6e to 6h) . In 
all cases the results were unacceptable and simply show the 
inadequacy of the designed B P 0 C with bounds on the con- 
trols and with or without error limiters, to handle situations 
approximated by the application of the pulse force. 

All of the subsequent test runs in this chapter simulate 
motion of the submarine using the CMC Only three kinds 
of external excitations were considered, the pulse force at 
the forward trim tank, shown in Fig. IV-5a, a step of the same 
magnitude at the auxiliary trim tank, and the same step at the 
forward tank (Fig. IV-7a) . 

1. Application of Step Force at Auxiliary Trim Tank 
Since it was desired to check the operation of" the 
CMC after the saturation of the forward planes, the magni- 
tude of the step was selected big enough to cause saturation. 

The CMC was formed by the combination of a S 0 P O C 
and a B P 0 C (bounded) having the same weighting factors. 

Under the above conditions, two controllers were 
tested. The first one with weighting parameters B=800, C=10, 
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E=l; the second with B=164, C=0.001, E=0.001. There is a 
distinct difference in the resulting controller orders, which 
was the reason for the selection of the above two sets of 
values. The optimal controllers, when the approximation of 
infinite duration control interval is used, result in time 
invariant state feedback given in the program in the form (for 
the order to stern planes) : 

■'ZOfd vN^i/uc •'ZoOOT +V5I-PERR j/u c • P ' OoT 

where Zo€P = Dep-t-K evvov 
20 00T - DeptU valf 

PepR = P'\VcK -evYor 
= P\Vc\\ raVe 

and '^ 21 ' ^31' ^41 linear combinations of the gains 

in the K matrix. 

The values of etc. in the tested controllers are 

given in Table IV- 1. 

Using the closed loop C.E. derived in Appendix C, it 
is seen that of the first controller, if used as the order 

to the planes in a stern planes only case, is expected to 
result in a stable submarine (roots with negative real parts) . 

of the second controller in a similar use is expected to 
result in an unstable submarine (roots with positive real 
parts) . 

For each of the combined mode controllers, different 
values of X were tried. The first controller showed improved 
response in terms of steady state depth error, as X was in- 
creased from 0.005 to 500. Results are plotted in Fig. IV-7 
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to IV- 15 and arithmetic values of important response charac- 
teristics are collectively presented in Table IV-2. From this 
table it is seen that for X=0.3 the steady state depth error 
was approximately -6.4 ft and the steady state pitch error 
-2.9 X 10 ^ rad (1.6°). For X=500, the s.s. depth error was 

approximately -4.52 ft and the s.s. pitch error again -2.9 x 
-2 

10 rad. Both responses are considered acceptable and so, 
for this controller the steady state error criterion does not 
impose any strict limits on X. Control effort differences 
are not easily distinguishable and simple inspection of the 
stern plane angle plots is not enough for conclusions. 

When the value X=0.5 was used in the CMC, the stern 
plane angle attained its minimum s.s. value. 

The second controller (B=800, C=0.001, E=0.001) was 
tested for the values of X between 0.5 and 6.0. Numerical 
results are given in Table IV-3. For values of X greater than 
3.75, the system was unstable. 

In Fig. IV-16, the response characteristics of the 
submarine with the CMC controller and X=6 are plotted. The 
plots are representative of responses resulting from the CMC 
with values of X greater than 3.75. Responses for X less than 
3.75 are plotted in Fig. IV-17, 18. For X=0.5, the stern 
plane s.s. angle reaches a minimum value. On the other hand, 
the s.s depth error is minimized for X values near 3.75. 

The region of acceptable values of X is considerably 
reduced in this controller. Again, steady state error and 
control effort criteria are not enough to distinguish any 
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"best" value at the exponent. Responses to positive or nega- 
tive step forces were similar (opposite sign) . 

2 . Application of Step Force at Foirward Trim Tank 

The step force was applied to the forward trim tank 
and each CMC tested again. Application of the force to the 
aft tank was not considered since it will saturate the stern 
planes first. 

Use of the first controller (B=800, C=10, E=l) with 
values of X greater than 0.9 resulted in unacceptable responses. 
The submarine followed a sinusoidal path between a depth of 
-30 to -45 ft. 

In the region 0.005 0.9 and among the tested 

values, X=0.7 gave minimum s.s. depth error of -25.15 ft, 
steady state pitch error of -0.14 rad (8°) and minimum stern 
plane s.s. angle (11.16°). 

Table IV-4 presents numerical results of the tests, 
while the corresponding responses are plotted on Fig. IV- 19 
to IV-23. 

In Fig. IV-19, the value of X=1.0 has been used but 
the curves are characteristic of submarine responses for values 
of X larger than 0.9. 

Use in the second controller (B=164, E=0.001, C=0.001) 
of X=0.5, gave the first indications of its stabilizing capa- 
bilities. After a transient oscillating period the sxibmarine 
reached a s.s. value of -28 feet in depth and -8° in pitch. 

Values of X greater than 0.5 resulted in a sinusoidal 
response of the boat around an average depth of -40 ft. 
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The response characteristics for X=0.9 are shown in 
Fig.- IV-24 and are representative of responses when values of 
X greater than 0.5 are used. 

From Table IV-5 and Fig. IV-25 to IV-28, it is seen 
that values of X around 0.4 result in relative better response 
characteristics. The value of X=0.005 is destructive for the 
boat, sending it quickly to the surface at an increasing posi- 
tive angle. It is noted that when the step was applied at the 
forward trim tank, permissible values of X became smaller and 
were restricted also in smaller interval. 

Application of the step on the forward trim tank 
results in higher values for depth and pitch rates than in the 
case where the step was applied on the auxiliary tank. This 
can be used as a crude indication that higher rates require 
smaller values of X, but in any case no indication was found 
that there is a unique "best" value of X. 

Since application of step input forces and moments 
and use of minimum s.s. errors were not considered satisfactory 
for the representation of near surface effects and determina- 
tion of the "best" value of the exponent X, the next step 
tested application of pulse forces (approximating impulse 
ones) at the forward trim tank. 

3 . Narrow Pulse Force at Forward Trim Tank 

The magnitude of the pulse force used was the same as 
for the step, i.e., about 70 tons. Its time duration was 
10 s. (Fig. IV-5a) . Both controllers with B=800, C=10, E=1 
and B=164, C=0.001, E=0.001 were tested. Table IV-6 presents 



some of the important characteristics of submarine response 
when the first controller is in use. Corresponding curves 
are shown in Fig. IV-29 to IV-36. For values of X larger than 
0.6, the submarine is unstable. Responses are similar to those 
indicated in Fig. IV-29, which correspond to X=0.7. Values 
of X from 0.005 to 0.6 gave acceptable results. There are 
small differences in depth and pitch overshoots but these 
differences are not sufficient to justify conclusions about 
optimality of a specific value of X. Values of X between 0.45 
and 0.6 result in slightly improved response characteristics 
in comparison with the rest of the permissible values, when 
depth overshoot is considered. 

Tests with the second controller gave, in general, 
unacceptable results. For values of X larger than 0.2 the 
boat was unstable. At the value of X=0.005, once again the 
controller fails and the submarine moves unhelped to surface. 
The results are shown in Fig. IV- 'bl to IV--H2, and in Table 
IV- 7. 

Judging from the tested values of X, it appears that 
the best to be expected after the pulse is removed, is a 
slightly damped sinusoidal response or even a small amplitude 
undamped one. For this to happen, the value of X is expected 
to be between 0.1 and 0.2. 

4 . Simulation with an Alternate Formulation of the CMC 

Existing information from previous experimental work 
by other sources revealed that the value X=0.5 could be "best" 
for a specific design of combined mode controllers where 
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(a) the both planes component-controller was designed 

optimal . 

(b) the stern planes component-controller was designed 
to give the same closed loop C.E. as in (a) . 

(c) (^gi. of the both planes-component-controller if used 
in a stern planes only submarine, would slightly stabilize it. 
The optimal both planes controller with B=800, C=10, E=1 ful.- 
fills condition (c) and the corresponding C.E. found from pro- 
gram 6. was used in program 4. to give feedback gains in SsZ 

of the CMC. The controller used in Ref. 1 with D=3000, 

C=10, E=1 also satisfies the condition (c) and was tested 
in C M C too. 

To complete the picture, the controller with B=164, 
C=0.001, E=0.001 was used in similar tests, although it was 
known that its order SsL would result in an unstable stern 
only mode submarine. 

In Fig. IV-43 to IV-47, the response characteristics 
of the submarine to a pulse force at FT are presented. The 
CMC uses a B P 0 C with B=800, C=10, E=1 and a S 0 P C 
resulting from the same closed loop C.E. 

Numerical results are collectively presented- in Table 
IV-8. Compared with the results in Table IV-6, they are 
inferior. The value of X=0.5 doesn't appear to be a "best" 
one. Rather a value of X=0.1 would be preferable for use in 
a controller. 

The submarine was also simulated with D=3000, C=10, 

E=1 in B P 0 C and S 0 P C with the same closed loop C.E. 

All the tested values of X (i.e. 1.0, 0.8, 0.6, 0.5. 0.4, 0.15) 

resulted in unstable operation. 
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Figures IV-48 to IV-50 and Table IV-9 present the results 
when in CMC , BOPC and S 0 P 0 C wilK the same weight- 
ing C=10, D=3000, E=1 are used. 

The last tests are with CMC that makes use of 
B P 0 C with B=164, C=0.001, E=0.001 and of S 0 P C with 
the same closed loop C.E. For values of X, 0.1, 0.3, 0.5, 
the submarine is quite unstable. 

In Table IV- 10, the permissible values of X for use in 
CMC are collected and presented in a condensed form. It is 
seen that only in the case of C M C using B P 0 C and 
S 0 P C are there enough cases with "best" value of.X around 
0.5. Even in these cases, the dependence of the parameter on 
the input force and on the design of the particular controller 
(i.e., weighting factors) is apparent. 

The strong indication that the "best" value of X lacks 
generality and that only strict specifications would result in 
the acceptance of a certain value (for example X=0.5) as opti- 
mum in some sense, led in Chapter 5 to a further investigation 
of the proposed CMC. 
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TABI,{£ IV ‘1a, Gains of designed controllers 
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TAIU 1. JV“1 C Goinf, uf (5e .iqn, :j centre 



X 


S . S .Depth 
error 
(ft) 


S.S. Pitch 
error 
( rad ) 


Time to ± 5 ?^ 
of the S.S. 
depth 
( s) 


Stem plane 
angle at S.S 
(deg) 


1000 


- 4.26 


-0.027212 


105 


-7. 


500 


- 4.522 


-0.02875 


170 


-6 . 564 


■16 


- 4.587 


-0.02876 


165 


-6.563 


8 


- 4.682 


-0.02876 


162 


-6 . 564 


6 


-^.737 


-0.02876 


165 


-6.563 


4 


- 4.835 


- 0.02877 


160 


-6.554 


3.5 


- 4.688 


-0.02751 


■ 175 


-6.55^ 


2.5 


- 4.981 


-0.02878 


160 


-6.553 


1.5 


- 4.993 


-0.02751 


175 


-6.555 


1 . 


- 5.2 


-0.02751 


175 


-6.563 


0.9 


-5.263 


-0.02751 


180 


-6.553 


0.7 


-5.649 


-0.02877 


167 


-6.553 


0.5 


-5.719 


-0.0275 


177 


-6 .125 


0.3 


-6.393 


-0.02875 


160 


-6.988 


0.005 


- 14.126 


-0.029 


240 


-6.988 



TABLE IV -2 



Simulation results of the application of a ? 0 -ton 
step force at the auxiliary trim tank. The submarine was 
"in trim" and the CMC used optimal des'-^n in the combined 
controllers .with B= 800 . ,C= 10 . ,E =1 . 
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X 


S . S .Depth 
error 
(ft) 


S . S . Pitch 
error 
(rad) 


Time to ±5^ 
of the S.S. 
depth 
(s) 


stern plane 
angle at S.S 
(deg) 


6 


Unstable su 


Dma.rine 


- 


- 




Unstable su 


bmarine 


- 


- 


3.75 


Unstable su 


bmarine 


_ 




3.5 


-5.462 


-0.02886 


170 


- 6.825 


2.0 


-5.48 


-0.02886 


170 


- 6.6 to - 7.2 


0.5 


-6.059 


-0.02886 


170 


- 6.608 



TABLE IV-3 

Simulation results of the application of a ? 0 -ton 
step force at the auxiliary trim tank. The submarine was 
"in trim" and the CMC used optimal design in the combined 
controllers with B=l64 . , C=0 . 001 ,E=0 . 001 . 
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X 


S . S .Depth 
error 
(ft) 


S . S .Pitch 
error 
( rad ) 


Time to + 5 ?^ 
of the S.S. 
depth 
( s) 


Stern plane 
angle at S.S. 
(deg) 


500 


Sinusoidal 


variation 


- 


- 


50 


Sinusoidal 


variation 


- 


- 


10 


Sinusoidal 


variation 


- 


- 


2.5 


Sinusoidal 


variation 


- 


- 


1.5 


Sinusoidal 


variation 


- 


- 


1.0 


Sinusoidal 


variation 


- 


- 


0.9 


-25.5 


- 0 .14034 


130 


11 .16 


0.8 


-25.79 


- 0.14036 


170 


11 .598 


0.7 


-25.15 


- 0.14033 


110 


11 .16 


0.5 


-27.25 


-0.14036 


140 


11 .598 


' 0.05 


- 48.78 


-0 .i 4 o 47 


150 


11 .16 


0.005 


-92.36 


- 0 .14032 


210 


11.377 



TABLE IV -4 

Simulation results of the application of a ?0 ton 
step force at the forward trim ta.nk.The submarine was 
"in trim" and the CMC used optimal design in the combi- 
ned controllers with B= 800 . ,C= 10 . ,E =1 . 
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X 


S . S .Depth 
error 
(ft) 


S . S .Pitch 
error 
( rad ) 


Time to 
of the S . S . 
denth 
( s) 


Stern plane 
ang le at S . S . 
(deg) 


3.5 


Sinusoidal 


variation 


- 


- 


0.9 


Sinusoidal 


variation 


- 


- 


0.5 


-28. 


-0 .142 


475 


12.31 


0 


- 29.362 


-0.l4l 


145 


11.86 


0.3 


- 30.694 


-0 .l4l 


145 


12.296 


0.05 


-113.04 


-0.l4l 


160 


11 .419 


0.005 


Unstable su 

. ^ 


bmarine 
i 


- 


- 



TABLE IV-5 

Simulation results of the application of a. ?0-ton 
step force at the forward trim tank. The submarine was 
"in trim" and the CMC used optimal design in the combi- 
ned controllers with B=l64. ,C=0 .001 ,E=0 ,001 . 
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X 


Maximum Depth 

over shoot( f t) 


Maximum Pitch 

overshoot( rad ) 


Time to positive 
depth pea,k( s) 


- 


+ 


- 


4 - 


10 


Unstable 




submarine 




- 


- 


0.9 


Unstable 
1 


submarine 


_ 


- 


- 


0.7 


Unstable 


submarine 




- 


- 


0.6 


7.56 


1 .26 


0.053 


0.012 


32 


0.55 


7.56 


0.826 


0.0523 


0.0094 


32 


0.5 


7.52 


0.826 


0.0496 


0.0096 


32 


0.45 


7.6 


1 .l 4 l 


0.0525 


0,012 


33 


0.3 


7.75 


1.377 


0.0531 


0.0137 


33 


0.05 


8.26 


1 .181 


0.0367 


0.0132 


33 


0.OO5 


8.42 


1.377 


0.0291 


0.0154 


33 



TABLE IV -6 

Simulation results of the application of a 70-ton x 
10 s pulse force at the forward trim tank. The submarine was 
"in trim" and the CMC used optimal design in the combined 
controllers with B= 800 . , C= 10 . ,E =1 . 



^3 



X 


Maximum Depth 

overshoot( f t 


Ma-Ximum Pitch 

overshoo t( rad ) 


Time to positive 
depth p6aJc ( s ) 






0.5 


Unstable submarine 


- 


0.3 


Unstable submarine 


- 


0.2 


Residual oscillatic 

8.26 2.913 1 


pns 

|o .0472 0.00425 


27.2 


0.13 


Residual oscillations 

8.42 3.94 |0.039 0.0164 


29.6 


0.1 


Residual oscillations 

8.89 3.86 [0.0269 0.0055 


32.8 


0.05 


Residual oscillations 

9.37 4.33 |o. 00157 0.0192 


39. 



TABIxE IV-7 

Simulation results of the application of a 70 ton x 10 
pulse force at the forward trim tank. The submarine was 
" in trim ” and the CMC used optimal design in the combined 
controllers with B=l64. ,C=0 .001 ,E=0 .001 



X 


Maximum Depth 

overshoot( ft' 


Maximum Pitch 

overshoot ( rad ) 


Time to positive 
depth peak( s) 


- 




- 






0.9 


Unstable submarine 


- 


0.7 


Unstable submarine 


- 


0.5 


11.14 


1.023 


0.075 


0.0125 


27 


0.4 


10.74 


0.157 


0.061 


0.0094 


28 


0.3 


10.31 


2.36 


0.049 


0.0244 


37 


0.1 


9.13 


1 .42 


0.0031 


0.0157 


33 



TABLE IV -8 

Simulation results of the application of a 70 ton x 10 s 
pulse force at the forward trim tank. The submarine was 
" in trim ” and the CMC used combined controllers 

with the same CE.(BP 0 C with B= 800 . , C= 10 . ,E= 1 . ) 



Y . 


Maximum Depth 

overshoot( ft) 


Maximum Pitch 

overshoot(rad) 


Time to positive 
<;iepth peak(s) 


A 


- 


-r 


- 


H- 




0.4 


7.637 


0.945 


0.0456 


0 .14 


30 


0.5 


7.637 


0.771 


0.0456 


0.018 


30 


0.6 


7.48 


1 .18 


0.047 


0.0145 


30 



TABLE IV -9 

Simulation results of the application of a ?0 ton x 10 s 
pulse force at the forward trim tank. The submarine was 
'' in trim ” and the CMC used optimal design in the combined 
controllers with D= 3000 . ,C= 10 . ,E= 1 . 



C M G 1 CMC 




o 



CO 

Cj 

o 

a 

> 
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table IY-10 



V. FURTHER INVESTIGATION OF THE 



COMBINED MODE CONTROLLER (CMC) 



A. ANALYSIS OF CONTROLLER OPERATION 

In Chapter IV a controller was devised for the near sur- 
face submarine. Its stern plane order was formulated by the 
combination of the stern plane orders of two separately 
designed controllers, one for a stern planes only submarine 
and one for a submarine having both planes available. 

The order took the form 



In one case, the two separate controllers were the result 
of optimal design with the same weighting factors for both. 

In the second case, while the both planes controller was kept, 
the stern planes one had its gains adjusted to give the same 
closed loop C.E. The controllers were used on a number of 
simulation runs and in order to compare their effectiveness 
under a sea way, step and pulse forces were applied on the 
submarine . 

It was desired also to determine a value or region of 
values of the parameter X that would correspond to improved 
performance of the system. If existence of such a value or 
region of values was apparent, then a theoretical justifica- 
tion would probably follow. 

Examination of the results shows that such a value or 
region of values is not likely to exist in general. On the 
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other hand, for a certain controller design and under certain 
inputs, values of X around 0.5 seem to offer better results. 

More specifically, as seen from Table IV-10, this was the 
case for the CMC using optimal separate controller design 
with weighting factors B=800 , C=10 , E=1 or D=3000, C=10 , E=1 
and when a narrow pulse force was applied at the forward trim 
tank . 

When the controller with B=164, C=0.001, E=0.001 was used, 
the results, not acceptable in general, showed improvement 
for values of X around 0.15. It appears that "optimality" of 
X is directly connected with the design of the separate con- 
trollers and the type of applied force. This then implies that 
generalization of the results is not possible and a specific 
value of X only accidentally can result as the "best" in a 
CMC with arbitrary separate controller designs. 

To further investigate that view, the following procedure 
was adopted: A linear relation is assumed between 

for example where f\^ is a constant. Then 



cation of the external forces, a stern planes only submarine 
configuration has resulted. To each value of M, a stern planes 
only controller corresponds. 




So under the assump- 



tion of a linear relation between the orders > the 

result is that i^scoyt) can be expressed as a linear combination 
of 6$2. only. Sscott] is used in C M C when due to the appli- 
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From this family of controllers the one or the ones with 
the most desirable response characteristics to certain inputs 
can be found by tests on the stern planes only submarine. 

Then the "best" value of X in the original expression of 
will be the one capable of keeping the values of M near the 
selected one, for all the expected values of S-^CO^i . In this 
way, the dependence of X on the controller designs and the 
applied inputs will have been shown. 

The above can be considered as a short description of what 
now is presented in detail. 

From optimal design with the assumption of infinite dura- 
tion process and in general from the case of state feedback 
design, Ss2. results in the form 

^si-SoPi-Xi 4-So%-K2. V SOC-X3 +‘boO'Xii where 

X^, X 2 / X^, X^ are the states and SOA, SOB, SOC, SOD are the 
gains. The closed loop C.E. as found in Appendix C is given 
by: 



s' + 






An 



1^12) ' -5015- 1!)2-Soo 







(au- A 9 ?- R)2- -&t - Pi2i fe-a • soc 

4 - (Pl H -(52- R 1 2 . So D 



bi ■ fAa I + • So P\ SOC 



(5fUc-P\.2 - tSs -Uc* SOPi 



9S 



Then the order <^sco’nn= M • (^ sQ. / if substituted in, will result 
in variation of the position of the eigen values of the above 
equation, since the gains are multiplied by the parameter M. 
For each set of values of feedback gains corresponding to a 
certain design of stern only mode controller, the characteris- 
tic equation takes the form: 

+C2.')- s' t ('c-nr) -rC-^-S f 

For the given submarine at the velocity of 25.33 ft/sec: 

= 0. 1982 ■ 

C 2 = 9.7316 X ^0■^ 




A, B, C, D calculated in program 8 are functions of the 
feedback gains. This form is suitable for root loci methods. 
The loci were plotted for different controller designs. It 
must be noted here that root loci methods are not in general 
helpful when impulse responses are examined, since there 
doesn't appear to be any simple correlation between pole motion 
and the corresponding transient response of the system. The 
trouble arises because the pole -zero concepts essentially 
refer to linear systems. As soon as the poles are in motion, 
the usual relations between pole and zero positions and over- 
shoot for example, are no longer valid. What was done in that 
case was to try to use probable existing simple correlation 
between pole motion and transient response, which would ex- 
plain accidental "optimality" of the exponent 0.5. 



loo 



I 




The following cases were examined: 

(1) SOPOC and BPOC with weighting B=800, C=10, E=l. 

(2) SOPC and BPOC (B=800, C=10, E=l) with the same 
closed loop C.E. 

(3) SOPOC and BPOC with weighting B=164, C=0.001 and 
E=0.001. 



(4) SOPC and BPOC (B=164, C=0.001, E=0.001) with the 
same closed loop C.E. 

(5) SOPOC and BPOC with weighting D=3000, C=10, E=l. 

(6) SOPC and BPOC (D=3000, C=10, E=l) with the same 



closed loop C.E. 

To establish an approximating relationship of the form 
6s> - / the following simplifying assumptions, 

valid only for the given submarine, were made: 

(1) Average depth rate during application of impulse at 

FT tank = (-0.7) - (-0.8) ft/s. 

(2) Pitch rate = to 



The above values are averages derived from tests made on 
the given submarine by applying forces of about 70 tons in one 
of the auxiliary tanks. To saturate the forward planes, the 
forces have to be applied in the forward or the central (auxil- 
iary) trim tanks. 



Then: 



^ F 5 R ■ x< + FSft-X?, + fsc- so • 



h ^\ = 6^cott^ = F6F1-X, + FB1&-X2 + F6C-x^ + 



<^S2. = SOFl'Xi + S06-X^ 4- SOC- X3 1- 'SOD • 



“Detl X I = X2.-1 , X5= X^^• t oqci ">^ A 



^s\ •= 



FSR -t 



FS 5 + rsg 
tv / 






Xi 






- 1 j 50H + SOC J . I V l^sotb 4 Sod 



• X 2 . 



6^co^=l ( ^^ytt^FBS+F^ 



•X^ 



xa 

rv. 
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1. SOPOC and BPOC with weighting B=800, C=10, E=1 



From 


optimal controllers 


From program 8 


FSA 


= 


-0.20694 X 


A = 0.77984 


FSB 




-0 . 4043 


B = 0.45489 


FSC 


= 


15.321 


C = 0.10600 


FSD 


= 


58.51 


D = 0.49318 X lO*'^ 


FBA 


= 


0.3155 




FBB 


= 


1. 1073 




FBC 


= 


-32.57 




FBD 


=r 


-73. 85 




SOA 


= 


-0. 31623 




SOB 


= 


-1.665 




SOC 


= 


53.71 




SOD 




158. 71 




Then, 


using N = 175 





^ (6.68x10 ^ -t - 0.07)- X 2 

^Sl (-9.32x10“^ -t -0.758)- X 2 

"X (0.129-t + 0.685) -X 2 and with ^2 x -0.8 
X -0.103‘t*-0.55 

For values of t between 0 and 12 s. , (time duration of the 

pulse force plus 2 s. for which the submarine is considered to 

continue its path with the same characteristics ), ^ ^ i. varies 

between 0.55 and 1.79. The ratio r = 7 ! - will vary from 

1 

0.79 to 0.244. The average value of r will be 0.433. 

Also, 6si _ xxc5~‘ • t - 0.07 and 

-9.59X \o'^ • I -0.15S 
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t 


0 


2 


4 


6 


8 


10 


12 






-O.OSQ. 


-0.1^2 


-o2S5 


-o, '^G 


-OA\S 


-0.^6'^ 



is approximately -0.245. 

and for 



The average value of f 
From 6scoTr)=6s2. Jr 

<^SCOrT| — (^5^ • (^ I - -1. 2-^5 

- b>S2 • ^ 



The corresponding root locus is plotted in Fig. V-1. For a 
system with a dominant pair of complex roots, values of 
in the region of = 0.4-0. 6 would be suggested. In this 
case, the order was used in simulation runs of 

a stern only mode submarine with taking the values 0.3, 

0 . 6 , 0 . 8,1 , 1 . 2 . 

To compensate for the effect of the forward planes, the 
magnitude and the time duration of the applied force were 
reduced to about 35 tons and 6 sec. correspondingly. The 



results 


are shown in Fig. V-2 to 


between 


0.3-0. 6 = 0.22- 0.4) . 


0. 3 


1 - 1. 245K 0.6 


0.562 \ K \ 0.321 


Using 


K = f 9.-±lp}\ = (0.433) 


0.562 


(0.433)^ ^ 0.321 


0.69 


/x / 1.35. 



and 



(Oi| 



If in the the above calculations the extreme values of 



= 0.092 and = -0.55 are used, they result in 

1. 12 X </ 3. 52 

^51 I ^ 

If the values -r = -0.469 and €)+» = -1.79 were used, then 

asa 

0.524 ^ X 0.92 



Finally, if the 
using average values 
651 = -0.165- hsz 



value N=200 



for 






had been used, 
and r, in: 



it would result. 



r =0.390 
0.54 X ^ 1.13. 

2. SOPC and BPOC (B=800, C=10, E=l) with the same closed 
loop C.E. 



From BPOC 

FSA = -0 . 20694x10 
FSB = -0.4043 
FSC = 15.321 
FSD = 58.51 
FBA = 0 . 3155 
FBB = 1.1073 
FBC = -32.57 
FBD = -73.85 
From program 4 

SOA =-0.31653 
SOB = 1.996 
SOC = 64.31 
SOD =-5.076 



From program 6 

C.E. = S^+1. 0044S^+0 . 49083S^ 
+0.1126S+0. 49 366x10"^ 



From program 8 
A = 0.80618 
B = 0.48109 
C = 0. 1126 
D = 0.49365x^d■^ 



^05 



Again, using N=175 

^ -O.Ol>^2. 

6^1 ■= ("- O. \ 0?> >: t - o. 55 

662= 5.0^5x 

\'^\oo = 

(6s>/6s 7^ ■= o-'^'^ 



Then 

65CO'on = ■ ( \ - o.^^\ o 

- 6s 5 • rq 

From the corresponding root locus in Fig. V-5a, b and 
for ^-0.2-0.4,rr^ is found between 0.34 - 0.64 . So, 

0.34 1 - 0.841 • K 0.64 

0.784 ^ K 0.428 

0.288 ^ ^ ^ 

3. SOPOC and BPOC with weighting B=164, C=0.001, E=0.001 
From optimal controllers From program 8 



FSA 




0.24 


A = 4.0547 


FSB 


= 


0.3459 


B = 1.7757 


FSC 


= 


-4.2 


C = 0.34961 


FSD 


= 


370.9 


D = 0.015596 


FBA 


= 


0.9705 




FBB 


= 


2.240 




FBC 


= 


-62.6 




FBD 


= 


-173.37 




SOA 


= 


O 

r — 1 

1 




SOB 


= 


-5.9127 




SOC 


= 


193.24 




SOD 


= 


700.3158 





\06 



with N=175 



^' 2,. '^6 5 ^) • '%% 

65 , < 0 ,€I 2-1 4-1.249) -XI =-0.99-l-i (xx-o,^-) 



Vlau 



o.^3G \ 



o. ' ^'2- 




ou 



5.5G 



For ^ — 0.2-0. 4, the root locus (Fig. V-6a, b) gives 
= 0.035-0.092 

Then 

6 SCOTT) = is 2 +- 5.56is2-6s2l 

- 6s2 C \ - A.'bG ■ V"] 



0.035 < 


V ‘ 


4.56 < 


/ 0.092 


0.965 < 


\ -K ' 


4.56 ^ 


^ 0.908 


0.211 


> ^ 


\ 


;> 0.200 


0.85 


/ X 




0.89 



4. SOPC and BPOC (B=164, C=0.001, E=0.001) with the same 

closed loop C.E. 

From BPOC 

FSA = 0.24 
FSB = 0.34592 
FSC = -4.200 
FSD = 370.9 
FBA = 0.97051 
FBB = 2.240 
FBC =-62.6 
FBD = -173.37 



<07 



C.E. from program 6 



4.74253^ + 3.64433^ + 1.37133 + 0.015405 

From program 4 
30A =-0.98777 
30B = 668.35 
30C = 1857.2 
30D = -22229 . 8 
From program 8 
A = 4.543 
B = 3.6345 
C = 1.3713 
D = 0.015405 

Again, (!)S\ — + 2. AGO 

,^§1 ~ o.G\2-t -V 1.2^3 

Q-u - 0 . 1 G 2. 

Also, 9,G2.‘l +^91.^ 

and 

From the root locus (Fig. V-7a, b) and for ^ = 0. 2-0.4, 
iH = 0 . 12-0. 22 

0.22 
0.785 
0.133 






0. 12 

0.886 

0.066 



/ 1 - K 
K 
X 



0.99 



3 < 



10fl 



5. SOPOC and BPOC with weighting D=3000, C=10, E=l. 

From optimal controllers 

FSA = -0.030196 
FSB = -0.25864 
FSC = 24.688 
FSD = 67.0746 
FBA = 0. 31478 
FBB = 1. 13873 
FBC = -31.983 
FBD = -73.249 
SOA = -0. 31623 
SOB = -1.7386 
SOC = 57.59 
SOD = 165.97. 

From program 8 
A = 0.81664 
B = 0.49528 
C = 0.10839 
D = 0.49318E-02 

Using the above, 

= (o-Lll-t -to. 

6zi r (o.0^285•t -0.1'^ ') -X2 
= -0.1056*t -0.576 

Then 

'f'ao — 0.417 




1 ^Sl\ - - 1 . 16 



8scorr^ - t K ^ * 6 6 SX -6 S7') 



- 6sq-( I - 2.16 =^-h>S2 



From the root locus (Fig. V-8a, b) and for ^ = 0.2-0. 4, 



rn. ^ 0.3-0. 6 
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Then 



0.3 

. 0.324 

1.23 

6. SOPC and BPOC (D=3000, C=10, E=l) with the same closed 
loop C.E. 

From BPOC 

FSA = -0.030196 
FSB = -0.25864 
FSC = 24.688 
FSD = 67.0746 
FBA = 0. 31478 
FBB = 1.13873 
FBC = -31.983 
FBD = -73.249 

From program 6 

C.E. = + 1. 15723^ + 0. 655273^ + 0.178723 + 0.018267 

From program 4 

SOA = -1.1713 
SOB = -120.9616 
SOC = -148.4758 
SOD = 5606.32 

From program 8 

A = 0.95901 
B = 0.64553 
C = 0.17872 
D = 0.018267 



> 



1-2.16 K 
K 
X 



^ 0.6 
0.185 
1.84 



UO 



Again, A- O. 

- (^O.l'i^-'t V 0.1S6) ' X2. 

^ -O, \056-i -0.51G 
r"oo - 0.^1 "7 

Also, 6 s 2 ^ C' t 

and fv^au " 1 x lo ^ _> c^sco^ = Ss 2 n - 

^ - 6sa- rn . 

From the root locus (Fig. V-9b) and for 'T^ =0.2-0. 4, 
= 0.46-0.63 

0.536 ^ K ^ 0.367 




It was not expected that use of the previously, almost 
heuristically found values of exponents would result in each 
case in "best" response to impulse inputs. Only for the 
controller using weighting B=800, C=10, E=l, simulations 
were done for the stern only mode submarine to find the value 
of rn, iri the order fn-<^S2. that would result in "best" 
response. From that point, the condition = 0.2-0. 4 was 
arbitrarily used in the rest of the controllers, although from 
the corresponding root loci distinct differences in the posi- 
tions and movement of roots are apparent. 

If the results would approximate that of the simulation 
in Chapter 4, it could be said that the best value of the 
exponent is the one that corresponds to a stern only mode system 
with order rn.-<^S2 and taking the value resulting to the 
better response of a stern only planes submarine to an impulse 
input. The form of the order allows smooth switching, but the 
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only thing that determines the value of the exponent and can 
result in a different value each time, is the design of the 
separate controllers in CMC. 

For easy comparison, the ranges of X variation for each 
case of controller, derived in that section under the assump- 
tion of linear relation of the form 65 '= P\-<^SX , are incor- 
porated and shown encircled in Table IV- 10. Furthermore, the 
above analysis indicates that using the relation ^ 2 . ~ A 

in an inverse way, then: 



troller will not make use of a separate stern only mode design 
but there will be some difficulty in defining the values of ckr 
■8 . The parameter M^will depend on the specific design of the 
BPOC and the specific submarine. 

To demonstrate and check the above, root loci for 6 si- n 
were plotted for the three BPOC (Fig. V-10). 

For the controller B=800, C=10, E=l, i. is selected from 
values that it is expected to give roots satisfying the cri- 
terion of dominant complex pair. 




I n 




The above suggests that a controller having 

for K 1 could be equivalent to the proposed CMC. This con- 




5^^L0WER = 0.4 = PHl 
UPPER = 2.1 = 



Then 



or 



1 



n 



I 








U 



S I 



For (\ = -4 /Indicated from the ratio j — = -0.245 



I found before 



Ss2. 



(n^-rO.XS) ■ ^ ^ X V 0.^5) • % 

Using for K as before the average (o.433) , it is found: 

&o<^ o. ■eo<^ o.A5'b 



and 

0.781 ^ X ^ -0.755 

Thus, the control order to be tested will be 

^iCOrr, * 

Values of X between -1.0 and 0.75 were tested and the 
results are shown on Fig. V-11 to V-14. The submarine 
responses are acceptable. Numerical results are collectively 
presented in Table V-1. It is seen that the maximum depth 
overshoot with X = -1 is lower than any one achieved with 
the previously tested CMC schemes. The controller has the 
same degree of generality as the proposed one when the expo- 
nent X is fixed to a specific value. 

It is functioning effectively only with a specific system 
of submarine and controller and already failed when tested 
with 3P0C having B=164, C=0.001, E=0.001. 



B. EVALUATION OF CMC 

The conclusions concerning the proposed CMC are: 
a. It functions satisfactorily for certain designs of the 
separate BPOC and SOPC or SOPOC that compose it. 
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b. It is desirable that 6si in BPOC, if used in a SOPC 

as the order to the planes, will result in a stable submarine. 

c. In the cases of satisfactory operation there appears 

to exist a relation between the orders which also 

can result in a specific "optimal" value of X, not necessarily 
0.5. 



d. In any case, the CMC using BPOC/SOPOC or BPOC/SOPC 
(same C.E.) schemes may not be the near surface optimum one. 

As for example, when in one of the simulation runs (BPOC with 
B=800, C=10, E=1 and SOPC with the same C.E.) the value of 

in hs'i was mistakenly written 5.056 instead of the correct 
50.56, the resultant submarine response (Fig. V-15) was much 
improved compared with that corresponding to the correct value 
(Fig. IV-44) ! 

The last thing that was examined was the possibility that 
for a given submarine, use of optimal theory for a BPOC and 
SOPOC would always result in the same relation between , 

possibly a linear one, independently of the values of the 
weighting parameters. 

A combining function of the proposed form could then 
have elements of generality, at least for a family of 
controllers. This general relationship doesn't appear’ to 
exist when the optimal design is closely examined. 

The steady state solution of the matrix Riccati equation 
can be found in terms of the eigenvectors of the 2n x 2n 



partitioned matrix M = 






I Q 



The eigenvalues and corresponding eigenvectors of M are found: 






■ O 
I — 



I 






where the subscripts denote the signs of the real parts of the 
eigenvalues, and the corresponding eigenvectors are 









-T- 









is partitioned and the solution to the problem is 

K = V V"^ 

- -T+ -B+ 



The relationship between the values of K's in the two combined 

-1 T 

controllers is a function of B R B of the controllers, if 



A, Q are kept the same. 



(S 0 P 0 C has & - 
(B P 0 C has C> = 



o 

an 

O 
6 lU 






:) 



o o 
e»a 699 



ft- 



r 



C ol 
o c 



since this is the only term that changes in M. 

Now if Q and R are changed, resulting in another set of 
optimal controllers with new weighting factors, it becomes 
evident that only by chance could the same relationship be- 



tween S 0 P 0 C. and B P 0 C be found. 



X 


Maximum Depth 

overshootC ft) 


Maximum Pitch 

overshoot( rad ) 


Time to positive 
peak( s) 


- 




- 




0.5 


10. 


1 .496 


2 . 74x1 0'^ 


1 . 83x10"^ 


32 


-0.5 


8.2 


1.338 


3 . 34x1 0 '^ 


1 .54x10'*' 


32 


-0.75 


8. 


1 .42 


3.31x10"^ 


1 . 62x1 0~* 


32 


-1 . 


7.32 


1.57 


2.4x10"^ 


1 . 66 x 10 '^ 


32 



TABLE V-1 



Simulation results of the application of a 70 ton x 10 s 



pulse force at the forward trim tank. The submarine was 

" ■ \x 



'in trim* and the controller had Sscom =Si.i -<1.2.5 • 



0 .^ 3)6 









co'^n j 



and B=800. ,C=10. ,E= 1 . 



VI. DISCUSSION, CONCLUSIONS AND RECOMMENDATIONS 



The problem of optimum near surface control of a submarine 
is relatively new and not an easy one to solve. Much of the 
experimental data that would help in gaining a better under- . 
standing of this subject has been classified. Optimal 
theory has made the setting up of a large class of optimal 
control problems almost routine. Again, this is not the case 
here . 

Inequality constraints on the control variables suggested 
approximate methods and exact approaches for the solution, with 
which the author was not familiar. 

It becomes necessary in a number of important fields of 
engineering to devise controllers for plants which exhibit 
very substantial changes in their characteristics, due in most 
instances to a changing environment. Such a devised controller 
was tested and analyzed since design of the optimum near sur- 
face control system would require both the classified experi- 
mental data and advanced knowledge of optimal control theory. 

The proposed scheme of control order for use after the 
saturation of the forward planes introduces, in fact, an inten- 
tional nonlinearity, modifying the system characteristics. 

The theory of intentional introduction of nonlinearities into 
the system (linear model) is still decisively incomplete. A 
suitable approach is not evident for a straightforward analysis. 
A more extended examination could perhaps include describing 
function methods, phase plane description using a second order 
model or even consideration of the minimum time problem. 

AH 



As K appears in c^SCO'oq , it is a variable weighting factor 
that in a way selects corresponding parts of the <^S^ and 
orders, after the forward planes are saturated. Because of 

XT • J J • XT f ^ 

the way K is defined, i.e. K = 






COrq 



for 



/ its values of interest range between 1 and 0. 

The value of corresponds to either a value of 
infinite for small values of X, or to a large value of 

X ( 10^ - 10^) , depending on the value of the ratio 9 * . 

M co>g. 

For a stabilizing controller, c^-JcotT] after saturation, is 
expected to take values in a limited range. The exponent X 
directs or scales the value of K somewhere in the region (0,1). 

There is no indication nor even possibility that an expo- 
nent X=0.5 could automatically result in a value of K select- 
ing the right proportion of 6s-i and ^S2, orders, to obtain 
the best response for every given controlled submarine and to 
any inputs . 

The simulation runs show that even if the value of X=0.5 
results in an order (bscorr| which satisfies criteria of "opti- 
mality" in response to a certain input, changing the input 
will reveal disadvantages in other respects, which could be 
softened by another value of X. 

The submarine is already a nonlinear system. It may behave 
quite differently with different input functions. The pro- 
posed controller was tested for a limited number of step and 
impulse inputs at a constant velocity. The logical design of 
a nonlinear system requires a complete description of the 
input signals. 



As demonstrated in Ref. 6, the design of a linear system 
to minimize the mean square error with random input signals 
requires a knowledge of the second probability distribution 
functions of the signal and noise components of input. In 
contrast, the design to yield the nonlinear system which is 
optimum on the basis of the mean square error criterion 
theoretically requires knowledge of all probability 
distributions . 

To the most elegant design problem, involving the deter- 
mination of nonlinear elements appropriate for insertion in 
either a linear or nonlinear system to improve overall dynamic 
performance, perhaps the major difficulty lies in the problem 
of phrasing the specifications in a form amenable to analyti- 
cal design techniques. The marked dependence of the perform- 
ance characteristics of nonlinear systems on the particular 
input signals means that in general the specifications must 
include precise description of input and desired output. For 
the specifications developed during this work, the most prom- 
ising scheme of Combined Mode Controller appears to be the one 
that combines two separate optimal designs. It is then prefer- 
able that^sL also corresponds to the order of a stabilizing 
Stern Only Plane Controller. 

The exact tracking controller can also be used in the 
Combined Mode Controller. 

The Combined Mode Controller can be tested under a ficti- 
tious seaway simulated by the parallel application of step 
and sinusoidal forces. 



Finally as an alternate solution, a model reference 
method can be worked. In that case, an improved linear 
model is required. The model can be used in an stochasti 
estimator which will provide the necessary decision logic 
to modify the controller parameters. 



APPENDIX A 



ON LINEARIZATION 



In Appendix F the general equations of motion, given 
in terms of a set of axes fixed in the boat, are repeated 
from Ref. 1. These are six nonlinear scalar equations, 
representing the components of dynamic equilibrium in 
each of the six degrees of freedom. Also, six additional 
equations, relating the motions of the boat in body axes 
to the orientation and motion of the boat with respect to 
fixed axes are given. 

The above mentioned set of twelve nonlinear equations, 
characterize the submarine and are sufficient to determine 
its response to an arbitrary set of time dependent forces 
and moments. Where only vertical motions are of interest 
considerable simplification may be achieved. Although 
this constrains the applicability of the results, somehow 
the constraint will not be too severe, since submarines 
spend the vast majority of their underway time traveling 
on a straight course between the proverbial points A and B. 

Further, for the class of problems in which the sub- 
marine perturbations from equilibrium are small, the equa- 
tions of motion may be linearized and thus lead to even 
greater simplification with acceptable accuracy. 



Basically a Taylor's series expansion of all variables 
is made about an equilibrium condition. After substituting 



these expansions into the equation of motion, terms of 
higher than first order, are omitted. For small perturba- 
tions about the equilibrium condition, the actual equations 
are very closely approximated by the linearized equations. 

As the magnitude of the perturbations increases the accuracy 
of the linearized equations is degraded. 

The equilibrium condition used, will be for the sub- 
marine traveling horizontally at a steady speed and with a 
steady pitch angle (assumed zero for simplicity) . As part 
of the equilibrium condition the "in trim" condition of 
weight and longitudinal center of gravity is included. 

The initial reference values of all the variables are 
denoted by the subscript i and the small perturbations are 
indicated as follows : 

- changes in U, V, W etc. are indicated by small 
letters, i.e. 

U = u^+u, P = p^+p, etc. 

- changes in control surfaces are denoted by 




In this paper the submarine is considered as a rigid 



body and in Appendix B, derivation of the required equations 
from first principles is presented. (Ref. 9) 
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The kinematic and d 3 mamic equations of motion derived 



in Appendix are collected below for ready reference. 
The XZ plane of the coordinate system is taken at the 
plane of S57rnmetry of the boat so that 

^yz “ ^xy ^ 

Then 

X = M(U+QW-RV) 

Y = M(V+RU-PW) 

Z = M(W+PV-QU) 

( 1 ) 

K = I^P4QR(I^-ip 
M = IyQ+RP(I^-I^) 

N = I^R+PQ(Iy-I^) 

The details of linearization of equations for Z,M,X are 
carried through. 

The hydrod 3 mamic forces and moments depend on the 
orientation, configuration and motion variables. Thus: 

Z=Z(U,V,W,U,V,W,P,Q,R,P,Q,R,e,4>,H^,\3s,Or,Ob) . 
X=X(U,V,W,U,V,W,P,Q,R,P,Q,R,0,4^,S^,Os,Dr,Db). (2) 

M=M(U,V,W,U,V,W,P,Q,R,P,Q,R,0,q^,M^,Ds,Or,Ob) . 



For steady equilibrium the reference values of the variables 
are taken as constants. Making the appropriate substitutions 



Z - rr^l^co 4-(pvt V - (Q^l + C^') (UV +u"i] 

n = + (3) 

X = 4-('c^\ + cy^(^CAj; + co'^- 

where Z,M,X are given by: 



Z = Z Vivv^ . 0b+6b^ 

M =M(^ Uv>U,Vv-vV, Ob-^^b'l 

X = xCuv^u . - - ■ 'Ol^Sh) 



The Taylor expansion of Eqn. (4) is 



Z ='|^u(uOu +vDv (uOuvNiOv^-.c^bO^D ^h.o.t 






M =|^\ 4-(\-a0'ji V V Ov a- .. - ■ +^'o + (u.Du^vDv-(--- ':^b D<^.'^ 4 .u.o. i y\(utYl ■ 



X =|\ vfuOu. vOv+--- +^bDitoV (bXCXAvVPv^- SbD^b. ' ^ Q T^^vr. /..w.-.^b) 

• 



(5) 



Since we are dealing with small perturbations, from the 
equilibrium condition, the second and higher order terms 
may be omitted, leaving only the first two terms of each 
of the above series. 
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Writing the D operators in their partial derivative 
notation 



where the subscript i denotes that the partial derivative 
is taken at the initial equilibrium condition. 

The use of only first order terms, implies that the 
forces and moments vary linearly with the disturbance 
variables for small enough disturbances. 

On substituting Eqns (6) in Eqns (3) one obtains: 











(7) 
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Since Zj_, are equal to the unperturbed terms on 

the right hand side of Eqns. (7) respectively, they may be 
cancelled out. Retaining only first order terms and 
writing for etc. gives: 

2uVA+ZvV -^Zw -v-'iiZw f 2pp 4-ZoyCy + Zyr' + 

pZp Jr T Zy 4-©'Z© -V '-pZv + 6sZ^^ 4- S\^ 'Z-Sy, -hC^rZ^v — 

-u -Q^ui V pW-t-pvi'^ 

(Xu- VA V Y.;j V -I- vXvjU V- rr) \f{^j y x; _rv-| O/iT^ \j0 + X p p -h 

V XV Z fZ.pp + CXo^--op O0\*^p;^ ^ ^y'+'rp vi"^ V ^ Xo 0 + Xcj, t 

+ X^\3C>b + S s + X^v6v ^ 



( 8 ) 

^uU V V V M«.u-) 4 -'!Uu’U V MvV \ Mcu CO Jr M y,p lMo,cy t MyX + 
^ ^pp ^ y + 



i ^6 



For the straight, level symmetrical travel, these 



equations may be considerably simplified. 

For a truly symmetrical boat it is clear that the side 
force Y, rolling moment K, yawing moment N, and rudder 
moment are zero. Thus the derivatives of the assymmetric 
moments and forces with respect to the symmetric variables 
U, W, are all zero. 

Furthermore, v^, p^^, q^, equal to zero. Since 

we use stability axes =0 and is the reference boat 
speed V\.Q. Since we are considering level flight 0jL=O. In 
addition the following approximations are made: 

1. We may neglect the derivatives of the S37mmetric 
forces and moments X, Z, M with respect to the as 3 nnmetric 
variables v,p,r,(^ ,(£v' 

2. We may neglect all the acceleration derivatives 
except Zw, Zq, Mw, Mq. When the above simplifications are 
made, the linearization of Eqns. (1), for calm sea operation, 
results in 





-'M.u- 




U7 



Also 



q = 0 

(Depth) = cO - Oo 0 

The above equations are now compared with the corre- 
sponding Equations of Ref. 4 and those of Ref. 1, where 
linearization was accomplished by dropping all nonlinear 
terms. To make the comparison easier the equations are 
written in similar form. 

I . Derived 

^ -= 'ZojCO •^-'ZeOj-t-coqkc^' (orce:5 
- r + Mc^Cy -V- T Mg©~j fconhof focces 

II. Reference 1 

7^00 V vj co-^V^oe jovcei 

- ^ V C - )=L^ +- ^ Co© I Tof^ -fcNce^ 

III. Reference 4 

(^-7o^;cb - '7.0ijC\ - 'Zojco -t- C2.Q/+ + Co© ItoC ^ovcei 

-v<^cucb V ^co© Woftorcey 

Except for the fact that some equations make use of 
dimens ionalized and others of nondimens ionalized coefficients, 
the terms in the equations are almost identical. 

las 



The hydrostatic restoring term lAo© appearing in the 
set of equations in I and II above, exists because the 
introduction of a pitch angle G on a submerged submarine 
results in the creation of an hydrostatic moment , 
equal to 67 - in ft-tons, opposing that angle. 

This gives rise to an term in the pitch equation 

having a nonzero value. The ways that the existence of 
influences stability in the vertical plane are 
defined in Ref. 4. 

In the nonlinear equations of motion about the Body 
Axis System y-Axis in references 1 and 3, the term appears 
as Bz sin© and is neglected during the linearization pro- 
cess of reference 1, although it could have been included 
as Bz*0 , since small angles are considered. 

The terms Zg0 , Z^u, M^u do not appear in reference 1 
or in 3, from where the nonlinear equations were originally 
taken. Finally, the term muq (or mu^q) is neglected in the 
linearized model of reference 1. 

This thorough, yet not complete, investigation of the 
already derived and used linear model, was soon justified, 
when the realization in the form of simulation of the 
design of the stem only planes control system was attempted. 
Initially the model of reference 1 was used unchanged. 

After a few tests it became apparent that the model was 
insufficient for the present purpose. 
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The term "muq" was included in the equations without 
any improvement in the simulation results. 

Use of the linearized equation 



(Depth) = w - u-0 



instead of 



(Depth) = w which in fact was 

used when the original vector was rewritten in the 

, gave acceptable results in the respect 



state form 



A, 



>5 

X4 



of the linearized model mathematical representation and 
simulation results. 

Insertion of the term Bz-0 was not considered in order 
to keep the linear model as simple as possible and because 
a constant speed was used in all test runs. 

Details on changes due to the application of the 
described steps are given below for future reference, and 
for a constant forward speed of 15 knots' (25.33 ft/sec). 

1. As already mentioned, in reference 1, the problem 
was initially treated as a linear tracking problem. The 
state equations were; 



X = A X + B u, and the performance measure to be 



minimized was: 



T 







+ 



Va 



\a 



21 C-t') Q \ ^(^)~ w ^ 



d\ 



where V represents the ordered depth and pitch. 

The following approximations were also made for small 
perturbations : 

state variable form of equations 

0 0 



^1 


1 


0 


1 


0 


0 




^1 


X 2 




0 


^11 


0 


A21 


• 


^2 


x *3 




0 


0 


0 


1 




X3 






0 


Ai 2 


0 


CM 

CM 

< 




^4 

/ 



+ 



Bii ^21 



0 0 



®12 ®22 



where 



= Depth 

X2 = x^ = Rate of depth change 



x^ = Pitch 



X 



;^ = x^ = Rate of pitch change, 



1 



Using f = 



ll" ZoYd 
o 

©ovd 
. o 3 



for the command inputs then 



Xi - Zqj.^ = Depth error = 
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X2 = Depth rate error = E2 
X2 - 0 = Pitch error = 
x^ = Pitch rate error = E^ 



and the performance measure became 



'3='/'^ d t ( iov y_ = Q 



Also : 



0 



0 0 



0 0 



0 0 A2j_ 



0 



0 ^12 ^ ^22 



L ®4 



+ B 



Dt 



* -IT 

The optimal solution li = "R ® E E in the simulation 
program of reference 1 gave the equations: 



DEAD = 

DBAD = 5.70fQ + Ns}2 2 oDOT/m + 7^2- 47^rPiOO^ / ^ 



In the described model: 



= -4.3769>^Ao-2 

Pi pri^ 

^ 21 = ^ lAf E -lev + A tc- M - u • E »-v = -6.2278 




i^Z 



0.15443 



^22“ 



"ll= 



®12 



®21“ 



®22~ 



PiCi-2c^ A- 


. 'wl * 


-205 + RK 1 




r ACE'-Zos V Acc • ^ 


i • /f ’ 


h . C " 


1 


[Af£-'ZobA 


• Inv 


-Zob + 



and Inv = ACC-AEE - ACE -AEG = 6.3713 y 10^ 



= -0.42783 

= -9.4019^ 

= -0.24839 
= 2.0465 XACJ”^ 



2. If the term u*m*q is included in the equations of 
the linear model then 



A 2 ^= (AEE*(Zq + M) + AEC * ^ ^ = 7.4422 

A 22 = (ACE*(Zq + M) + ACC • ^ • U • INV = -0.1623 



3. If the term BZ*0 is included then setting BG 
(BZ=B*Zg), then 



_BZ 



2 



w 


' 


M-Zw 


0 


^Zq 


0 

1 




0 


1 


0 

1 


1 

q 1 


! 

c 


-Mw 

P 


0 


ly-Mq! 

J 



r 








U- Zw 

e 


0 


U(Zq+m) 




0 


0 


1 


1 


1 

1 UMw 

1 e 


BG 


• 





w 



0 

q 



1 55 



ds 



4 - 



M-Zw 


0 


-2-Zq 


-1 


a^.zds 


U^-Zdb 






e 


£ 


0 


1 


0 




0 


0 


-Mw 


0 


ly-Mq 




a^Mds 


ll^Mdb 

e 


\ 




/ 






/■ 



and finally 



y" .. 

• ^ 




w 




« 

0 




q ^ 





11 

0 

Ai2 



BG-^-Zq 

Determ 

0 

BG • 

Determ 



^21 

1 

'22 



w 

e 

q 



+ 



B 



B 



11 

0 

12 



B 



B 



21 

0 

22 



db 
< / 



ds 

db 



where Determ = (M-Zw) (ly-Mq) -Mw* Zq = 



= 1 

6.3713 X 10^ 

and ^ 21 ^ ^12’ ^22 » ^11 > ®21> ®12> ®22 ^^^e the same 

values as before. With BG —2.44 x 10^ then 



c( = ? i Zq ^ 0.0013 

Determ 

% = _ 0.0025 

Determ 



To study the effect of the changes in the values of the 
elements of matrix A, the eigenvalues of the plant matrix 
were found and are given at the end of this Appendix. 






To check for controllability the third order system 



(w, 0 , q) was decoupled using linear transformation. 



Kalman criterion was also 


used. In 


all cases 


was shown 


to be controllable. 




a. 


Original A matrix 








0 


1.0 


0 


0 




0 


-0.43769 E-01 


0 


-6.2278 




0 


o 

• 

o 


0 


1.0 




0 


0A1121 E-03 


0 


-0.15443 




Eigenvalues 








0.0, 0 


.0, -0.89659 E- 


01, -0.108539 


b. 


Matrix 


A after the addition of " 


-u-0" 




0 


1.0 


-25.33 


0 




0 


-0.43769 E-01 


0.0 


-6.2278 




0 


0.0 


0.0 


1.0 




0 


0.47727 E-03 


0.0 


-0.15443 



Eigenvalues 
Same as above. 



In the state model used in that paper the above matrix was 
substituted. 






c 



Original matrix A with terms "m«u«q" added; 






0 


1.0 


0 


0 


0 


-0.43769 E-01 


0 


7.4422 


0 


0.0 


0 


1.0 


0 


0.47727 E-03 


0 


-0.16123 










Eigenvalues 







0.0, 0.0, -0.18826 E-01, -0.18617 



d. Above matrix A with term *'-u»Q " included 



0 

0 

0 



1.0 




.25.33 


0 


-0.43769 


E-01 


0 


7.4422 


0.0 




0 


1.0 


0.47727 


E-03 


0 


-0.16123 



Eigenvalues 

0.0, 0.0, -0.18826 E-01, -0.18617 
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APPENDIX B 



DERIVATION OF EQUATIONS OF MOTION OF RIGID BOAT 
1 . Rigid Body Equation 

In this Section the submarine is considered as a 
rigid body and a derivation of its equations of motion 
from first principles is presented. 

From Newton's second law, the force F acting through 
the center of gravity c of the boat, in vector notation is 



where m is the mass of the boat and Uq its velocity. The 
moment of the boat G is equal to the time rate of change 
of its ngular momentum and is given by 



where the angular momentum h is , in terms of a particle of 
mass dm located a vector distance r from the center of 
gravity. 



dt 



( 1 ) 



= ^ 
dt 



( 2 ) 




dm 
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Since 



where the integration is taken over the whole boat. 

U is the same for all particles and J r d m = 0 the above 
equation becomes, after expanding the vector triple product, 

dm. 

Let r and cu be expressed by 



h = 



uo ( r • r) - r (uj . r 



r = xi + yj + zk 

^=Pi+QT+Rk 



(3) 



where X, y and z are the scalar components of r, P, Q and 
R are the scalar components of cO and i, j and k are unit 
vectors in the direction of x, y and z.^ Substituting 
these expressions in Equation 3 gives for the scalar 
components of h, in terms of the moments and products of 
inertia A, B , ... F , 



h 


= IP-I Q-I 


R 




X 


X xy ^ zx 






by 


~"^xy P ^y Q " ^yz 


R 


(4) 


^z 


~"^zx ^ “ ^yz Q '*■ ^z 


R 





^All rotations (control surfaces, angular velocity 
components, etc.) and moment components are positive 
according to the usual right handed convention; i.e. from 
the x-axis to the y-axis to the z-axis to the x-axis . 







■ 









' s 




■ 










where 



^x = 


ji/ 


+ 




dm. 


^yz 


/ yz 


dm. 


1) 

M 


/ 


+ 


x2) 


dm. 


^zx ~ ; 


^ zx 


dm. 


M 

N 

II 


7 (x2 


+ 


y^) 


dm. 


and = j 


^ xy 


dm. 



If , y and z are taken as the coordinates of dm in a 
non-rotating frame of reference, with origin at the center 
of gravity of the boat, C, then it is clear that in general 
the moments and products of inertia as well as the angular 
velocity components P, Q and R will vary with time as the 
boat rotates. This is an unnecessary complication which 
can be avoided if the coordinate system 0 x y z is fixed 
in the boat and allowed to rotate with it. Though this 
introduces additional terms the resulting equations are 
much simpler since now the inertia terms will remain con- 
stant. Thus Equation (1) becomes in terms of the velocity 
components of (U,V,W) and angular velocity uJ 

?=m^(ur+vJ + Wk) 

dt 

= m (U r+Vj + Wk) +(U^ + V^ + W^) 

I dt dt dt i 



since 



di 

dt 




61 X Uq 



(5) 



ojx i. 



t) t 



COX j 



and — 
dt 



X k. The 
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operator has the definition implied by Equation (5). 

Ot 

Similarly, Equation (2) becomes 



G — + oO X h. 

d t 



( 6 ) 



Combining Equations (5) and (6) with Equation (3) gives for 
the scalar components of F in the i, j and k directions 

= m (U + QW - RV) 

Fy = m (V + RU - PW) 

F^ = m (W + PV - QU) (7) 

and for the scalar components of & in the i, j and k 
directions 



II 

X 


+ Q 




- R 


hy 




M = hy 


+ R 


hx 


- P 


^Z 


(8) 


N = h^ 


+ P 


hy 


- Q 


hx 





Equations (7) and (8) are the Euler equations of motion of 
the submarine . 

3 . Motion of Boat Relative to Fixed Coordinate 

Solution of Equations (7) and (8) gives the linear 
velocity components U, V, W and angular velocity components 
P, Q, R relative to the 0 x y z axes fixed in the boat. To 



HO 



obtain the motion of the boat center of gravity it is 
necessary to express the linear velocities relative to a 
fixed coordinate system. 

We therefore define a fixed orthogonal coordinate 
system 0 x^ y^ in which the x^ y^ plane is fixed 
parallel to the equilibrium plane of the free water surface 
and the direction is positive downward. The orienta- 
tion of the boat axes (x,y,z) relative to the fixed axes 
(Xq, y^, Zq) is shown in Figure B-1. It is assumed that 
at first the two reference frames are parallel. Then the 
orientation of the boat is determined by considering the 
following three rotations, in the order indicated, where 
all rotations are in the positive direction. 

i. Rotate the Xq,Yq axes about Zq through the angle of 
yaw to the position (x]_,y^). Then the direction cosines 
between and (x]_ , y]_ , Z]_) are given by 





^1 


yi 




^0 


cos ^ 


-sin ^ 


0 


yo 


s in 


cos 


0 


^0 


0 


0 


1 



Thus a vector X]_ i]_ + y]_ + Z]_ k]_ in the X]_, yi, Z]_ 

system has the scalar components 



in matrix notation 



cos ^ sinM^ 

x^ sinM^ + cosH^ 

"i 

in the (x^jy^jZ^) system. Expressed 
these equations may be written as: 




^0 




cos V 


-sin 


0 




'xi 


^0 


= 


sin y 


cos 


0 




^1 


^o 


1 


0 


0 


1 

1 


1 


^1 



(9a) 



Thus Equation (9a) represents the first rotation. 

ii. Rotate the (x-|_,z^) axes about the y^ axis through 
an angle of pitch 0 to (x2,Y2j^2^ ’ may be expressed 

in matrix notation by 



^1 




cos 0 


0 


sin 0 




X2 




^1 


= 


0 


1 


0 




^2 


(9b) 


^1 

k 


i 


-sin 0 


0 


cos 0 




Z2 





iii. Rotate the about the y ^2 axis through 

an angle of roll <|) to (x,y,z) the actual orientation of the 
boat. This may be expressed by 



lAa 



1 



0 



0 



X2 




^2 


= 


^2 


i 1 



0 

0 



X 



coS(^ -sin<^ 

sincjp cos(p 



y 



(9c) 



z 



Combining Equations (9a) , (9b) , (9c) gives for the resultant 
of all three rotations 



X 

o 



^o 



z 



o . 



cosvp 

sinvp 

0 



-sin vp 
cosvl^ 
0 



0 




cos 9 


0 


sin 9 


1 


0 


0 


1 


jx 


0 




0 


1 


0 


0 


cos 


-sin 4? 




t 


1 




-sin 0 


0 


cos 9 


0 


sin($) 


cos4) 




z 



( 10 ) 



By carrying out the operations indicated by Equation (10) , 
the direction cosines between they,y,z axes and the 
y ,Zq axis are obtained. These are tabulated below. 





X 


y 


z 


^o 


cos 0 COSY 


1 

-cos<p sinvp 
+ sin 9 sin 4) cos'Y 


s in <p s in \p 
+ sin 0 cos<t> cosM^ 


3^0 


cos 0 sinf 


coscp cos\p 
+ sin 9 sin 4? sin 4? 


-sin 4? cos vp 
+ sin 0 cos<b 9invp 


^o 


- sin 9 


cos 9 sin<p 


cos 9 cos<|:) 



Resolving the U, V, W velocity components in the ^q’^o’^o 
directions gives for the velocity components of the center 
of gravity of the boat in fixed coordinates 



= U cos 6 cos + V (sin 0 sincip cos\^ - cos<^ sin\p ) 

+ W (sin4> sin\p + sin 0 coscj) cos\j;) 

dy^ 

— 2. = U cos 0 sin^p + V (cosc^ cosvj) + sin 0 sin 0 sin^jj) 

+ W (sin 0 coscp sinvjj - sinc^ cosvp>) 

= -U sin 0 + V cos 0 sin<^ + W cos 0 cos p (11) 

To obtain Xq Zq from integration of these equations, 
as a function of time, in the most general case would 
clearly require a very considerable effort. However, for 
most problems involving the motions of submarines these 
equations may be considerably simplified. 

4 . Angular Orientation of the Submarine 

The angular orientation of the submarine (0, ) may 

be expressed in terms of the angular velocity components 
(P,Q,R) by expressing the angular velocity of the boat in 
terms of the angles (0,'$^,V). Thus the resultant rotation 
of the boat from orientation (0,*^^ ,V) lo orientation 
(0+d0 ,<ii+d43,'^+d'v) in time dt, using the same order and axes 
of rotation as in the previous section, may be represented 
by 

d^ = (H^ + dy ) + (0+d0) + (d=+d4>) i' 

“ ^ - e 1 



where the subscripts indicate the axes along which the unit 



vectors I , j and k are taken and 



. j and i ' 
•^1 



1 as 



dt 



0 



dO=dVk+d0i +d<4>i 

o -"I 

and 

53 = 4^ "k + 9 K + 4>T - Pi + qT + R~k 

dt o 1 

Resolving k^ and © along the body axes x,y,z by the 
use of Equations (9b) , (9c) and (10) gives 

P=c^-\p sin© 

Q = 9 cos ($> + \^ cos 9 sin (13) 

R = ^ cos 9 cos 4^ - 0 sin<t> 

Solving Equation 13 for 9, and vp gives 
0 = Q cos - R sin<4> 

ci? = P+ Qsin<iD tan 9 + R cos cp tan 0 (14) 

V = (Q sin ct) + R cos ) sec Q. 
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APPENDIX C 



INVESTIGATION FOR VALUES OF WEIGHTING 
FACTORS IN THE COST FUNCTION 

A. STERN PLANES SUBMARINE 

In the design of the optimal controller for the stem 
planes only submarine, one of the problems faced with was 
the proper selection of weighting parameters that would 
result in the desired submarine response. At that point 
the idea of setting the submarine "in trim" had not been . 
considered yet. So, when the values C=10, D=3000, E=l, 
taken from Ref. 1, were tested, the resulting gains and 
finally the controller orders failed to keep the submarine 
at the ordered depth. 

A laborious trial and error procedure gave unacceptable 
results. The need for a more systematic search for the 
proper values of the weighting parameters led to the method 
described below. 

In short, this method uses linear feedback of states 
and the closed loop CE for root placement at selected 
positions. Then a linear system of equations is formed 
from corresponding relations in the optimal design and the 
simple state feedback design. Solution of the system 
gives values for the unknown factors. 



The problem is not solved exactly but the method helps 
to indicate magnitude relations between the factors and 
range of values from which they can be selected. 

The solution to the linear regulator problem gave the 
optimal control law 

u* = -r’^- K • E 
where R ^ = Cl = 1/C , or 



U* = -Cl • j^O B^ 0 B2J * 


^11 


K2I 


1 1 

ro 


^41 




Ki2 


K22 


^32 


K42 




Ki 3 


K23 


K33 


K43 1 

1 




Ki 4 


K24 


K34 


K44 



where for simplicity B^j^, B^2 were replaced by B^ . 

Then 

. r ' ' ' ~i 

U?'=-CI ' 1 B-]^ • K^2'^B2 * ^X4 ' B^ • K22"fB2 * ^24* ® 1 * ^32*^® 2 * ^34 1 ® 1 * ^42"^® 2 * ^44 * 

^21 ^31 ^4l] 

Deviating from optimal control methods, the problem of 
controlling the linearized model by use of linear feedback 
of states is examined. Ds is replaced by 



1^7 



SOA-E, + SOB-E^ + SOC-E^ + SOD-E, 
12 3 4 



where SOA, SOB, SOC, SOD are the feedback gains 



Then 





0 


1 


-u 


0 




0 


E= 


0 


Hi 


0 


^21 


•E + 


^1 




0 


0 


0 


1 




0 




0 


Ai2 


0 


Hi 




^2i 






|SOA • E]_+SOB • E 2 +SOC • E 3 +SOD • E^ 



or 



E = 



0 



SOA-B^ A^^+SOB-B^ 



-11 



0 



0 



0 



SOC-B^ A^^+SOD-B^ 



0 



S0A'B2 A^2‘^S0B*B2 S0C*B2 A22d’SOD*B2 



To find the CE the following matrix is formed 

^ U 0 



S 

-SOA-B. 

0 

-SOA-B^ 



-1 



S - A ^^- SOB - B ^ 
0 



-A^ 2 “SOB ♦ B 2 



■SOC*B^ -A 2 -[_-SOD-Bj_ 
S -1 



•S0C*B2 S-A22“S0D*B2 



Setting the determinant of the above equal to zero gives 



A 8 



HI* ‘.iMi 













+ (-A|^]_-A22-SOB*Bi-SOD-B2) -S' 



4- /A]_ 2. * ®2"^^22 * ® 1 

■Ai2-A2i-Ai2'SO 
-S0C*B2-S0A-B]_ 



“ A |^2 *^ 21~^12 ” ® 1“^21 * SOB • B2 



+ / -SOC • B]_ • A2_2"^'^X 1 * * ®2 

+SOA • Bj^ * A 22 ~^2 1 * SOA ♦ B 2 
+H.SOA.B2 



*S 



+ (U-SOA-B1'Ai2-U-SOA-B2*Aii) 
= 0 

where 

A2_2_ = -0.43769 E-01 
Ai2 = 0.47727 E-03 
A 21 = -6.2278 
A 22 = -0.15443 
B2_ = -0.42783 

B 2 = -0.94019 E-02 



The original positions of the plant eigenvalues are con- 
sidered and new positions selected according to some 
specified criterion (for example selection of a pair of 
relative dominant complex roots with negative real parts) . 
The new CE is formed and comparison with (1) gives the gains 
SOA, SOB, SOC , SOD (program 4). 






for roots: 



= -0.065 
S 2 = -0.51 

83^4 = -0.00102 + jO.486 

the CE is + 0.57524 S^+ 0.2702 + 0.13533 S + 0.77979 E-02 

Then SOA = -0.5 
SOB = -7.03 
SOC =32.5 
SOD = 360.0 

The order DSAD = -0 . 5 -ZOER-7 . 03 • ZODOT+32 . 5 • PERR+360 • PIDOT 
gave the submarine the satis fac'tory behavior shown in 
Figure C-1, for zero depth and pitch ordered. 

Another test was made, 
for roots 

s^ = -0.082 

S 2 = -1.998 

83^4 = -0. 0044 7+j 0.0 642 

and CE 8^2.08976 S^+0. 18699 S^ + 0.0010087 S+0. 679973 E-03 
the feedback gains were: 

SOA = -0.0436 
SOB = -3.49 
SOC = 0.0523 
SOD = 360.0 
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It has to be noticed here, that even if the values for s^, 
S 2 , S 3 , s^ in the examined cases appear to contain some 
magic relative to their selection, in fact they resulted 
from a limited trial and error circular procedure. The 
path was from roots to the CE, to the feedback gains, to 
simulation, to modification of the gains to improve con- 
troller performance, to new corresponding CE and finally to 
the given values of roots . 

Figure C-2 corresponds to the use of a feedback con- 
troller with the gains found in the second test. The 
results are acceptable. 

Finally for roots s-|^=0.1, S 2 =~ 1 . 0 , S 3 ^=-0.1 + jO.015 
and CE S^+1.3 + 0.33 +0.031 + 0.0010225, the gains 

FBA= -0.0656, FBB= -1.47, FBC= 26.2, FBD= 184.1 were found 
and tested. (Fig. C-3) 

The main reason for these repeated tests was to check 
how well assumptions concerning stability of the linear 
model were valid on the simulated nonlinear submarine.. 

During runs at the speed of 15 knots with both planes 
available, it was found that for a submarine "in trim" the 
following quantities representing ballast, had to be added 
in the three trim tanks . ^ AT= -6.4 x 10*'^, AU= 8.8 x lO”'^, 
FT= 6.4 X 10 ”^. 

^To give weight in pounds the above numbers are 
multiplied by (Length of Submarine)^. 
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Figures C-4 to C -6 show the behavior of the submarine 
for each of the examined control laws as before but with 
the required ballast added. Judging from the SoS„ error 
values the performance of all controllers has now been 
improved. 

The next step was an effort to correlate the results 
from classic feedback control with that of the optimal 
theory. 

For the feedback gains 
SOA = -0.5 
SOB = -7.03 
SOC =32.5 
SOD = 360.5 

the control order was: 

DSAD= -0 . 5 • E^-7 . 03 • E2+32 . 5,* E^+360 -E^ 

From the optimal law 

DSAD* = -Cl *|^F E]_+F 2 ]_ * E 2 +F 3 * E 3 +F^]_ • E^J 

where use of the definitions in main program of Ref. 1 is 
made . 

Comparison of the two orders show the following 
correspondence : 

-Cl • F^^^ = -0.5 = SOA 

-Cl ♦ F2 ;l " “7.03 = SOB 
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-Cl • = 32.5 = SOC 

-Cl ' F^i.= 360.0 = SOD 

This gives four equations, linear, algebraic, in K's, i.e. 

’^14~ “C*S 0 A = C • 0.5 

F 21 = B;j_*K22+B2-K24^ -C-SOB = C • 7.03 

F 31 = B]_*K23+B2*K34= -C-SOC =-C • 32.5 

F 41 = B^-K24+B2'K^4= -C'SOD =-C « 360.0 

Also from the nonlinear algebraic equations in matrix 
form 

K = 0 = -K'A-a'^-K - 2. + K'B*R"^.b'^«K 

the following ten equations written explicitly are used: 
(F 3 _P^ • CI - E = 0 

(F 11 .F 21 ) •CI-(Ki;l'^3_2-Ai1+Ki4'Ai2) =0 
(Fii'F3i) 'CI + U • Kii = 0 
(Fii'F 4 i) ' CI-(K3^2 •'^21+^13+^14'^22> 

(F 21 ) ^ ' CI-2 ' (Ki2+A;l1 ‘^22+^12 '^ 24 ) "A=0 
(F21 'F3]_) • CI-(K]^3+A]_]_ 'K23+K34 ' A]^2) ■*■0 'K2]_=0 
(F 2 i*F 4 ]_) 'CI-(K]^4+A]_]^'K24+K44'A]_2+A2i*K22+K.23+A22 ‘K 24 ) =0 
(F3i)2-CI+2'U'K2_3 - D = 0 
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(F3i*F4i) *CI-(A2]^*K32+K33+A22 •K34)+U'K]^^=0 

(F4i)2-CI-2-(K24-A2i+K34+K44‘A22)-B =0 

Substituting in the above: 

Fii = -C-SOA 
F23_ = -C-SOB 
F32^ = -C • SOC 
F 4 ]^ = -C-SOD 

and defining also: 
x(l) = 

x( 2 ) = K^2 
x(3) = 
x(4) = K ]^4 
x(5) = K 22 
x( 6 ) = K 23 
x(7) = K 24 
x( 8 ) = K 33 
x(9) = K 34 
x( 10 ) = K 44 



x(ll) = E (Depth error weighting factor) 
x(12) = D (Pitch error weighting factor) 
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x(13) = A (Depth error rate weighting factor) 
x(l4) = B (Pitch error rate weighting factor) 

A system of 14 linear algebraic equations is formed, from 
which given the specific feedback control gains SOA, SOB, 
SOC, SOD, the steady state constant K values and the 
corresponding values for the weighting factors in ^ of an 
equivalent optimal controller can be found. The system 
written in the matrix form A X = B is ; 
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0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 


0 




(FBA) 


1 


AA 


0 


AB 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 




FBA*FBB*C 


-u 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 




FBA*FBC*C 


0 


AC 


1 


AD 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 




FBA*FBD*C 


0 


3 


0 


0 


m 


0 


2AB 


0 


0 


0 


0 


0 


1 


0 




(FBB) ^*C 


0 


-U 


1 


0 


0 


AA 


0 


0 


AB 


0 


0 


0 


0 


0 




FBB*FBC*C 


0 


0 


0 


1 


AC 


1 


(AA-+AD) 


0 


0 


AB 


0 


0 


0 


0 


•x= 


FBB*FBD*C 


0 


0 


-2U 


0 


0 


0 


0 


0 


0 


0 


0 


1 


0 


0 




(FBC) 


0 


0 


0 


-u 


0 


AC 


0 


1 


AD 


0 


0 


0 


0 


0 




FBC*FBD*C 


0 


0 


0 


0 


0 


0 


2AC 


0 


2 


2AD 


0 


0 


0 


1 




(FBD) 2*C 


0 


BA 


0 


BB 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 




-FBA*C 


0 


0 


0 


0 


BA 


0 


BB 


0 


0 


0 


0 


0 


0 


0 




-FBB*C 


0 


0 


0 


0 


0 


BA 


0 


0 


BB 


0 


0 


0 


0 


0 




-FBC^'^C 


1 

: 0 


0 


0 


0 


0 


0 


BA 


0 


0 


BB 


0 


0 


0 


0 




-FBD*C 

1 


c, 


the control 


effort weighting parameter. 


is con 


sidered 



known in the above. Its value is determined by the designer. 



For the gains in (1) and C=l, solution of the system 
gave E=0.25, D=-0.6686 x 10^, A=-0.1137 x 10^, B=0.233 x 10^ 
The negative values of A and D are not acceptable in the 
optimal analysis. So C was made equal to 0.001 and as a 
result- of the form of equations (right hand side multiplied 
by C) , the above values were changed to E=0. 00025, D=-6.686, 
A=-0.1137, B=233. Having reduced their values, A and D were 
then neglected and the optimal method tested with values of 
E, B, C around the found ones. 

The values of C=-.001, E=0.001, B=164 were selected for 
subsequent use and results of their use are shown in Fig. 

C-7 and C-8 with the submarine "in trim" and "out of trim," 

The S.S. errors are comparable with that of the feed- 
back controllers and acceptable. 

As for a first indication of the behaviour of the con- 
trollers and submarine under a simple harmonic seaway, the 
sinusoidal forcing function: 

AU = 2 • E-05*sin( t) 

was introduced. Results are presented in Figures C-9 to 
C-12. 

The described method has the disadvantage of the possibl 
negative weighting factors that can be found. Variation of 
the value of C can be used to lessen the magnitude of these 
negative values and indicate a region of the positive ones. 



A good linear model will improve the results and allow 
use of the CE (1) for a systematic search for the "best" 
placement of roots either by parameter plane or root locus 
methods. Using the described method the controller C=10, 
B=800, E=1 was also found, corresponding to a closed loop 

CE 5^+0.9780375^+0. 4646185^+0. 106002S+0.493183E-02 and 
roots 

s^ = -0.06077 
S2 = 0.40235 

S3^4 = -0.25745+j0. 36798. 

Results of its use are shown in Figure C-10 with the 
sinusoidal input at AU. This controller, having acceptabl 
performance, was also used in the formulation of the CMC. 

B. BOTH PLANE5 5UBMARINE 

Extension of the method developed in the case of the 
both planes submarine was proved very tedious and of doubt 
ful usefulness. 

For this, only a few results will be presented that 
would be useful in similar studies and the intermediate 
steps will be omitted. 

Use of 

D5AD = Yi3_-E]^+Y2i/U*E2+Y3i-E3+Y43_/U-E4 
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DEAD = 

in E = A E + B U, results in the closed loop CE given 
below. 

"^11 » ^21, Y31, Y4;l> ^12, Y22, Y32, Y42 were defined in 
Ref. ( 1 ) and the following substitutions were also made 
for the general case of feedback controller. 

FSA = Yj_j_ FBA = Y ^2 

FSB = Y21/U FBB = Y22/U 

FSC = Y^i FBC = Y32 

FSD = Y41/U = "^42/^ 

Then 

^ • FSB - 622 ^ * FBB -A 22 “®22 * 



"^A 2^2, • ^22”^® 11 * ^22 * FSB+B2]_ ’ ^^2 * FBB+A2_]_ • B^^2 * ^ 

S^W ”I"B 2 ]_ * ^ 2 ^ 2 * ® 22 * 1 1 * ® 22 * 

- A21 • A ]_2 ”‘^12 * ® 11 * ”^12 * ® 2 1 * FBD - A2 2 _ * B 2_2 * 

“B21 • B 2_2 * FBD • FSB“A22_ * ®22 *BBB-B2_2^ * ®22 * * FBB 

-Bii • FSA-B2 2_ * FBA- B 2^2 * ”®22 * 



+ S 



^^22 * ® 11 * ^^^”^®22 * ® 11 * BSA 4 A 22 * ® 2 1 * 

+B12 ' FSD • B22_ * FBA-A2 2_ • B2_2 * FSA-B22_ * FBD * B 2^2 * FSA 
-A22_ * ®22 * FBA-B 2_2_ * FSD • B22 * FBA4A2_2_ * B 2_2 * FSC 
+B2 1 * FBB * B 2^2 * FSC4A2_2^ * ® 22 * FBC+B 2_2_*FSB*B22* FBC 
-A 2^ 2 *B2_2^*FSC-B22 * FBB * B2_2_ * FSC -A2_2 * ® 2 1 * FBC 
-Bi2*FSB-B2i-FBC+U-B2_2-FSA+U-B22'FBA 
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^B22 * FBC • B]_]_ • FSA+B]_2 • FSC • B2i • FBA-B21 • FBC • B]_2 • FSA^ 
-B]_]_ • FSC • B 22 • FBA+U- A]_2 ‘ B]_]_ • FSA-U • A]_]_ • B]_2 • FSA 
+U«B22*FBB*B]_]_*FSA-U*B2i*FBB*B]_2*FSA 

+U • A]_2 • B2 1 * FBA-U • A]_]_ . B^ 2,’ 

\+U-B]_2*FSB'B2i-FBA-U*B]_]_*FSB-B22*FBA 



= 0 

The simplest case will be to accept the relations: 

FBA = FSA 
FBB = FSB 
FBC = FSC 
FBD = FSD 

Then the resultant CE will be of the form 
+ (bi-FSB + di-FSD + el) -S^ 

+ (a2‘FSA + b2'FSB + C2*FSC + d2*FSD + e2)‘S^ 

+ (a3*FSA + b3*FSB + C3*FSC + e3> 

+ (a4*FSA) 

Parameter plane methods will help now in selecting positions 
for root placement and give values for gains. The next in 
difficulty but more interesting case will be to accept 
relations of the form 



I 



FBA = -K-FSA 



FBB = -K-FSB 
FBC = -K-FSC 
FBD = -K-FSD 

The characteristic equation can result, by the applica- 
tion of some stability test, into conditions on FSA, FSB, 
FSC, FSD for stable roots in the form of a system of 
inequalities or again parameter plane methods can be used 
to give values of the feedback gains. 

Finally desired values can be given to four of the 
feedback gains and the other four used as parameters. 

To correlate the results with that of optimal control 
analysis on altered form of the DSAD, DBAD equations is 
used, i.e. 




E 



r 



-1/Cj^(Fii-Ei+F21-E2+F3i-E3+F4i-E4) | 

"1/C2(Fi 2 ' El+^ 22 * ^2+^ 32 ' E3+F42 * E4) ^ 



Then 



FSA = -1/C;l*Fii 
FSB = -1/C;l*^21 



FBA = -1/C2*F;l2 
FBB = -1/C2*F22 
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FBC = -1/C2-F32 



FSC = -l/C^L'^Sl 

FSD = -1/C^*F4]_ FBD = -1/C-2-F42 

With linear feedback gains selected for a desired CE 
and using the defining equations o ^ ^11“®11’^12‘^®12'%4’ 
Fi2=B2i ’K 12 +B 22 '^14 > etc. together with the Riccati 
algebraic equations for the steady state solution, a system 
of eighteen linear equations in sixteen unknowns (C]_, C 2 , 

A, B, e, D, K 22 , K 33 , K 44 , Kx2, Kx 3 , K 23 , K 14 , K 24 , 

K 34 ) results. 

This requires the introduction of two additional non- 
diagonal weighting factors as unknowns in the 2 . matrix and 
consequently in the algebraic equations. Solution of this 
system of eighteen linear equations will give the corre- 
sponding optimal control weighting factors, where once 
again the possibility of negative values is not excluded. 
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APPENDIX D 



SIMPLIFIED TRIM ANALYSIS 



As the submarine changes speed, hydrodynamic forces 
and moments change in proportion to velocity squared. 

When the submarine is at a neutral bubble, these changes 
are in equilibrium and is neither necessary to change the 
weight of the submarine nor to shift the center of gravity 
Necessary changes of variable ballast are accomplished 
with the trim system. In this paper three trim tanks are 
used and shown in Fig. IV -4. 

To determine the required ballast at the speed of 
15 knots, (the constant speed of all simulation runs), 
the following approximate method was used. 

The equation of motion along the body axis system 
z-axis when w, v, p, q, r and their derivatives are set 
equal to zero gives : 



W. 




b cos 0 • cos 



W 0 =^>ballast to be added 





blown 



where 






For 0 = dp =0, Zuu = -0.1x10"^ (Ref. 3), 2 = 415, ^/2 
Wb = 17.22 • 



1 



This parabola is plotted in Fig. D-1. 

For U'= 25.33 ft/s = 11050 Ibf. 

From the equation of motion about the body axis system 
y-axis 



X = B-Zg‘sin 0+P/2*^^-Muu-U^ 
W‘cos 0 



where 



W 

B 

Zb 

For 0 



^ = 



= X coordinate of the CG 

= m*g = weight of submarine = 0 . 625xl0^x32=2xl0^1bf 
= buoyant force 

= separation of submarines c.b and c.m 



= 0, Muu = 0.4 X 10”^ (Ref. 

e^.Muu.U^ - . 1.43 K 10-4-U2 



w 



This curve is plotted in Fig. D-2 
For U = 25.33 XG X 0.1 



The above calculations suggest addition of 11050 lbs in 
the auxilliary trim tank to compensate for Wg . 
ballast has to be added since X^ = 0.1 0. 



Additional 



For FT = -AT 



W-Xg = FT* (175.5+215.5) 

FT = 

391 

Setting W = 2x10^ + 11.05x10^ 

FT 5118 Ibf and AT -5118 Ibf. 

If these values were used in the program, they would appear 
as AT=-7. 16x10-5, FT=7 . 16xl0"5 , AU=1 . 546x10"^ 

The values used were 

AT=-6.4xlO"^, FT=6.4xlO“^, AU=8.8xl0'^. 






iK r ”T’ KT fifli" 



iimii 



Kw- 

^ Efc ‘MMT 





APPENDIX E 



In an effort to find limits on the exponent x in 

K » (0^)^ , 

O fcom 

different computer programs were developed. 

The programs use almost the same approximations but 
different criteria each time lead to different regions 
of variation of x. 

A short description of each one is given next and 
intermediate steps are omitted where the analysis is 
straightforward. Because of the approximations used it 
was not expected that the results would be much restrictive, 
even sometimes acceptable. Further refinements are possible. 
In all cases the root loci analysis already described in 
Chapter V gave better results. 

1. The first two programs (programs 9 and 10) result from 
an analysis of the response to a step input of the sub- 
marine with the CMC and optimal design in the combined 
controllers. The criterion used was minimization of a 
function of steady state depth and pitch errors. 
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From the designed optimal controllers 



= FSA*Xi+FSB*X 2+FSC*X3+FSD*X4 | order 

^ s^ = SOA*Xt+SOB*X^+SOC*X^+SOD*X, I order 

^2 1 2 3 4|in SOPOC 

S = FBA*Xj+FBB*X 2 +FBC*X 3 +FBD*X 4 ^^^'^^^gp^^^ order 



First approximation 



X^ = X^ * t 



X 3 = * t 



Substituting in ^ and solving for t results in: 
S £1-FBB‘'-X2-FBD-'-X4 



t = 



FBA*X2+FBC*X4 



Second approximation 



X£ = P * X 4 (P experimentally found between 150-200 

for the given submarine and acceptable 
response) 

Then 

t = ^ fl-(FBB+FBD/P)-^X 2 

(FBA+FBC/P)*X 2 




FSA+FSC/P 

FBA+FBC/P^ 



+X2' 



( -FSA-FSC/P)^'(FBB+FBD/P) 
FBA+FBC/p 



+FSB+FSD/P 

J 
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h S2=hf-\* 



SOA+SOC/P 




^FBA+FBC/P. 


L 



(-SOA- J O | Z | ^ . ( ^ ^ ^ ^^ +SOB+SOD/P 



J 



Since 



^ scorn ~ ^ s? ^ ((Ss]^-(^S 2 ) then at the steady state 



scorn O s2 
with X 2 = 0 

scom~ <^f 1 



' ^SOA+SOC/pA +j.,v f FSA+FSC/P-SOA-SOC/p \ 
FBA+FBC/P/ i FBA+FBC/P ) > 



Using the limiting condition 



-0.436 



c) som 



0.436 



and setting 



H 






fi' 



FSA+FSC/P-SOA-SOC/P 

FBA+FBC/P 



Then for H\0 and Q-4.36 . / 1 

/ |6fcom^ 



log 



'-0.436-^fl*( S g| +S 0 C | P )' 



H 



log 

[ofcom 



I'm Ijr /SOA+SOC /P^ 

0.436-df]L''(pg^^Pg2yp) 



>> 



log 



H 



log ( 



0.436 

[ c3fcom 



where fcom = C^fl 



For H <^0 the signs of the inequality are inversed. 
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In order for the above unequalities to be of some use, 
the value of to be substituted in, must be specified. 

From the linear model 







0 


1 


-U 


0 




1 — 1 

X^ 




0 


0 




0 


X2 




0 


^11 


0 


1 — 1 
CM 

c 




X2 




hi 


®21 


U+ 


F 


• 












0 




+ 




X3 




0 


0 


0 


1 




X3 




0 


0 




1 0 


• XI 




0 


Ai2 


0 


^22 


1 




1 


®12 


®22 

i 


1 

i 

1 ! 


M 

1 

1 



where F, M are terms due to external step forces and moments, 
and with U = + (l-K)* ^ ^2 

= ^K*FSA+(I-K)*SOA^ *X^ + 

J^k*fsb+(i-k)*sob] *X2 + 

[k*fsc+(i-k)*soc^ *X3 + 

[k*fsd+(i-k)*sod'] *X^ 

it is finally found 

s.s depth error = K“Rl+R2 

K*R3+R4 

s.s pitch error = 

K*R3+R4 

R^^, R2, R3, R4, R3, Rg are given explicitly in programs 9 
and 10. 
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The function to be minimized is 



2 

(s.s depth error) + (p * s.s pitch error) 



Taking derivatives with respect to K gives K as a solution 
of the quadratic equation 



K^+K* 



[riR4+RiR2R3-2P^R3R5R61 + 

[R?R3-p^R3R5] 



[R1R2R4- p^RaRe] = q 

lj^R^R3-p2R3R^ "I 



( 1 ) 



At steady state 



6 focm = + FBC 



It was also found that 



( 2 ) 



^Iss 



R^^K + R2 

R3K + R4 



^2ss 



= ^5^ + Re 

R 3 K + R 4 



(3) 



2 2 

and that the value of K minimizing X^^gg + (P*X 3 gg) is 

the solution of (1). Substituting this value in (3) and 
using (3) in (2) results in the desired value in the x 
inequalities to give the region of x variation. This 
region is expected to be reduced if the values +0.436 in 
the limiting condition are lowered depending on the magni 
tude of the applied moments and forces. Also, it is 
possible instead of using in the above the value of K 
minimizing the function of depth and pitch errors, to use 
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a second order approximating system and express the steady 
state errors as functions of specified (for example 
0.4-0. 6), settling time (for example 120 sec) and input 
forces (for example 70 tons). These values of Xj_gg, ^ 3 gg 

will subsequently give the value of ^fiss used in 

the X inequalities. 

2. The application of step inputs to the system leads to 
the saturation of the forward planes and a different from 
zero angle of the stem planes, before the CMC is energized. 

For a more elaborate analysis the effect of the above 
must be considered. This is done in programs 11 to 13. 

In program 11 the effect of saturated forward planes is 
translated to initial depth and pitch errors. 

Programs 12 and 13 are almost the same as 9 and 10 but 
properly modified to accept results of program 11. 

3. The next approach was to put limits on the exponent 
of K considering the response of the system to impulse 
inputs . 

The starting point was to consider a second order 
approximating system, find the closed loop CE with 
U=K,S^^+(1-K)^ g 2 > accept a value of and use the form 



no 



2 T 2 

S +2 Wj^S-i-W^ and the relations between the power coeffi- 
cients to get a corresponding value of K suitable for use 
with the inequalities from the step response. (Programs 
9-13) 

In program 14 the found CE of the form S^+A(K) S+B(K) 

2 T 2 

is set equal to S +2 \ ^ relation derived which 

gives values of K with ^ as a parameter. 

Finally in program 15 further approximations are made 
to result in an expression giving K with parameters P(150-250) , 
^ , and X 2 (average value ±0.8 depending on side of response) 
and fcom independent variable . Using this program, 

the CMC with separate optimal controllers having weighting 
factors B=800, C=10, E=l, also - 0. 6-0.8, P=220-250 
and X 2 - -0.2, values of K varying as function of fcom 
were found. The results are shown in Fig. E-1 to 6 where 
also ^ 0.436 >^ 

c5'fcom is plotted for different values of x. 

In these figures the region of comparison is for 
between 0.74 and 0.9. The curve with X=1.2 is the one 
that better approximates the K-curve found by the program. 
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4. Additional informations for the use of the programs. 

a. When forces or moments are desired to be introduced 

in the right side of the linear model state equations, the 

3 5 

quantities already divided by either L or L have then 
to be multiplied by the matrix 



XMA 


XMC 


I ^ 


61.72 


-5288.51 


XMB 


X^® J 




-0.0031 


1035.02, 



For example application of the step -2.208x10”^ at the FT 
results in the terms 



61.7a 


-5288.51 


1 

1 


-2.208x10"^ ^ 


J 


- n 

FRC 


-0.0031 


1035.02 J 


1 


-2.208x10 ^ • Xpij 




TRQ 






L- 4152 


. 


L__ ^ 



b. 



SCA, 


SCB, 


see, 


seD 


are 


the 


same 


as 


FSA 


, FSB, Fse, 


FSD 


BCA, 


BCB, 


Bee, 


BeD 


are 


the 


same 


as 


FBA 


, FBB, FBe, 


FBD 


^11’ 


^21’ 


^12’ 


A22 


are 


the 


same 


as 


AA, 


Ae, AB, AD 




^11» 


^21> 


^12» 


B22 


are 


the 


same 


as 


BA, 


Be, BB, BD 
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APPENDIX F 



EQUATIONS OF MOTION 



The following set of equations are referred to a body 
fixed systen of axes which are coincident with the prln 
clpal axes of Inertia of the body. The origin of this 
axis-system is located at the assumed center of nasa of 
the bodv 



Equation of Motion Along the Body Axis System x-Axls 



m(u - vr + wq) r ^ 




* + X * + X 

^ rr 





u + X * vr + X 
vr 



wq 





* + X * + X 

vv 



ww 





+ I X 



vvn* 



(n> - 1) 




wwn 



(n' - I ) 



+ X 



Ssdsn' 



' (a' - I) 




6 r6 ra’ 



' (n' - 1) 6^^ 



sin 9 




1T3 



Equa-tlon of Motion Along the Body Axis System y-AxLa 



m(v - wp 



ar) .ft- [r; ■ f * Y . ■(,] 



ft- [y 



pq 



pq + Yp|p| > pIpI] 



+ [y. *v+ Y^ • wp+ Y,|,|,V||(v^.w*)h Irl] 






[""r 'ulrUr. Yp -up] 



+ 1^' ^rn- ’ 






+ Y . '6 

2 6r r 



+ -f Y , , ’ (n’ - 1) 4 

2 6ra* ' ' 3 



+ 1 Y , ' (n' - 1) uv 

2 vn* ' 



+ 2^'^v|vU' ’ 



+ -^ Y * wv# 

2 wv 



+ (F ) 

2 y V3 



(-w) Sin, cot 



U 



+ D W • sin (|) cos 0 



? Multiplied by 



U 

for large angles of 
attack near -90^ 
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Equation of Motion Along the Body Axis System z-Axls 



m(w - uq + vp) = 















+ 



+ 



+ 



+ 






+ 



+ 



^ t* Z . • q 
2 q ^ 





fz 

L rr 


' rM + Z 

rp 


' rp#] 

Noce 1 




0 


[^W ’ 


w + Z • vri!l + Z ' vp + iZ ’ vp#1 

vr vp ^ vp J 




- 

z 


[% ' 


"'1*^ |,|6. 


• ujql 6s f Z I ,'ifj(v* + 
wlq] Iwf 


i»i 111] 




^ 1 ’ 
qn' 


' (n' - 1 ) uq 




. 


I'M 


/ • ‘ 


u^+Z^ 'uw+Z^I^I ' wl(v* + w^)^| J 




-I* fz, 

2 ^ 1 |w| 


u I w| + Z ’ 
‘ • ww 


|w (v^ + 1 + Z^^' v^# J 




|iV 






' 6b ] 




I'M 


wn' 


* (n* - 1) uw 


'(a*. l)w|(v^ . 




2 6sn* 


' (n- - 1) 63 






(F".) ... 


V 


Sin cot 






d* V 3 


U 






EW. cos 9 cos 


H Multiplied by 


- 



u 



for large angles of 
attack near -90^ 

Note 1 

v-'hen not multiplied by 
u add to Z ' 

D vp 
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Equation of Motion About the Body Axis System x-Axls 



-Iy)qr [k. ‘qr + K. ’/ + Kp(p, ,p,p 



I z* [k ' up + K ^ ' ur + K . ' V + K ^ ' wp] 



+ ^ Z^ [k ^ ' u* + K ^ ' uv + K ' vl(v* + w* )^l] 



* S K ' vw 
2 vw 



|z" u 






6r 



+ Bz g ein <+) cos 8 
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Equation of Motion About the Body Axis System y-Axls 



* 


[“^ ' 


q + M 


1'* 1 


[“, ' 


uq + M 


l‘‘ 1 


Fm . ' 

L w 


' w + M 


-f X* M , ' 
2 qn* 


(n> - 1) 


i‘“ 1 


[m, • 


+ M 


i‘‘ 1 


1 

L 1 w| 


' u j w j + 


[m 


6s 






Note 1 
I * r 



vp 



uw+M^l^l ' wl(v* 

' I w (v^ + 1 -f M ' 1 

W • • VV j 

6b ' ®b ] 



f ^ wn' ' 



+ f wlwln- ’ 






+ Bz ^ sin 0 



- EW ^ cos 8 cos ^ 



a Multiply by g* for 
large ancles ot attack 
near -90^ 

Note 1 

when not multipUed 
by ^ add to M ' 

U ^P 



m 



i:quaClon of Motion About the Body Axis System z-Axls 



^ ^ ‘V * V = I r ' ^ ^ ^ pq ' ^ ^ p ’ P] 

+ |l* [n^ ''^'■+^lrl6r ' N |v[r ' ^ 

■*■ I p ' “P ^ ^ ^ ^ wp ' '^] 



+ -|jl* N ' (n> - I) ur 



+ [n ^ + N ^ • uv + N ' v](v^ + v/^)^ l] 






+ ^ ^ , , ' (n* - 1) 6 

2 6rn' r 



+ -? X® N , ’ (n' - 1) uv 

2 vn* 



+ I-'" ‘'^vMn' - n vl(v^+wM‘ 



+ N • wv# 

2 WY 



+ Ew. cos 0 sia^^ 



VI 

S Multiply by ~ for large 
angles of attacK near -90'^ 
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AUXILIARY EQUATIONS 



e 

<> 

Xo 



yo 



io 



U 



p + sin0 

(q - Icosd sin({) ) / cos4> 

( r + 0 sin <t>) / cos 0 COs4^ 

UCO80CO3^J) + V {aini? sin 0 coslf) - cos({) sin 0 ) 
+ w (sin 4 ) sini^ + cos({) sin 0 cos i)) 

u cos 0 sin 0 + V ( cos4> cos + sin4) sin 9 sin <) ) 
+ w(cos<t)sin0sin0- sin 4) cos ^ ) 

- u sin 9 + V COS0 sin<{) + w cos 0 cos 4> 

( 



(F^)p = f 1 


a,' + ag' n' + as' ^ J 


when kj < n' 


= 


^ b^' + bj' n' + bj' J 


when kj < n’ < ki 


= fl-u> 


[ Cl’ + C 3 ’ n' + C 3 ' n'* ] 


I when k 3 < n’ < kg 


- 


[ di’ + dg' n' + dj' n'^ 


j when n‘ < k 3 



< I I 

i* I • I 

^li ^2f ^3 
• I I 

Cl, Cp, C3 

I 1 * 

^2* ^3 



Sets of non-d l^.ens iof.a I cocfriclcnts Msed In the 
pulsion equation above, set which, will be In e: feet 

at any tlrr.e durlnr, a slnulaled r.ar.euver vUL depend on the 
value of n‘ and the numbers 
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NOMENCLATUPJ^: 



All symbols used In the equations of motion and In the 
auxiliary e^^uacions and relationships which appear In this 
report are defined below. Any dimensions involved will be 
consistent with the f oo t -p otin d - sc c ond system of units. All 
angles are in degrees. The Fortran variables corresponding 
to these symbols are shown in Appendix B 



SYMBOL 



B 



m 



Si 

0 



u 



V 



V 



DEFINITION 

A dot over any symbol signi- 
fies differentiation with res- 
pect to time. 

Buoyancy force which Is posi- 
tive upwards . 

Mass of the submarine including 

the water In the free 

spaces. 

Overall length of the submarine 

Linear velocity of origin of 
body axes relative to an earth- 
fixed axis system. 

Component of U along Che body 
x-axls • 

Component of U along the body 

y-ax Is . 

% 

Component of U along the body 
t - a X 1 s 4 
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Command speed: A steady value 

of u for a given propeller rpm 
whenc^jj?) and control surface 
angles are zero. Sign changes 
with propeller reversal. 

Longitudinal axis of the body 
fixed coordinate axis system. 

Transverse axis of the body 
fixed coordinate axis system. 

Vertical axis of the body fixed 
coordinate axis system. 

Distance along the Xq axis of an 
earth-fixed axis system. 

Distance along the y^ axis of an 
earth-fixed axis system. 

Distance along the z*^ axis of an 
earth-fixed axis system. 

Component of angular velocity 
about the body fixed x-axls. 

Component of angular velocity 
about the body fixed y-axls. 

Component of angular velocity 
about the body fixed z-axls. 

The z^coordlnate of the center 
of buoyance (C3) of the subma- 
rine. 



Ang Ic of attack. 



o< 


Ang Ic of attack. 


Sb f Dfc 


Angle of drift. 

Deflection of bowplane ( or 
•allplane ) 


Sr ; Df 


Deflection of rudder. 


Ss j Dj 


Deflection of sternplane. 


1 

Q 

e 

$ 

4> 

C 


The ratio u^/u, 
c 

Pitch angle. 

Yaw angle. 

Rol 1 angle. 

Mass density of sea water. 




Weight of water blown from a 
particular ballast tank ident- 
ify! ed ty the integer assigned 
to the index i. 


6) 


An gu la r velocity. 


t 


Time. 


*tl 


Location along the body x-axis 
of the center of mass of th# 1^^ 

ballast tank when this tank is 
filled with sea water. 






I 



X 



Propulsion force (see auxi- 
liary equaclons and relation- 
ships) . 

Moment of Inertia of a sub- 
marine about the x-axls. 



I 



J 



Moment of Inertia of a sub- 
marine about the y-axls. 



I 



z 



Moment of inertia of a sub- 
marine about the z-axls. 



* 


^p* '^p|pf>^qr^> 


s' 


''r'* K*: 


Ky'. 


|v| ' ^ VW » 


M. * , M * , 


M ' , AH ’ , H 


q * rr ’ 


rp rp 1 


Ml 1 ^ M. ' 
1 w| q ’ u 


, H ' , M ' , M ' . 

vr vp qn 


M ' , M 1 , 


! M , 1 ’ , M ’ . M, ' 


w V 1 w 1 


1 W 1 ww vv 


Mxv ' 1 M , 

Ob wn 






5 s 



Non-dimensional constants each 
of which Is assigned to a parti- 
cular force terra In the equation 
of motion about the body x-axls. 



Non-dimensional constants each 
of which Is assigned to a parti- 
cular force term in Che equation 
of motion about the body y-axls. 
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«r’* "Mir’- 



V r 



”p * • ‘^wp * ‘^rn* ’ * * ’ * 



Non-dlnena lonal constants each 
of which Is assigned to a parti- 
cular force tera In the equation 
^v|v|*’ ^6r * ^6rn* * ^'vn* * ^v|v|n* * of motion about the body z-axli. 



N ’ 
wv 



X *,x *.x *,x,*,x *,x *. 

qq rr rp u ’ vr ’ wq * 

X * X * X ^ X * X * Non-dimensional constants each 
uu ’ vv * vw * 6 ri r * 5 s6 s * 



X * X * X * X * cular force term in the equation 

6 b6 b * vvn * wwn * ’ 6s5sn* * 



^<S z6 rn * 



of which is assigned to a parti- 
cular force term in the cquatlot 
of motion along the body x-axIs. 



r p pq p I p I y wp * 



"vlrl’’ ■'UUr'- %’• ''rn''. 



Y*'. 


> 

> 


Y » Y * 

’ 5r ’ ^6rn’ ’ 




V 1 n ' ’ 


Y * , (F ) 

UV ■/' V3 


^r’ 




Z. ’ , Z , Z ' 
w * V r vp 


AZ ' , Z 
^ vp ’ q 


’’ 'hi 


5s* 1 q 1 ’ 


Z 1 . z*' 

qn ’ * 


• "w’* 


wlw| |wl 


Z ’ , Z 
vw vv 


' Z ' 


7*7 * 

’ ^5b ’ vn* * 


^w| w|n* • 


^4sn' ’ 


(F ) 
z vs 



, Non-d imens 1 ona 1 constants each 
of which is assigned to a parti- 
cular force term in the equatlo' 
of motion along the body y-axis 



Non-dimensional constants each 
of which is assigned to a parti- 
cular forcp term In the equatio' 
of motion along the body z-axls 
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A 


Coefficient matrix of the state vector 


B 


Coefficient matrix of the input vector 


X 


State vector 


E 


Error vector 


u 


Input vector 


R 


Control effort weighting matrix 


C,C1 ,C2 


Control effort weighting parameters 


Cl ,CI1 ,CI2 


Control effort inverse weighting parame- 
ters 


Q 


Weighting matrix in quadratic performance 
index 


A 


Depth rate weighting factor 


B 


Pitch rate weighting factor 


D 


Pitch error weighting factor 


E 


Depth error weighting factor 


BG 


BZ/1 =B Z /I 


X 


Longitudinal CG of the submarine 


CE 


Characteristic equation 


SOPOC 


Stern only planes optimal controller 


SOPC 


Stern only planes controller 


BPOC 


Both planes optimal controller 


CMC 


Combined mode controller 


FSA , FSB , FSC FSD 
SCA.SCB.SCC ,SCD 


Feedback gains in 


FBA , FBB , FBC , FBD 
BCA.BCB.BCC .BCD 


Feedback gains in 


SOA,SOB,SOC ,SOD 
FA , FB , FC , FD 


Feedback gains in 


EOM 


Eqs of motion 



\85 



stern plane command in both planes submarine 

Fairwater plane command in both planes sub- 
marine 

Stem plane command in stern plane only sub- 
marine 
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Fig. I-l . Coordinates of waves 




Fig. 1-2. Near surface "venturi" effect 
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PLANT 




Fig. II-l. Pla.nt and optimal feedback 
controller . 
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Fig. Ill-la. Ramp force at AU 
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Fig. I II- lb. Depth vs. Tinie. Response to a ramp 

force at AU. Bounded controller with 
error limiters. Submarine not "in trim." 
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Fig. III-lc. Pitch vs. Time. Response to a ramp force 
at AU. Bounded controller with error 
limiters. Submarine not "in trim." 
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Fig. Ill-ld. Stern Plane Angle vs. Tiir.e. Response to 
a ramp force at AU. Bounded controller 
with error limiters. Submarine not 
"in trim. " 
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^ig. Ill-le. 



Fair\>?ater Plane Ancle vs. Time. P.esponse 
to a ramp force at AU. Bounded controller 
with error limiters. Submarine not "in trim." 
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Ordered 

Planes 



Figure III- 2 
Error Limiters 
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Fig. III-3a. Sinusoidal force at AU 







. Depth vs. Tine. Response to a sinusoidal 
force at AU. Bounded controller without 
noise filters. Submarine not "in trim." 






Fig. III-3b 
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rig. III-3C. Pitch vs. Time. Response to a sinusoidal 
force at AU. Bounded controller without 
noise filters. Submarine not "in trim." 
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Fig. III-3d. stern Plane Angle vs. Time. Response to 
a sinusoidal force at AU. Bounded con- 
troller without noise filters. Submarine 
not "in trim." 
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Fig. III-3e. Fairwater Plane Angle vs. Time. Response 
to a sinusoidal force at AU. Bounded 
controller without noise filters . Submarine 
not "in trim." 
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Fig. III-4a. Depth vs. Time. x^esponse to a sinusoidal 
input at AU. Unbounded controller with 
noise filters. Submarine not "in trim.” 
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Response to a sinusoidal 
input at AU. Unbounded controller with 
noise filters. Submarine not "in trim." 
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Fig. III-4b. Pitch vs. Time. 
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Fig. III-4c. Stern Plane Angle vs. Time. Response to a 

sinusoidal input at AU. Unbounded controller 
with noise filters. Submarine not "in trim." 
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Fig. III-4d. Fairwater Plane Angle vs. Time. Response 
to a sinusoidal input at AU. Unbounded 
controller with noise filters. Submarine 
not "in trim." 
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Depth vs. Time. Response to a depth order 
of 10 ft. Unbounded controller with error 
limiters and noise filters. 
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Fig. III-5a. 




X-SCALE= 40. s/inch 
Y-SCALE= S-XIO""* rad/ inch 



Pitch vs. Time. Response to a depth order of 
10 ft. Unbounded controller with error 
limiters and noise filters. 
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Fig. III-5c. Stern Plane Angle vs. Time. Response to a 
depth order of 10 ft. Unbounded controller 
with error limiters and noise filters. 
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Fig. III-5d. Fairwater Plane Angle vs. Time. Response 
to a depth order of 10 ft. Unbounded 
controller with error limiters and noise 
filters . 
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Depth vs. Time. Response to a depth order 
of 10 ft. Bounded controller without noise 
filters . 
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Fig. III-6a. 
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Fig, III-6b. Pitcli vs. Time. Response to a depth order 
of 10 ft. Bounded controller without noise 
filters . 
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Fig. III-6c. stern Plane Angle vs. Time. Response to 

a depth order of 10 ft. Bounded controller 
without noise filters. 
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Fairwater Plane Angle vs. Time. Response 
to a depth order of 10 ft. Bounded 
controller without noise filters. 
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Fig. III-6d. 
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Depth vs. Time. Submarine "in trim." 
Bounded controller without noise filters. 
Zero ordered depth. 
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Fig. III-7a. 
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Fig. III-7b. Pitch vs. Time. Submarine “in trim." 

Bounded controller without noise filters. 
Zero ordered depth. 
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Fig. III-7C. Stern Plane Angle vs. Time. Submarine 
"in trim." Bounded controller without 
noise filters. Zero ordered depth. 
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Fig. III-7d. Fairwater Plane Angle vs. Time. Submarine 
"in trim." Bounded controller without 
noise filters. Zero ordered depth. 
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Fig. III-8a. Depth vs. Time. Submarine not "in trim.’’ 
Response to depth and pitch orders. 
Tracking controller without error limiters 
noise filters. 
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Fig. IIl-8b. Pitch vs. Time. Submarine not "in trim." 

Response to depth and pitch orders. 
Tracking controller without error limiters 
and noise filters. 
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Fig. III-8c. Stern Plane Angle vs. Time. Submarine not 
"in trim." Response to depth and pitch 
orders. Tracking controller without error 
limiters and noise filters. 
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II-3d. Fairwater Plane Angle vs. Time. Submarine 
not "in trim." Response to depth and pitch 
orders. Tracking controller without error 
limiters and noise filters. 
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Depth vs. Tiraa. Submarine not "in trim." 
Response to depth and pitch orders. 
Sinusoidal input (AU). Tracking controller 
without error limiters and plane noise 
filters . 
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Pitch vs. Time. Submarine not "in trim." 
Response to depth and pitch orders. 
Sinusoidal input (AU). Tracking controller 
without error limiters and plane noise 
filters . 
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Fig. III-9c. Stern Plane Angle vs. Time. Submarine not 
"in trim." Response to depth and pitch 
orders. Sinusoidal input (AU). Tracking 
controller without error limiters and plane 
noise filters. 
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-9d. Fairwater Plane Angle vs. Time. Submarine 
not "in trim." Response to depth and pitch 
orders. Sinusoidal input (AU). Tracking con 
troller without error limiters and plane 
noise filters. 
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Depth vs. Time. Submarine "in trim 
Response to a depth order of 50 ft. 
Tracking controller (final) without 
error limiters and noise filters. 






Fig. Ill-lOa. 
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Fig. Ill-lOb. Pitch vs. Time. Submarine "in trim. 

Response to a depth order of 50 ft. 
Tracking controller (final) without 
error limiters and noise filters. 






n 

uO 




I 



I 



0 G. 



.H 

(de<5) 



• .<,1Q . 

’ 3^1HC^ 

" S ^ T r 



Stern Plane Angle vs. Time. Submarine 
"in trim." Response to a depth order of 
50 ft. Tracking controller (final) with- 
out error limiters and noise filters. 






Fig. III-lOc. 
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Fig. Ill-lOd. Fairwater Plane Angle vs. Time. Response 

to a depth order of 50 ft. Tracking con- 
troller (final) without error limiters and 
noise filters. 
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(a) Linear regulator problem 
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( bilinear tracking problem 



Fig. III-ll. Block diagrams for optimal feedback 
controllers . 
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Excises 
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Fig. IV-1 . Block diagram for stern planes only controller 
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Fig. IV-2. Block diagram for Both planes controller 
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Fjg. IV-3. Block diagram for the combined controller 
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Fig IV- 5b. Depth vs. Time. Response to a pulse fore 
at FT. Bounded controller without error 
limiters. (B=800, C-10, E=l) 
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Pitch vs. Time. Response to a pulse force 
at FT. Bounded controller without error 
limiters. (B=800, C=10, E=l) 
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Fig. IV-5c. 
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stern Plane Angle vs. Time. Response to 
a pulse force at FT. Bounded controller 
without limiters. (B=800, C=10, E=l) 
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Fig. IV-5d 
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Fig. lV-5e. Fairwater Plane Angle vs. Time. Response to 

a pulse force of FT. Bounded controller without 
error limiters. (B=800, C=10, E=l) 
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Fig, IV-5f. Depth vs. Time. Response to a pulse force at 

FT. Bounded controller without error limiters. 
(B=164, C=0.001, E=0.001) 
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Fig. IV-5g. Pitch vs. Time. Response to a pulse force at 

FT. Bounded controller without error limiters. 
(B=164, C=0.001, S=0.001) 
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stern Plane Angle vs. Time. Response to a 
pulse force at FT. Bounded controller with- 
out error limiters. (B=164, C=0.001, E=0.001) 
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Fig. IV-5h. 
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Fairv/ater Plane Angle vs. Time. Response to 
a pulse force at FT. Bounded controller with- 
out error limiters. (B=164, C=0.001, E=0.001) 






Fig. IV-5i. 
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Fig. IV-6a. Depth vs. Time. Response to a pulse force at 
FT. Bounded controller with error limiters. 
(B=300, C=10, E=l) 
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Fig. IV-6b. Pitch vs. Time. Response to a pulse force at 
FT. Bounded controller with error limiters. 
(B=800, C=10, E=l) 
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. Stern Plane Angle vs. Tinae. Response to a 
pulse force at FT. Bounded controller with 
error limiters. (B=800, C=10, E=l) 
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Fig. IV-6d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. Bounded controller with 
error limiters. (B=800, C=10, E=l) 







I 

i 



J 

* "i t 

Fig. 




Ki'J 







iV“5e. Depth vs. Time. Response to a pulse force at 
FT. Bounded controller with error limiters. 
(C=10, D=3000, E = l) 
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^^*3" 6f . Pitch vs . Tiins. Rssponss to a pulss fotcs 

at FT, Bounded controller with error limiters 
(C=10, D=3000, E=l) 
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Stern Plane Angle vs. Tine. Response to a 
pulse force at FT. Bounded controller ';ith 
error limiters. (C=10, D=3000, E=l) 






Fig. IV-6g. 
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IV-6h. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. Bounded controller with 
error limiters. (C=10, D=3000, E=l) 
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Fig. IV-7a. Step force at AU or P'T 
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Fig, IV-7b. Depth vs. Time. Response to a step force at AU. 

CMC uses optimal design in combined controllers 
with B=800, C=10, E=l. Parameter X=16 
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IV-7c. Pitch vs. Time. Response to a step force at AU. 

CMC uses optimal design in combined controllers 
with B=800, C=10, E=l. Parameter X=16 
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Fig. IV-7d. Stern Plane Angle vs. Time. Response to a step 
force at AU. CMC uses optimal design in com- 
bined controllers with B=800, C=10, E=l. 
Parameter X=16 
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PairT'^ater Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design in 
combined controllers with B=800, C=10, E=l. 
Parameter X=16 
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Fig. IV-8a. Depth vs. Time. Response to a step force at .VJ . 

CMC uses optimal design in combined controllers 
with B=800, C=10, E=l. Parameter X=8 
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Fig. IV-8b. Pitch vs. Time. Response to a step force at AU. 

CMC uses optimal design in combined controllers 
with B=800, C=10, E=l. Parameter X=8 
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Fig. IV-8c. Stern Plane Angle vs. Time. Response to 
force at AU. CMC uses opti.mal design in 
bined controllers with B=800, C=10, E=l. 
Parameter X=8 



a step 
com- 



Z5T 






r 



t I 




T 



^5CRie 



,<1Q 



UN IT 5 ' iHUH 



'4 •• ijU 




Fig. IV-8d. Faiin/\rater Plane Angle vs. Time. Response to = 
step force at AU. CMC uses optimal design in 
combined controllers with B=800 , C=10, S=l. 
Parameter X=8 
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Fig. IV-9a. Depth vs. Time. Response to a step force at AU. 

CMC uses optimal design in combined controllers 
with B=800, C=10, E=l- Parameter X=4 
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Fig. IV-9b. Pitch vs. Time. Response to a step force at AU. 

CMC uses optimal design in combined controllers 
with B=800, C=10 , E=l. Parameter X=4 
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Fig. IV-9c. Stern Plane Angle vs. Time. Response to a step 

force at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=4 
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Fig. IV-9d. Fair^vater Plane Aiigle vs. Time. Response to a 
step force at AU. CMC uses optimal design in 
combined controllers with B=800, C=10, E=l. 
Parameter X=4 
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Ga. Depth vs. Time. Response to a step 
CMC uses optimal design in combined 
with B=800, C=10, E=l. Parameter 
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Stern Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design in 
combined controllers with B=800, C=10, E=l. 
Parameter X=2 . 5 
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Fig. IV-lOc. 
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Fig. IV-lOd. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=2.5 



SS6 



LI . 



r 




:: cs^ j 

Fig. IV-lla. Depth vs. Time. Response to a step force 

at AU. CMC uses optimal design in combined 
controllers with B=8C0/ C=10, E=l. 

Parameter X=1.0 






Fig. IV-lli 
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Pitch vs. Tiras. Response to a step force 
at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=1.0 
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Fig. IV-llc. 


Stern 


Plane Angle vs. Time. 


step 


force at AU. CMC uses i 




in combined controllers with 




E=l. 


Parameter X=1.0 
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Fig. IV-lld. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10 
E=l. Parameter X=1.0 
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Fig. IV-12a. Depth vs. Time. Response to a step force 

at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.9 
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Fig. IV- 12b. Pitch vs. Time. Response to a step force at 
AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.9 
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Fig. IV-iac. 



Stern Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0.9 
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Fig. IV-12d. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0 . 9 
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Fig. IV-13a. Depth vs. Time. Response to a step force 

at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.7 
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Fig. IV- 13b. Pitch vs. Time. Response to a step force 

at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0 . 7 
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Fig. IV- 13c. 3 tern Plane Angle vs. Time. Response to- a 

step force at AU. CMC uses optimal design in 
combined controllers with B=800, C=10, E=l. 
Parameter X=0 . 7 
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Fig. IV~13d. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0.7 
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Fig. IV-14a. Depth vs. Time. Response to a step force 

at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.5 









I 



< 

• T, 









Fig. IV- 14b. Pitch vs. Tirae. Response to c step force 

at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.5 
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Fig. IV- 14c. Stern Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0.5 
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Fig, IV-14d. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=800 , C=10, 
E=l. Parameter X=0.5 
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Fig. IV- 15a. Depth vs. Time. Response to a step force at 
at AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0 . 3 
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Fig. IV- 15b. Pitch vs. Time. Response to a step force at 
AU. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.3 
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Fig. IV-15c. 
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Stern Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0 . 3 
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Fig. IV-15d. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=800 / C— 10, 
E=l. Parameter X=0.3 



as 6 



(n 



i r\j 



L ^ 



I' 

> I 






I 






i I 



; .> 
— fv 



■ ! 
I ’ 



!! 1 ' il 



» i I 



) I 



i I 



! I 



I I 



I i 

i ! 



1 ' i 1 

■; . i 

'! !i 

i, ji 
■ : ! ' 



t 1 



I I 



1 



1 ! 



! 



! I 



/ f 






rc,, 



■ u 






. CC (»■) 



<i0 



^ ' -ill 

- r , * i 1- 

3/ - 



Fig. IV-16a. 



Depth vs. Time. Response to a step force at 
AU. CMC uses optimal design in combined 
controllers with E=164, C=0.001, E=0.001. 
Param.eter X=6 



28T 



CP C.P 



n 

c' 

•>!» 



o 

LC 



o 




Pitch vs. Tir".3. Response to a step force at 
AU. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=6 
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Fig. IV- 16c. Stern Plane Angle vs. Time. Response to a step 
force at AU. CMC uses optimal design in com- 
bined controllers with B=164, C=0.001, E=0.001. 
Parameter X=6 
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Fig. IV-16d. Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=164, C=0.001, 
E=0.001. Parameter X=6 
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Depth vs. Time. Response to a step force 
at AU. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=3.5 
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Fig. IV-17b. Pitch vs. Time. Response to a step force at 
AU. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=3.5 
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Stern Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design 
in combined controllers with B=164, C=0.001, 
E=0.001. Parameter X=3.5 






Fig. IV-17C. 
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Fig. IV-17d. Fairi'/ater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=164, C=0.001, 
E=0.001. Parameter X=3.5 







Fig. IV- 18a. Depth vs. Time. Response to a step force 

at AU. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=2. 
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Pitch vs. Time. Response to a step force at 
AU. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=2 






Fig. IV- 18b. 
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stern Plane Angle vs. Time. Response to a 
step force at AU. CMC uses optimal design 
in combined controllers with B=164, C=0.001, 
E=0.001. Parameter X=2 
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Fairwater Plane Angle vs. Time. Response to 
a step force at AU. CMC uses optimal design 
in combined controllers with B=164, C=0.001, 
E=0.001. Parameter X=2 
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Fig. IV- 19a. Depth vs. Time. Response to a step input at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=1.0 
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Fig. IV- 19b. Pitch vs. Time. Response to a step input at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l, 

Parameter X=1.0 
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3tern Plane Angle vs. Time. P.esncnse to 
a step input at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=1.0 
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IV-19d. Fairwater Plane Angle vs. Time. Response 
to a step input at FT. CMC uses optimal 
design in combined controllers with 3=800, 
C=10, E=l. Parameter X=1.0 
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Depth vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0 . 9 
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Pitch vs. Time. Response to 
FT. CMC uses optimal design 
controllers with B=800, C=10 
Parameter X=0.9 
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Fig. IV- 20b. 
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Fig. IV-“20c. Stern Plane Angle vs. Time. Response to 

a step force at FT. CMC uses optimal design 
in combined controllers with B = 800, C=10 , 
E=l. Parameter X=0 . 9 
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Fig. IV-20d. Fairwatar Plana Angle vs. Time. Response to 
a step force at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0.9 
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Depth vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.8 
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Fig. IV-21a. 
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Fig. IV-21b. Pitch vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0 . 8 
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iig. IV-21c. StsiTx Plane Angle vs. Time. Response to a 
step force at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0.8 
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Fairv.’atar Plane Angle vs Tine 
a step force at FT. CMC uses 
in combined controllers with 
E=l. Parameter X=0 . 8 
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Fig. IV-22a. Depth vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l, 

Parameter X=0.7 
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Pitch vs. Tiir.e . liespor.se to a stat force 
at FT. CMC uses optimal design in combined 
controllers with B=800, C-10, E=l. 

Parameter X=0.7 
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Stern Plane Angle V3 . -rir-.a. Response to a 
step force at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0 . 7 
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Fig. IV-22d. Fairwater Plane Angle vs. Time. Response to 
a step force at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 
E=l. Parameter X=0 . 7 
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Fig. IV-23a. Depth vs. Time. Response to a step force 

at FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.05 
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Fig. IV-2 3b. Pitch vs. Tirris. Response to a step force 

at FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Parameter X=0.05 
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Fig. IV-23C. Stern Plane Angle vs. Time. Response to a step 
force at FT. CMC uses optimal design in com- 
bined controllers with B=800, C=10, E=l. 
Parameter X=0.05 
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Fairwater Plane Angle vs. Time. Response 
to a step force at FT. CMC uses optimal 
design in combined controllers with B=800, 
C=10, E=l. Parameter X=0.05 
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-"ig. IV-24a. Depth vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=0.9 
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rig. IV-24h. Pitch V 3 . Tirf.e. Response to a step force 

at FT. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=0.9 
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seem Plane Angle vs. Tire. Aesvensa to 
a step force at FT. CAC uses optimal design 
in combined controllers with B=154, C=0.001, 
E=0.001. Parameter X=0.9 
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Fairwater Plane Angle vs. Time 
a step force at FT. CMC uses 
in combined controllers with I 
E=0,001. Parameter X=0.9 
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Fig. IV-24d. 
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Fig. IV-25a. Depth vs. Time. Response to a step force at FT. 

CMC uses optimal design in combined controllers 
with B=164, C=0.001, C=0.001. Parameter X=0.5 
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Pitch V3. Time. Response to a step force at 
?T. CMC uses optimal design in combined 
controllers with B=164, C=0.001, E=0.001. 
Parameter X=0.5 
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Fig. IV-25b 
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Fig. IV-25C. Stern Plane Angle vs. Time. Response to a 
step force at FT. CMC uses optimal design 
in combined controllers with B=164, C=0.001, 
E=0.001. Parameter X=0 . 5 



325 



Fairwater Plane Angle vs. Time 
a step force at FT. CMC uses 
in combined controllers with B 
E=0.001. Parameter X=0 . 5 
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Fig. IV-2Ga. Depth vs. Time. Response to a step force 

at FT. CMC uses optimal design in combined 
controllers, with B=164, C=0.001, E=0.001. 

P ar ame ter X= 0 . 4 
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Fig. lV-26b. Pitch vs. Time. Response to a step toroe 

at FT. CMC uses optimal design in combined 
controllers, with B=164, C=0.001, E=0.001. 
Parameter X=0.4 



^^8 




r* :*C £ ,r ,d<^t“VX 

6*P4 '. ■ . .li»i*5«f iMtff ■ - 5-1 a "’■■■' .T1 i% 

• i Vif.. - _ C»0 < i 41 ' ' « w * ■;, w ^ ' X : ? SCO ' 




1 

I J 




1 

I 




J 







Stern Plane Angle vs. Time. Response uo 
a step force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0 . 4 
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Fig. IV-26c. 
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g. IV-26d. Fair\^?ater Plane Angle vs. Time. Response to 
a step force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.4 
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IV-27a. Depth vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers, with B=164, C=0.001, E=0.001. 
Parameter X=0 . 3 
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Pitch vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers, with B=164, C=0.001, 2=0. GJl. 
Parameter X=0 . 3 
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Stern Plane Angle vs. Time. Response to a 
step force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.3 



Fig. IV-27C 
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Fairwater Plane Angle vs. Time. Response to 
a step force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.3 




Fig. IV-27d. 
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Fig. IV-28a. Depth vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers, with B=164, C=0.001, E=0.001. 
Parameter X=0.05 
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Pitch vs. Time. Response to a step force at 
FT. CMC uses optimal design in combined 
controllers, with B=164, C=0.001, E=0.001. 
Parameter X=0 . 0 5 
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Fig. IV-28b. 
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Stern Plane Angle vs. Time. Response to a 
step force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.05 






Fig. IV-28C. 
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IV-28d. Fair'/vater Plane Angle vs. Time. Response to 
a step force at FT. CMC uses optimal design 
in combined controllers, with B=164, C^O.OOl, 
E=0.001. Parameter X=0.05 







Fig. IV-29a. Depth vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.7 
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Pitch vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0 . 7 
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Fig. IV-29b. 
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stern Plane Angle vs. Time. Response to a 
pulse fotce at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0 . 7 
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Fig. IV-29C. 
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Fig. IV-29d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.7 
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rig. IV- 30a. Depth vs. Tine. Desgonse to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, S=l. 
Parameter X=0 . 6 
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Fig. IV-30b. Pitch vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.6 




Fig. IV-30C. Stern Plane Angle vs. Time. Response to 

a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0 . 6 
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Fig. IV-30d. Fairwater Plane Angle vs. Tir.e. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 

E=l. Parameter X=0.6 
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Fig. IV- 31a. Depth vs. Time. Response to a pulse force 
at'”FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.55 
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Fig. iv-31b. Pitch vs. Tine. Response to a pulse force at 
FT. CMC uses optimal design in combined 
controllers, with B=800, C=10 , E=l, 

Parameter X=0.55 
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Fig. IV- 31c. Stern Plane Angle vs. Tine. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.55 
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Fig. IV-31d. Fair'/vater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers with B=800, C=10, 

E=l. Parameter X=0.55 
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Fig. IV-32a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined 
controllers with B=800, C=10, E=l. 

Param.eter X=0 . 5 
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Fig. IV-32b. Pitch vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.5 
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stern Plane Angle vs. Time. Response co a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.5 
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Fig. IV-32d. Fair^.'/ater Plane Angle vs. Tirr.e. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0 . 5 
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. Depth vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.45 
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Fig. IV-33b. Pitch vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, ^vith B=800, C=10, E=l. 
Parameter X=0.45 
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Fig. IV-33c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.45 



btT 






o 

'T 

'0 



i 

i 

o 

l:. 

i 




- 

Fair^.'/ater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800 , C=10 , 
E=l. Parameter X=0.45 
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cu') 

34a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in comJDined 
controllers, with B=800, C=10, E=l. 

Parameter X=0.3 
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34b. Pitch vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, S=l. 
Parameter X=0 . 3 
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Fig. IV-34c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0 . 3 
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Depth vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.05 
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Fig. IV-35b. Pitch vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.05 
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Fig. IV-35c, Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.05 
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Fig. IV-35d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.05 
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Fig. IV- 36a. Depth vs. Time. Response to a pulse force 
at FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 
Parameter X=0.005 
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Fig. IV-36b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined 
controllers, with B=800, C=10, E=l. 

Parameter X=0.005 
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Fig. IV-36c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.005 
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Fig. IV-36d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=800, C=10, 
E=l. Parameter X=0.005 
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Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.5 
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Fig. IV-37b. 


Pitch vs. Time. 



Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.5 
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Fig. IV-37C. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.5 
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Fig. IV-37d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0 . 5 
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Fig. IV- 38a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0 . 3 
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Fig. IV- 38b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.3 
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Fig. IV-38c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0 . 3 
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Fig. IV-38d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.3 
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Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.2 
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Fig. IV-39b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, vrith B=164, C=0.001, E=0.001. 
Parameter X=0 . 2 
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Fig. IV-39C. Stern Plane Angle vs. Time. Response to a 

pulse force at FT. CMC uses optimal design in 
combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.2 
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Fig. IV-39d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.2 
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Fig. IV-40a. Depth vs. Time. Response to a pulse force ar 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0 . 15 
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Fig. IV-40b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.15 



3SA 



lO 




Fig. IV-40C. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0 . 15 
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Fig, IV-40d. Fairv^ater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0,001, 
E=0,001. Parameter X=0.15 
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Fig, IV-41a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0 . 1 
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Fig. IV-41b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0 . 1 
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Fig. IV-4lc. Stern Plane Angle vs. -Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.1 
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Fig. IV-41d. Fairwater Plane Angle vs. Time. Response to 
a pulse force of FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0 . 1 
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Fig. iv-42a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.005 
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Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with B=164, C=0.001, E=0.001. 
Parameter X=0.005 
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Fig. IV-42c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.005 
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Fig. IV-42d. Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers, with B=164, C=0.001, 
E=0.001. Parameter X=0.005 
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Fig. IV-43a. 



Depth vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10, S=l, and 
SOPC with the same closed loop C.E. 

Parameter X=0.7 
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Fig. IV-43b. Pitch vs. Time. Response to a pulse 
FT. CMC uses BPOC with B=800, C=10, 
SOPC with the same closed loop C.E. 
Parameter X=0.7 
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Fig. IV-43c. Stern Plane Angle vs. Time. Response to a 

pulse force at FT. CMC uses BPOC with B=800, 
C=10, E=l, and SOP C with the same closed loop 
C.E. Parameter X=0.7 
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Fig. iv-43d. Fairwater Plane Angle vs. Time. Response to 

a pulse force at FT. CMC uses BPOC with B=800, 
C=10, E=l, and SOPC with the same closed loop 
C.E. Parameter X=0 . 7 
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Depth vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10, S=l, and 
SORC with the same closed loop C.E. 

Parameter X=0.5 



Fig. IV-44a. 
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Fig. IV-44b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10, E=l, 
and SOPC with the same closed loop C.E. 
Parameter X=0.5 
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Fig. IV-44C. Stern Plane Angle vs. Time. Response to a 

pulse force at FT. CMC uses BPOC with B=800, 
C=10, E=l, and SOPC with the same closed loop 
C.E. Parameter X=0.5 
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Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses BPOC with 
B=800, C=10, E=l, and SOPC with the same 
closed loop C.E. Parameter X=0 . 5 
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Fig. IV-45a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10, E=l, 
and SOPC with the same closed loop C.E. 
Parameter X=0 . 4 
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Fig. IV-45b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10, E=l, 
and SOPC with the same closed loop C.E. 
Parameter X=0.4 
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* Fig. IV-45C. Stern Plane Angle vs. Time. Response to a 

pulse firce at FT. CMC uses BPOC with B=800, 
C=10, E=l, and SOPC with the same closed loop 
C.E. Parameter X=0.4 
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Fig. IV-45d. Fairwater Plane Angle vs. Time. P.esponse to 
a pulse force at FT. CMC used BPOC with 
B=800, C=10, E=l, and SOPC with the same 
closed loop C.E. Parameter X=0.4 
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Depth vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10, E=l, 
and SOPC with the same closed loop C.E. 
Parameter X=0.3 
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Fig. iv-46b. Pitch vs. Time. Response to a pulse 
FT. CMC uses BPOC with B=800, C=10, 
and SOPC with the same closed loop C. 
Parameter X=0 . 3 
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Fig. IV-46c, Stern Plane Angle vs. Time. Response to a 

pulse force at FT. CMC uses BPOC with B=800, 
C=10, E=l, and SOPC with the same closed loop 
C.E. Parameter X=0.3 
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Fairwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses BPOC with 
B=800, C=10, E=l, and SOPC with the same 
closed loop C.E. Parameter X=0 . 3 



Fig. IV-46d. 
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Depth vs. Time. Response to a pulse force at 
FT. CMC uses BPOC with B=800, C=10 , E=l, and 
SOPC with the same closed loop C.E. 

Parameter X=0 . 1 
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Fig. IV 47a. 
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Fig. IV-ATo. Pitch vs. Ti:ne. Response to a pulse force at 
FT, CMC uses BPOC with B=800, C=10, E=l, 
and SOPC with the same closed loop C.E. 
Parameter X— 0.1 
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Fig. IV-47C. Stern Plane Angle vs. Time. Response to a pulse 
force at FT. CMC uses BPOC with B=800, C=10, 
E=l, and SOPC with the same closed loop C.E. 
Parameter X=0.1 
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Fairwater Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses BPOC with B=800, 
C=10, E=l, and SOPC with the same closed loop 
C.E. Parameter X=0 . 1 






Fig. IV-47d. 
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Fig. IV-48a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers with C=10, D=3000, E=l. 

Parameter X=0.4 
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Fig. lV48b. Pitch vs. Time. Response to a' pulse force at 
FT. CMC uses optimal design in combined con- 
trollers with C=10, D=3000, E=l. 

Parameter X=0 . 4 
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Fig. IV-48c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers with C=10, D=3000, 
E=l. Parameter X=0.4 
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Fig. IV-49a. Depth vs. Time. Response to 
FT. CMC uses optimal design 
trollers, with C=10, D=3000, 
Parameter X=0 .5 
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Fig. IV— 49b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with C=10, D=3000, E=l. 

Parameter X=0.5 
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in combined controllers, with C=10, D=3000, 
E=l. Parameter X=0.5 




Fig. IV-49C. 
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Faiirwater Plane Angle vs. Time. Response to 
a pulse force at FT. CMC uses optimal design 
in combined controllers with C=10, D=3000, 
E=l. Parameter X=0.5 
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Fig. IV-50a. Depth vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with C=10, D=3000, E=l. 

Parameter X= 0,6 
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Fig. IV-50b. Pitch vs. Time. Response to a pulse force at 
FT. CMC uses optimal design in combined con- 
trollers, with C=10, D=3000, E=l. 

Parameter X=0.6 






c/1 (/■) 




.,1 

J 

i ^ 1 I 

-j.:c -i.jo -j.o: le.. oj 2 c. jc 

CRLf;-40.00(s'^ UHIT5/IHCH 
CflLE- 3.00 (d€cj') UHITS^IHCH 

Stern Plane Angle vs. Time. Response to a 
pulse force at FT. CMC uses optimal design 
in combined controllers, with C=10, D=3000, 
E=l. Parameter X=0.6 



Fig. IV-SOc. 



n 

ct 




L 1 



-..jc o.oa "c ^ c.jc 

:kic > ^ 

KSCRLE-40-00 (s') UH!T5^IHCH 
VSCRLE-g.OO (de^l UHITS/IHCH 

^^9* IV~50cl. Fa.i 2 rwa.t 01 r Plana Ancla vs. Tiine . RGsponsG to 
a pulSG forcG at FT. CMC usgs optimal dasign 
in combinad controllars, with C=10, D=3000, 
E=l. Paramatar X=0.6 



42G 




Fig, V-la. Root locus for SOPOC with B=800, C=10, E=l. 
M > 0 
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Fig. V-lb. Root locus for SOPOC with B=800, C-10, E— 1. 
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Fig. V-2a. Depth vs. Time. Response to a pulse force at 
FT. SOPC equal to 0 . 3 • SOPOC (B=800, C=10, 
E=l) 
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Fig. V-2b. Pitch vs. Time. Response to a pulse force at 
FT. SOPC equal to 0.3-SOPOC (B=800, C=10, 
E=l) 
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Fig. V-2c. Stern Plane Angle vs. Time. Response to a pulse 
force at FT. SOPC equal to 0.3-SOPOC (B=800, 
C=10, E=l) 
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Fig. V-3a. Depth vs. Time. Response to a pulse force at 
FT. SOPC equal to 0 . 6 • SOPOC (B=800, C=10, 
E=l) 
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-3b. Pitch vs. Time. Response to a pulse force at 
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Fig. V-3c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. SOPC equal to 0.6"SOPOC 
(B=800, C=10, E=l) 
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Fig. v-4a. Depth vs. Time. Response to a, pulse force at 
FT. SOPC equal to 0.8 • SOPOC *(B=800 ,C=10 , 
E=l) 
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Fig. V-4b, pitch vs. Time. Response to a pulse force at 
FT. SOPC equal to 0.8-SOPOC (B=800, C=10 , 
E=l) 







Fig. V-4c. Stem Plane Angle vs. Time. Response to a pulse 
force at FT. SOPC equal to 0.8*S0P0C (B=800, 
C=10, E=l) 







Fig. V-5a. Root locus for SOPC v/ith CE same as for BPOC 
(B=800, C=10, E=l) M>0 
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Fig. V-5b. Root locus for SOPC with CE same as for BPOC 
(B=800, C=10, E=l) M<0 
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Fig. V-6b. Root locus for SOPOC with B=164, C=0.001, 
E=0. 001. M <0 







Fig. V-7b. Root locus for SOPC with CE same as for BPOC 
(B=164, C=0.001, E=0.001). M<0 
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Fig. V-8b. Root locus for SOPOC with D=30C0, C=10, E-1. 
M <0 
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Fig. V-9a. Root locus for SOPC with CE same 
(D=3000, C=10, E=l) . M>0 



as for BPOC 




Fig. V-9b. Root locus for SOPC with CE same as for 3?0C 
(D=3000, C=10, E=l) . M <0 
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Fig. V-lOa. 



Root locus for SOPC using Ssi 
(D=3000, C=10, E=l).Mi')0 
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Fig, V-lOb, Root locus for SOPC using Ssi of BPOC 
(D=3000, C=10, E=l) .Mi.<0 
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Fig. V-lOe. Root locus for SOPC using 6si of BPOC 
(B=164, C=0.001, E=0.001) .M^) 0 
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Depth vs. Time. Response to a pulse force at 
FT. Both planes submarine with oscoTr^-^si-fC^'^ • 
X=0. 5 
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Fig. V-llb. - Pitch vs. Time. Response to a pulse force at 
FT. Both planes submarine with oscov\ = bs^ • -fCH") • 
X=0. 5 
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Fig. V-llc. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. Both planes submarine 
with 6scorr\=,^Si--f( '<') • X=0 . 5 
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Fig. V-lld. Fairwater Plane Angle vs. Tine. Response to a 
pulse force at FT. Both planes submarine 
with O'oCO’t^^osi. • X=0 . 5 
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Depth vs. Time. Response to a pulse force at 
FT. Both planes submarine with ot)COTii = 65\ • 

X=-0 . 5 
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Pitch vs. Time. Response to -a pulse force 
FT. Both planes submarine with 6scorn = 6si • 
X=-0 . 5 
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Fig. V-12c. Stern Plane Angle vs. Time. Response to a pulse 
force at FT. Both planes submarine with 
S5CO'rq= • X=-0.5 
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a pulss fores at FT. Both planss sv±>marins 
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Fig. V-13b, Pitch vs. Time. Response to a pulse force at 

FT. Both planes submarine with a5COTn= 6 sl -fCO • 
X=-0 . 75 
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Fig. V-13d. Fairwater Plane Angie vs. Time. Response to 
a pulse force at FT. Both planes submarine 
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Fig. V-14a. Depth vs. Time. Response to a pulse force at 
FT. Both planes submarine with oiCo^^ ^osi • f ( k ) 
X=- 1 . 0 
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Fig. V-14b, 



Pitch vs. Time. Response to a pulse force 
at FT. Both planes submarine with 
^scOTT^= i • X=-1.0 
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Fig. V-14c. Stern Plane Angle vs. Time. Response to a 
pulse force at FT. Both planes submarine 
with <^acoT<i - 6si • f CK) • X=-1.0 
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Fig. V-14d. Fair^>/ater Plane Angle vs. Time. Response to 
a pulse force at FT. Both planes submarine 
with S^co-rg-'^si •fC'O • X=-1.0 
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Fig. V-15b. Pitch vs. Time. CMC 
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Fig. V-15c. Stern Plane Angle 
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Fig. V-15d. Fairwater Plane Angle vs. Time 
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Pig. B-la. Relation Between Body Axes and Fixed Axes 




Figure B-lb. Submarine Axes 




Figure B-lc. 
Axes Definitions 
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Figure B-ld . 
Axes Definitions 
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Fig. C-la. Depth vs. Time. Stern planes only submarine, 
not "in trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD=360) and 
zero pitch and depth ordered. 
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Fig. C-lb. Pitch vs. Time. Stern planes only submarine, 
not "in trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD=360) and 
zero pitch and depth ordered. 
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Fig. C-lc. Stern Plane Angle vs. Time. Stern planes only 
submarine, not "in trim" with feedback 
controller (SOA=-0.5, SOB=-7.03, SOC=32.5, 
SOD=360) and zero pitch and depth ordered. 
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Fig. C-2a. Depth vs. Time. Stern planes only submarine, 
not "in trim" with feedback controller 
(SOA=-0. 0436, SOB=-3.49, SOC=0.0523, 

SOD=360) and zero pitch and depth ordered. 



A90 




^SCRlE;^.- 2. 002-jHUN ITS/ INCH 

Fig. C-2b. Pitch vs. Time. Stern planes only submarine, 
not "in trim" with feedback controller 
(SOA=-0.0436, SOB=-3.49, SOC=0.0523, 

SOC=360) and zero pitch and depth ordered. 
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Stern Plane Angle vs. Time. Stern planes 
submarine, not "in trim" with feedback 
controller (SOA=-0 . 0436 , SOB=— 3.49, SOC=0 
SOD=360) and zero pitch and depth ordered 
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Fig. C-2c. 
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Depth vs. Tine. Stern planes only submarine, 
not "in trim" with feedback controller 
(SOA=-0.0656, SOB=-1.47, SOC=26.2, SOD=184.1) 
and zero pitch and depth ordered. 
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Fig. C-3b. Pitch vs. Time. Stern planes only submarine, 
not "in trim" with feedback controller 
(SOA=-0.0656, SOB=-1.47, SOC=26.2, SOD=184.1) 
and zero pitch and depth ordered. 
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Stern Plane Angle vs. Time. Stern planes only 
submarine, not "in trim" with feedback controller 
(SOA =-0 . 0656 , SOB=-1.47, SOC=26.2, SOD=184.1) 
and zero pitch and depth ordered. 
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Fig. C-3c. 
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Fig. C-4a. Depth vs. Time. Stern planes only submarine, 
"in -trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD=360) 
and zero pitch and depth ordered. 
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Fig. C-4b. Pitch vs. Time. Stern planes only submarine, 
"in trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD=360) 
and zero pitch and depth ordered. 







stern Plane Angle vs. Time. Stern planes only- 
submarine, "in trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD=360) and 
zero pitch and depth ordered. 
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Depth vs. Time. . Stern planes only submarine, 
"in trim" with feedback controller 
(SOA=-0.‘0436, SOB=-3.49, SOC=0.0523, 

SOD=360) and zero pitch and depth ordered. 
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Fig. C-5b. Pitch vs. Time. Stern planes only submarine, 
"in trim" with feedback controller 
(SOA=-0.0436, SOB=-3.49, SOC=0.0523, 

SOD=360) and zero pitch and depth ordered. 
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stern Plane Angle vs. Time. Stern planes only 
submarine, "in trim" with feedback controller 
(SOA=-0.0436, SOB=-3.49, SOC=0.0523, 

SOD=360) and zero pitch and depth ordered. 
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Fig. C-5c. 
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Depth vs. Time. Stern planes only submarine, 
"in trim" with feedback controller 
(SOA=-0 .0656, SOB=-1.47, SOC=26.2, 

SOD=184.1) and zero pitch and depth ordered. 
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Fig. C-6a. 
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Fig. C-6b. Pitch vs. Time. Stern planes only submarine, 
"in trim" with feedback controller 
(SOA=-0 .0656, SOB=-1.47, SOC=26.2, 

SOD=184.1) and zero pitch and depth ordered. 
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Fig. C-6c. Stern Plane Angle vs. Time. Stern planes only 
submarine, "in trim" with feedback controller 
(SOA=-0.0656, SOB=-1.47, SOC=26.2, 

SOD=184.1) and zero pitch and depth ordered. 
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Fig. C-7a. Depth vs. Time. Stern planes only submarine, 
not "in trim" with SOPOC (B=164, C=0.001, 
E=0.001) and zero pitch and depth ordered. 
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Fig. C-7b. Pitch vs. Time. Stern planes only submarine, 
not "in trim" with SOPOC (B=164, C=0.001, 
E=0.001) and zero pitch and depth ordered. 
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VSCflLE=1.00 UHITS/IHCH' 

Stern Plane Angle vs. Time. Stern planes only 
submarine, not "in trim" with SOPOC (B=164, 
C=0.001, E=0.001) and zero pitch and depth 
ordered. 






Fig. C-7c. 
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Fig. C-Sa. Depth vs. Time. Stern planes only submarine, 
"in trim" with SOPOC (B=164, C=0.001, E=0.001) 
and zero pitch and depth ordered. 
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Fig. C-Cb. Pitch vs. Time. Stern planes only submarine, 
"in trim" with SOPOC (B=164, C=0.001, E=0.001) 
and zero pitch and depth ordered. 
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Fig. C-8c. Stern Plane Angle vs. Time. Stern planes only 
submarine, "in trim" with SOPOC (B=164, 
C=0.001, E=0.001) and zero pitch and depth 
ordered . 
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Fig. C-9a. Sinusoidal force at AU 
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Fig. C-9b. Depth vs. Time. Response to sinusoidal force 
at AU. Stern planes only submarine, "in trim" 
with SOPOC (B=164,' C=0.001, E=0.001) and zero 
pitch and depth ordered. 
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Pitch vs. Time. Response to sinusoidal force 
at AU. Stern planes only submarine, "in trim" 
with SOPOC (B=164, C=0.001, E=0.001) and zero 
pitch and depth ordered. 
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Fig. C-9c. 
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Fig. C-9d. Stern Plane Angle vs. Time. Response to 

sinusoidal force at AU. Stern planes only 
submarine, "in trim" with SOPOC (B=164, 
C=0.001, E=0.001) and zero pitch and depth 
ordered. 
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Fig. C-lOa. Depth vs. Time. Response to sinusoidal force 
at AU. Stern planes only submarine, "in trim" 
with feedback controller (SOA=— 0.5, SOB=— 7.03, 
SOC=32.5, SOD=360) and zero pitch and depth 
ordered . 
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Fig. C-lOb. Pitch vs. Time. Response to sinusoidal force' 
at AU. Stern planes only siobmarine, "in trim" 
with feedback controller (SOA=-0.5, SOB=-7.03, 
SOC=32.5, SOD=360) and zero pitch and depth 
ordered. 
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Fig. C-lOc. Stem Plane Angle vs. Time. Response to 

sinusoidal force at AU. Stern planes only sub- 
marine, "in trim" with feedback controller 
(SOA=-0.5, SOB=-7.03, SOC=32.5, SOD=360) and 
zero pitch and depth ordered. 



607 




S5CRLF-80.0C (s') 
'.F5CRLE-C.05 ■ 0+') 



UNIiS/INCH 
UN ITS/ INCH 



Depth vs. Time. Response to sinusoidal force at 
AU. Stern planes only submarine, "in trim" 
with feedback controller (SOA=-0 . 0436, SOB=-3.49, 
SOC=0.0523, SOD=360) and zero pitch and depth 
ordered. 
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Fig. C-lla. 
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Fig. C-llb. Pitch vs. Time. Response to sinusoidal force 
at AU. Stern planes only submarine, "in trim" 
with feedback controller (SOA=-0 . 0436 , 
SOB=-3.49, SOC=0.0523, SOD=360) and zero pitch 
and depth ordered. 
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Fig. C-llc. Stern Plane Angle vs. Time. Response to 

sinusoidal force at AU. Stern planes only sub- 
marine; "in trim" with feedback controller 
(SOA=-0.0436, SOB=-3.49, SOC=0.0523, 

SOD=360) and zero pitch and depth ordered. 
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Depth vs. Time. Response to sinusoidal force at 
AU. Stern planes only submarine, "in trim" 
with feedback controller (SOA=-0 . 0656 , 

SOB=-1.47, SOC=26.2, SOD=184.1) and zero pitch 
and depth ordered. 






Fig. C-12a. 
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Fig. C-12b. Pitch vs. Time. Response to sinusoidal force at 
AU. Stern planes only submarine, "in trim" 
with feedback controller (SOA=-0 . 0656 , 

SOB=-1.47, SOC=26.2, SOD=184.1) and zero pitch 
and depth ordered. 
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Fig. C-I2c. Stern Plane Angle vs. Time. Response to 

sinusoidal force at AU. Stern planes only sub- 
marine', "in trim" with feedback controller 
(SOA=-0 .0656, SOB=-1.47, SOC=26.2, SOD=184.1) 
and zero pitch and depth ordered. 
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Fig. C-13a. Depth vs. Time. Response to sinusoidal force 

at AU. Stern planes only submarine, 

"in trim" with SOPOC (B=800, C=10, E=l) and 
zero pitch and depth ordered. 
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Fig. C-13b. Pitch vs. Time. Response to sinusoidal input at 
AU. Stern planes only _s.ub marine , "in trim" 
with SOPOC (B=800, C=10, E=l) and zero pitch 
and depth ordered. 
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Fig. C-13c. Stern Plane Angle vs. Time. Response to 

sinusoidal input at AU. Stern planes only 
submarine, "in trim" with SOPOC (B=800, 

C=10, E=l) and zero pitch and depth ordered. 
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Fig, D-2. X- vs. Speed 
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Fig. E-1. K vs. dfcom. and X=3.2 
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Fig. E-2. K vs. dfcom and X=1.2 
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Fig. E-3. K vs. dfcom and X=0 .T. 
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X-SCFILE--=2.00E-01 UNITS INCH. 
Y-SCRLE=2.00E-0i UNITS INCH. 



K vs . df com and X=0 . 5 






Fig. E-4. 
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Fig. E-5. K vs. dfcorr. and X=C.3 
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Fig. E-6. K vs. dfcom and X=0.05 
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